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Research  has  been  carried  out  in  the  areas  of  (l)  Laser  Spectroscopic 
Collisional  Studies,  (2)  Optical  Collisions,  (3)  Radiative  Collisions, 

(1»)  Laser  Assisted  Cooling  or  Heating  and  (5)  Photon  Echoes  and  Optical 
Ramsey  Fringes.  These  subjects  are  described  below.  In  cases  where 
preprints  or  reprints  are  not  yet  available,  more  detailed  accounts  are 
given. 

1.  Laser  Spectroscopic  Collisional  Studies  (P.  Berman,  J.  LeGouet) 

Three-level  systems  offer  an  interesting  and  novel  means  for 

studying  collisions  in  atomic  and  molecular  systems.  Our  work  in  this 

1 2* 

area  has  been  summarized  in  two  recent  review  articles  * with  the  latter 
of  these  articles  stressing  the  physical  principles  involved  in  three-level 
systems.  Current  work  involves  the  extension  of  these  calculations  to  allow 
for  different  collision  models,  to  include  the  effects  of  strong  radiation 
fields,  and  to  test  the  validity  of  generally  used  collision  kernels. 

The  three-level  system  under  consideration  is  shown  below 


A - 

A - S).  -u! 


* Stars  indicate  articles  published  or  submitted  during  the  report  period 
which  are  appended  to  this  report.  The  reader  is  referred  to  these  articles 
for  detailed  discussions  of  the  subject  matter. 
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(0  < a < 1 ; u = most  probable  speed  of  thermal  distribution)  giving  the 
probability  density  per  unit  time  for  a collision  induced  velocity  change 
v'  -*■  v is  valid  for  weak  collisions,  but  is  suspect  for  strong  or  hard 
sphere  collisions  if  the  perturber  to  active  atom  mass  ratio  is  large  and 

if  v*  ;>  u.  It  is  a kernel  that  has  been  commonly  used  owing  to  its 

analytical  properties,  but  has  received  limited  experimental  verification. 

Several  situations  in  which  the  above  approximations  became  invalid 
have  been  explored  and  are  discussed  below. 

Collisions  which  are  velocity-changing  in  their  effect  on  optical 

j— 

coherences  . When  the  collision  interaction  is  nearly  the  same  for  two 
levels  involved  in  a transition,  collisions  are  velocity-changing  (rather 
than  phase-interrupting)  in  their  effect  on  the  optical  coherence  associated 
with  those  levels.  We  have  carried  out  a theoretical  calculation  for  this 
case.  In  particular,  we  have  shown  that,  if  the  populations  in  levels  1 

and  2 are  nearly  the  same,  owing  to  some  incoherent  pumping,  then  the 

probe  absorption  line  shape  clearly  reflects  whether  collisions  are 
phase-interrupting  or  velocity-changing  in  their  effect  on  the  1-3  optical 
coherence.  An  experimental  test^  of  these  results  for  excited  state  rare 

-U- 


gas  systems  seems  to  be  consistent  with  phase-interrupting  collisions 
rather  than  velocity  changing  ones. 

Transit-time  effects.  In  analyzing  experimental  data^  in 
which  one  of  the  energy  levels  is  metastable,  the  importance  of  transit 
time  effects  became  evident.  In  this  case,  the  lifetime  of  the  atom  is 
determined  by  the  time  it  spends  in  the  laser  field  rather  than  by  the 
decay  time  of  the  levels  involved.  We  were  able  to  show  that  a simplified 
model  of  transit  time  effects  is  in  good  agreement  with  more  extensive 
calculations^  carried  out  in  the  weak  field  limit.  In  this  model,  the 
metastable  decay  rate  y is  replaced  by  y + v^/R  where  v^  is  a transverse 
velocity  and  R is  the  radius  of  the  light  beam.  The  transit  time  can 
significantly  affect  the  saturation  parameter  (important  in  strong  field 
calculations)  which  is  proportional  to  y~  and  this  simple  model  allows 
for  an  adequate  treatment  of  this  effect.  Inclusion  of  transit  time 
effects  led  to  a correction  to  the  value  of  the  saturation  papameter 
measured  by  linear  absorption  techniques  on  excited  states  in  the  rare 
gases. ^ 

Strong  pumr>  field.  In  many  cases  of  experimental  interest,  the 

pump  field  exciting  the  1-2  transition  is  strong  enough  to  invalidate  a 

perturbation  theory  approach.  Thus  it  is  of  important  practical  interest 

to  generalize  our  earlier  work  to  allow  for  strong  ptnnp  fields.  In  the 

7 

absence  of  collisions,  line  shape  formulas  have  been  derived  which, 
although  algebraically  complicated,  are  easily  evaluated  on  computer. 

Line  splittings  due  to  the  ac  Stark  effect  may  be  seen  in  some  cases. 

We  have  obtained  the  equations  giving  the  probe  field  absorption  in 
a three-level  system  with  a strong  pump  field  and  collisions  occurring 
in  the  system.  The  probe  field  absorption  may  be  written  as  (for  simplicity 
level  1 only  is  assumed  to  be  incoherently  pumped  - «=  0) 


^ V 


where  r^c  » rate  of  velocity-changing  collisions  in  level  i;  x.X'  “*e 
field  strengths  E , E'  in  frequency  units;  w ij(v'-»v)  is  the  collision  kernel 
for  state  i;  Yjj  is  decay  rate  of  Pjj  (i  4 j)  including  effects  of 
phase-interrupting  collisions;  y £ = decay  rate  of  level  ? back  to  level 
1. 


These  equations  are  virtually  impossible  to  solve  for  arbitrary 

collision  kernels.  A key  quantity  is  the  saturation  broadened  homo- 

B B 

geneous  width  y . In  most  laser  spectroscopy  experiments,  Y^2  is 

less  than  the  Doppler  width  and  Y?0/k  represents  the  width  of  the  hole 

burnt  in  the  velocity  distribution  by  the  pump  field.  If  kAu  <^Y^ 

(Au  = rms  velocity  change  per  collision)  collisions  are  relatively  weak 

and  a good  approximation  is  to  take  the  collision  kernel  of  the  form 

W(v'-v).  For  this  kernel  we  have  solved  the  equations  and  obtained 

solutions  in  terms  of  Bessel  functions.  As  yet  we  have  not  numerically 

evaluated  the  line  profiles  that  were  obtained. 


Another  case  of  interest  which  we  have  solved  is  that  in  which 

kAu  > y12»  1 <e-  collisions  which  are  strong  enough  to  remove  atoms 

from  the  velocity  hole  created  by  the  pump  field.  In  that  limit,  the 

terms  in  the  eouations  may  be  dropped  and  the  equations  for 

12 

reduce  to  the  form  of  the  weak  field  equations  * ' with  the  replacement 

g 

of  Y^g  Tkese  equations  were  solved  for  the  Keilson-Storer 

Kernel  and  the  resulting  probe  field  absorption  is 

j cc 
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where  e - 1 (copropagating)  or  -1  (counter-propagating  waves),  a is  the 
parameter  in  the  Keilson-Storer  kernel,  and  2 is  the  Plasma  Dispersion 
function. 


In  the  weak  field  limit,  the  first  sura  reduces  to  the  two-quantum 
contribution  to  the  line  shape  and  the  second  to  the  step  wise  contribution 
from  atoms  which  have  not  experienced  velocity-changing  collisions.  The 
nth  term  in  the  third  sura  provides  the  contribution  from  atoms  undergoing 
n velocity-changing  collisions  in  level  2 before  being  excited  to  level  3. 
We  have  not  yet  analyzed  the  line  profile  in  detail  to  see  the  effects  of 
collisions  on  line  splittings,  but  have  so  far  restricted  our  analysis  to 
particular  experimental  situations  to  be  discussed  below. 

Experiment  of  Liao  and  Bjorkholm.  Liao  and  Bjorkholm  have 

carried  out  a saturation  spectroscopy  experiment  on  a three-level  system 

in  Na  perturbed  by  the  rare  gases.  Using  pump  detunings  of  *4.0  GHz  and 

1.6  GHz  (compared  with  a Doppler  FWHM  of  1.7  GHz),  they  obtained  probe 

field  absorption  profiles.  The  analysis  of  the  data  is  at  the  preliminary 

stage,  but  some  conclusions  have  been  reached.  Typical  large  ( 1* . 0 GHz) 

detuning  data  is  shown  on  the  next  page.  At  these  large  detunings,  there 

is  no  velocity-selection  and  velocity-changing  collisions  are  unimportant. 

The  collision  induced  absorption  about  A'  = 0 is  the  "collisional  re- 

1 2 

distribution  term"  that  has  been  described  elsewhere  ’ and  is  also  re- 
ferred to  as  excitation  via  an  optical  collision.  A fit  of  the  experiment 
to  theory  with  no  free  parameters  (all  broadening  parameters  are  available 
from  other  experiments)  yields  very  good  agreement. 

At  smaller  detunings  (1.6  GHz),  velocity-selection  and  velocity- 
changing collisions  play  an  important  role  as  can  be  seen  in  the  data  on 
the  following  page.  The  "shoulder"  in  the  data  arises  from  velocity- 
changing collisions  and  increases  in  magnitude  with  pressure.  First 
attempts  to  analyze  the  data  seem  to  indicate  that  the  Keilson-Storer 
kernel  does  not  adequately  explain  these  results.  We  shall  test  whether 
a hard-sphere  model  gives  better  agreement. 

Thus,  this  experiment  offers  great  premise.  With  large  detunings, 
one  measures  collision  induced  absorption  and  the  various  width  and 
shift  parameters.  With  smaller  detunings,  the  nature  of  the  collision 
kernel  can  be  studied. 
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2.  Optical  Collisions  (S.D.  Yeh,  P.R.  Berman) 


In  "optical  collisions",  a collision  is  used  to  provide  energy 
to  aid  in  radiative  excitation.  We  have  termed  such  a process  "Collisionally 
Aided  Radiative  Excitation  (or  Emission)"  (CARE)  and  have  carried  out 

numerical  calculations  of  CARE)  in  two  level  systems  for  attractive,  re- 

8* 

pulsive  and Lennard-J ones  type  potentials.  Saturation  effects,  dressed 

states,  bare  states,  and  collision  induced  level  crossings  are  discussed 

0 

in  our  article. 

In  addition,  a model  problem  was  solved  with  the  goal  of  obtaining 
additional  insight  into  the  nature  of  the  results  in  CARE  calculations.  For 
an  interatomic  potential  of  the  form 

V(t)  = A(b)  e‘ltl/T{b) 

(b  = impact  parameter),  an  analytic  solution  in  terms  of  tabulated 
functions  may  be  obtained.  This  solution,  valid  for  arbitrary  field 
strengths  gives  an  excitation  probability  P 
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where  7b  is  the  perturber  density,  v the  relative  speed,  and  <|>(a,b,c)  a 
confluent  hypergecmetric  function.  We  have  not  yet  analyzed  this 
solution  in  detail;  the  asymptotic  properties  of  the  <(>  function  are  not 
fully  developed  in  the  literature.  However,  certain  properties  of  the 
solution  are  apparent. 

The  asymptotic  solution  for  larege  "red"  detunings  (attractive 

potential  - negative  detuning)  where  collision  induced  level  crossings 

o 

occur  agrees  well  with  results  using  power  law  potentials.  For  large 

"blue"  detunings  the  exponential  and  power  law  results  differ  in  two 

ways.  First,  the  potential  we  have  chosen  has  a discontinuity  in  slope 

at  t = 0 and  this  gives  rise  to  an  additional  asymptotic  contribution 

not  found  with  smooth  potentials.  Even  if  this  term  is  discarded, 

however,  the  remaining  asymptotic  solution,  on  integration  over  impact 

parameter,  has  an  inverse  power  dependence  on  detuning  rather  than  the 

9 

exponential  one  found  by  other  authors.  Further  analysis  of  our  results 

9 

as  well  as  those  obtained  by  others  is  needed  before  the  usefulness  of 
this  model  potential  can  be  determined.  A similar  potential  was  used  by 
Crothers^0  for  radiative  collisions  and  may  possess  similar  drawbacks. 

Work  is  continuing  on  CARE  in  three-level  systems  with  special 
emphasis  being  given  to  understanding  interference  effects  in  these 
systems. 

3.  Radiative  Collisions  (E.J.  Robinson) 

Considerable  experimental  and  theoretical  attention  has  been  paid 

in  recent  years  to  the  problem  of  "radiative  collisions",  inelastic 

processes  in  which  two  dissimilar  atoms  collide  under  the  influence  of 
H 10 

a laser  field.  One  atom,  A is  initially  excited,  the  other,  B 

in  its  ground  state,  and  the  radiation  is  (nearly)  resonant  with  the 
difference  in  excitation  of  the  initial  state  of  A and  a level  of  B. 

As  a result  of  the  combined  effect  of  the  collision  and  the  light-matter 

interaction,  A makes  a transition  to  its  ground  state,  B to  an  excited 

*« 

state,  and  a photon  absorbed.  Dynamically,  A*+B  + hV-*A+B  . 


-13- 


The  customary  theoretical  formulation  which  we  followed^  treats  the 
radiation  field  and  nuclear  motion  classically  (with  the  further  assumption 
that  the  latter  consists  of  straight  line  paths  for  both  atoms)  and  the 
electrons  ouantally.  In  addition,  it  is  possible  to  reduce  the  problem  to 
an  effective  two  level  system,  where  initial  state  jl>  is  the  composite 
(A*B)  and  final  state  |2>  is  composite  (A,B**).  Assuming  that  the  atom- 
atom  potential  may  be  adequately  represented  by  its  dipole-dipole 

H« 

asymptotic  form  (see  the  references  and  the  enclosed  preprint  ) and  that 
the  rotating  wave  approximation  holds  the  effective  Hamiltonian  is  . 

Na  + £7j-)  + *L3  r"3  5 i><»  e + 

, where  are  the  van  der  Waals'  and  Stark  co- 

efficients for  state  |i>,  and  B an  angle-averaged  constant  sum  over  off- 
diagonal  matrix  elements.  A is  the  detuning.  Atomic  units  are  employed 
here  and  below. 

The  time-dependent  Schroedinger  equations  for  the  state  amplitudes 
and  a2  become  the  simultaneous  linear  equations 

N,-  $|,  ^-e4±+  ( R ^ ) S (i.) 

(lb) 

A useful  alternative  form  is  given  by  a transformation  to  the  inter- 
13 

action  representation. 

. 

' lr!3  e*p  ^ 1 } (2a) 

1 frL  ■£  J 

- ffs  /79-<P04w]  ■ (2b> 


Either  Eqs.  (l)  or  (2)  are  to  he  solved  subject  to  the  initial 

condition  at  t = - a.  = 1,  a.  * 0.  In  terms  of  the  impact  parameter 
2 2 2 1/  ^ ^ 

p,  R(t)  = (v  t + p )/*■  , where  v is  the  relative  velocity  of  A and  B. 

Our  initial  plan  was  to  exploit  the  simple  structure  of  this 
effective  two  level  problem  by  solving  it  numerically  on  a machine,  and 
to  present  the  results  without  being  concerned  about  analytic  solutions. 

We  found,  however,  that  the  step  size  in  time  required  for  large  detuning 
and/or  small  impact  parameters  tended  to  produce  extremely  time-consuming 
and  costly  computer  Jobs.  We  also  found  that  the  asymptotic  methods  dis- 

llj 

cussed  by  the  Soviet  authors  were  quantitatively  unreliable.  Consequently, 
we  elected  to  seek  model  problems  which  could  be  solved  in  closed  form. 

A useful,  simpler  version  of  the  problem  is  the  weak  field  limit. 

7 

(in  the  context  of  the  present  problem,  "weak"  means  less  than  about  10 
watt/cm^. ) Geltman*'’  and  Knight*^  have  obtained  analytic  solutions  to  this 
problem  in  the  low  intensity  case  by  ignoring  the  level-shifting  term  in 
the  effective  Hamiltonian,  V =(Q  -Q0)/R^.  This  simplification  enables 
the  problem  to  be  reduced  to  an  integral  representation  of  a modified 
Bessel  function  K . However,  it  yields  some  unphysical  results,  e.g. , 

lineshapes  which  are  symmetrical  about  exact  resonance.  It  also  fails 
to  generate  oscillations  of  transition  probabilities  for  impact  parameters 
smaller  than  the  Weisskopf  radius.  These  omissions  are  quite  serious,  and 
lead  us  to  prefer  calculations  in  which  is  taken  into  account,  at 
least  approximately. 

An  immediate  generalization  of  the  calculations  of  Geltman*^  and 
Knight‘d  may  be  made  if  one  replaces  the  actual  asymptotic  level  shifting 
term  by  a delta- function  in  time,  where  the  coefficient  of  6(t)  is 
adjusted  to  the  actual  phase  induced  by  the  R term.  The  R-0  transition- 
inducing  term  is  included  exactly  in  this  treatment.  We  have  completed 
our  analysis  of  this  model  and  submitted  it  for  publication.**  The 
closed  form  nature  of  the  solution  of  the  earlier  investigators  is 
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retained,  but  the  spurious  symmetry  of  the  spectrum  about  A = 0 is  removed. 
The  oscillations  referred  to  above,  vhich  arise  from  interference  effects, 

are  present,  also.  Curiously,  the  impact  model  predicts  a line  shape 

12  13 

closer  to  experiment  than  does  "exact"  numerical  calculations.  We 

are  unable  to  account  for  this  apparent  "superiority",  and  it  might  be 

due  to  a fortuitous  cancellation  of  errors  applicable  only  to  the  specific 

atoms  studied  by  the  Harris  group. 

It  would  be  of  interest  to  extend  the  impact  level  shift  model  to 
arbitrary  field  strengths.  We  were  unable  to  accomplish  this  without 
some  modification,  which  we  now  discuss. 

The  interaction  representation  form,  Eqs.  (2), of  the  two  level 


equations  is  of  the  general  form 

i (ft 

(3a) 

/o!L-V(-£)e  o (3b) 

The  pair  of  simultaneous  first  order  equations  may  be  replaced  by  two 
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non-coupled  second  order  equations 

w,  - L (£f  + i A]  '5,  + IVI'-,  = O)  "*«) 

+ 1V|\(2  = q d,b) 


For  arbitrary  fields,  it  is  necessary  to  solve  Eqs.  (3)  or  (U)  without 
the  approximation  = 1.  Even  with  VLg  = 0,  I have  been  unable  to  find 
an  analytic  solution  to  this  equation.  (Since  there  does  not  seem  to  be 
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any  reference  to  the  solution  to  this  equation  in  the  literature, 

apparently  neither  has  anybody  else).  If  Eqs.  (*♦)  were  solvable  for 

VTO  = 0,  it  would  be  possible  to  obtain  an  immediate  generalization  to 

the  case  VTO  = A6(t),  a strong  field  impact  level  shift  model. 

Lo 

A problem  that  has  been  solved  exactly  with  V = 0 is  the  Rosen- 

Lo 

Zener  problem,  for  which  V = K sech  irt/T,  instead  of  the  present 

(v  t + p ) ' . (Presumably,  this  arises  from  the  milder  nature  of  the 

singularity  in  the  complex  t plane  for  the  former  compared  to  the  latter 

- - simple  pole  vs.  branch  point).  Along  the  real  t axis,  however,  there 

—3 

is  quite  a close  resemblance  between  R (t)  and  K sech  irt/T  - - that  is, 
there  should  be,  at  most,  only  a small  change  in  the  numbers  generated 
by  a power  law  coupling  potential  and  a hyperbolic  secant  coupling 
potential.  We  verified  this  in  the  weak- field  impact  level  shift  case, 
choosing  the  parameters  of  K sech  irt/T  so  that  the  two  potentials  are 
equal  at  t = 0,  and  gave  the  same  transition  amplitude  at  zero  detuning. 
(Since  we  work  in  the  time  domain,  the  constants  characterizing  the 
hyperbolic  secant  vary  with  the  impact  parameter).  It  turns  out  that 
not  only  does  the  hyperbolic  secant- impact  level  shift (HSILS) model  agree 

_3 

exactly  with  the  R -impact  level  shift  model  at  zero  detunings,  (this 
is  by  construction)  but  also  for  large  detunings,  with  the  greatest  dis- 
crepancy between  the  two  occurring  at  intermediate  values  of  A,  where  we 
found  cross-sections  differing  by  up  to  about  10/5.  We  consider  this 
situation  to  be  satisfactory. 


For  completeness,  the  weak-field  transition  amplitude  in  the  HSILS 
model  is  given  by  the  quadrature,  | or.  ( -►  «°)  | = 2K|cos  ^ I sech  — 

J o 

cosAt  dt  + sin  ~ I sech  ~ sinAt  dt|,  where  (h  is  defined  as  for  the  R~^ 
coupling  potential  (see  preprint),  and  the  integrals  given  by  j 

cosax  sech  Bx  dx  « ir/2B  sech  ~ } C sin  U*  S «£x.V|  Bx  cL  = - 


where  ip  is  the  digamma  function. 


-17- 


£ 


l 
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Thus  we  regard  the  substitution  of  the  hyperbolic  secant  for  R as 
a minor  additional  approximation.  Our  motivation  for  the  replacement  is 
that  the  HSILS  model  may  be  directly  extended  to  the  strong  field  limit, 
which,  for  t J 0,  obeys  the  same  differential  equations  as  apply  to  the 

pn 

classical  problem  of  Rosen  and  Zener.  Equations  (3)  become,  for  non- 
zero times 


iA't 

e 


e 


(5a) 


(5b) 


Making  the  change  of  variable,  z = (tanhirt/T  + l)/2,  Eqs.  (*<)  become^, 
for  t 4 0 (z  J 1/2) 


- aktf I = o, 

VT  1 At  1 At 

where  a = - b = — , = — - i — , C2  = — + i — . Primes  denote  differ- 

ention  with  respect  to  z.  This  is  the  hypergeometric  equation.  Eqs.  (6) 
are  subject  to  the  indicated  initial  conditions  at  t = a (z  = 0),  so 
that  for  z < 1/2,  the  solutions  are  identical  to  the  pure  Rosen-Zener 
problem.  At  z = 1/2,  the  level-shifting  delta  function  is  switched  on. 
Taking  the  delta  function  to  act  only  in  the  upper  state  (this  may  be  done 
without  loss  of  generality  since  only  the  relative  level  shift  is  of 
physical  significance),  it  changes  the  phase  of  a2  by  <f>,  while  leaving  the 
magnitude  unchanged.  The  amplitude  factor  a^»  and  hence  is  unchanged 
by  the  level  shifting  potential.  Therefore  the  latter  is  continuous  at 
z = 1/2.  These  provide  new  initial  conditions  at  z = 1/2  for  the  solution 
of  Eqs.  (6)  for  the  region  z > 1/2.  Our  desired  solution  is  the  value  of 


(6a) 

(6b) 
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a2  at  z = 1 ( t = + °°) , to  wit 


IT  b esc  b 

where  AX  - ~'3  ^3-^^  AAl 


(7) 


Tr(&i.R.O-^^))-  HiQJ 

6 - aTa-tCt-g^)  kf>t  Q-i 

IT  ( &tPx  Cl  - e®))-  ’ 

/Jl=l)  &l  = FCc,>b,ei)  •t) 

A = J=I  . U*(rT  C V (a-e-z + »,  \rc^ 1 >x  vk) , 


G,=  FCaA/^-W),  H,  = FC^I,  a-c.,-k), 

^ _ ok  F(a-tl)  t+i^dj.+  lji)  ( fi={-b  p(a-<!+') 
Wb  a-c.,-k)+  0-c-^±)C^  F(a-c^ 


The  F's  are  Gauss  hypergecmetric  series. 
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The  manerlcal  work  for  the  strong  field  HSILS  model  has  been  com- 
pleted and  is  being  analyzed.  Our  findings  indicate  that  there  is  no 

perceptible  difference  between  weak  field  and  strong  field  formulations 

& 2 

up  to  an  intensity  of  about  2 x 10  watt/cm  for  the  Ca-Sr  system 
studied  experimentally  by  the  Harris  group.  Above  that  power  density, 
we  find  a shift  of  the  spectral  maximum  to  the  red,  and  a small  change 
in  the  lineshape.  Also,  as  the  intensity  rises,  the  relative  contri- 
bution of  small  impact  parameters  diminishes.  Consequently , the  shift 
in  the  peak  seems  to  be  almost  entirely  due  to  the  A.C.  Stark  effect. 

That  partof  the  shift  due  to  the  exact  treatment  of  the  Ksech  irt/T 
coupling,  which  is  in  the  opposite  direction  toward  the  blue,  seems  not 
to  have  much  effect,  since  it  is  greatest  at  small  imoact  parameters, 
whose  relative  contribution  to  the  total  cross  section  becomes  less  and 
less  at  high  intensities. 

A paper  based  on  the  HSILS  model  should  be  completed  and  submitted 
for  publication  in  the  near  future. 

Thus,  the  simulation  of  an  inverse  power  law  potential  by  a hyper- 
bolic secant  seems  to  be  well-motivated  and  useful.  On  the  other  hand, 
the  replacement  of  R ^(t)  by  a delta  function,  while  yielding  an  exactly 
solvable  model  seems  to  be  less  convincing,  even  though  it  generates 
spectral  properties  known  to  be  consistent  with  experiment,  and  is  a distinct 
improvement  over  calculations  which  omit  level  shifts  altogether.  Con- 
sequently, we  were  led  to  explore  a model  in  which  the  potentials  are 

pure  hyperbolic  secant,  i.e.,  R replaced  by  Ksech  irt/T,  and  R~  by 
2 

Psech  irt/T.  In  the  weak  field  limit,  this  means  that  we  need  the 
integral 

I +**0  — kj  cli:  .secJiTL*:  £ iPSsecVi 

= eAle*p£  , 


(8) 


In  terms  of  dimensionless  parameters,  this  is  proportional  to 


ft  y r 'i  i 

~r—  $ Sec-We  ef.p\ifl'^‘lnhx\cix-  (9) 

-Wj 

The  integral  in  Eq.  (9)  nay  be  formally  written  in  terms  of  a 
contour  integral.  Consider  the  path  shown  in  Fig.  1.  We  first  consider 
the  case  of  \i  = 0 for  simplicity. 


Figure  1 

For  large  R,  the  contributions  to  the  contour  integral  arising  from  the 
vertical  sections  vanishes  exponentially.  The  two  horizontal  sections 
are  easily  shown  to  be  proportional  to  one  another.  Thus,  apart  from  a 
multiplicative  constant,  the  integral  of  Eq.  (9)  is  given  by  the  residue 
theorem.  For  p = 0,  there  is  only  the  simple  pole  at  iir/2,  and  the 

integral  becomes  X e^*dx  =X=  TT . 

A similar  procedure  should  apply  for  p ¥ 0.  The  same  contour  is 
used,  the  vertical  segments  of  the  contour  integral  vanish  as  M“,  and 
the  horizontal  segments  are  proportional  to  one  another,  leaving  the 
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desired  integral  expressible  in  terms  of  the  singularity  at  iu/2.  This, 
however,  is  no  longer  a pole,  but  an  essential  singularity.  To  evaluate 
the  residue  in  a straight f orward  way , we  would  require  the  Laurent 
series  for  sechxe*^tanhX  and  the  Taylor  series  for  e*®X  in  the  neighborhood 
of  the  non-analytic  point,  which  ordinarily  would  involve,  among  other 
things^ expanding  e*lJ'tanhx  e^xj  and  sechx.  This  gives  a mess,  and  it 
turns  out  that  a simpler  procedure  is  available.  Writing  the 
integral  in  Eq.  (9)  in  a series  representation  at  the  outset,  we  have 


X = Z.  ch  x e 'icrn\\ 


(10) 


The  zeroth  term  has  already  been  given.  The  coefficient  of  U may^be 
calculated  via  a parts  integration  H > (Jk 


i 0X  I 


— — se oh  * 


I ^ I 

^xe  cl*  ~ i PX>. 


— 


' ^06 


For  k > 2,  we  have  again,  via  parts  integration, 
r?0  k , . $>  1&X  » r>  , . . , ' 

^ setUx.  -tnnn  vy;  € =-  o v^’Ptli x "raw  Vu)r*«ih  < *?  Ok- 

-s^cVix^i te/x  e^X  4i  3 ^s<?tiurfrMih,c  x eQ,tcix 

4^-l)  $ k €2^*  dx  ^ »hK 

4-  S'eJt  x "to  h V)  fc~9'x  <2?X  cIk  — (jcH  1 *'-Uk  h c X e clxj 


or>;rfc*  ; PXfc_v  +(^-0Xfe.t-<5i-1)Xsfel  or 

^IR-I  +-  , 


Thus,  instead  of  exolicitly  working  out  the  residue  at  the  essential 
singularity,  we  have  the  series 


arity,  we  nave  we  series  r. 
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where  the  terms  are  obtainable  via  the  recursion  formula 

fc>2. 


This  appears  to  be  a simple,  tractable  procedure  and  numerical  work  is 
being  carried  out. 

It  is  noted  that  a similar  trick  works  in  the  case  of  optical 
collisions.  The  weak  field  transition  probability,  in  the  dressed  atom 
basis,  requires  the  evaluation  of  the  integral 

^ s.  dt*  R ^ e \ 1 ^ ^ ^ ^ • 


We  may  construct  a hyperbolic  secant  model  analogous  to  the  radiative 
caseby  replacing  the  R (t)  term  by  a sech^nt/T  in  both  the  exponential 
and  its  prefactor,  reducing  the  problem  to  an  integral  of  the  form 

^ • ! ’ |k*  i/vtonh* 

'T  — S dx  Sec  h * e e 

We  use  the  same  contour  as  previously  (Fig.  (l)).  Again,  the  vertical 
segments  give  a vanishing  contribution  to  the  integral,  and  the  hori- 
zontal segments  give  contributions  proportional  to  each  other,  leaving 
the  integral  to  be  expressed  in  terms  of  the  singularity  at  iu/2.  For 
U = 0,  this  is  now  a second-order  nole.  For  u 4 0,  we  have  the  series 
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and  may  avoid  the  need  to  evaluate  the  residue  at  an  essential^ingularity 
by  again  formulating  the  problem  recursively.  Thus,  T=  X «£LVto 

,(ix  , t , ,JS  , . . 

"T"^  — Sec  X X o'/.  &ov-£v~>r  <V 

~o&  ° ^ 

fij.  I,  X = TT  P 'X,*  dX 

« S;  (svJi  * -httbliK)  s««.«  X e c‘x  — — s «cl,  * e | _ c* 

— c>&  V 

^ , iftrt  I , I iftx.  I 

IftS  &ckLxe  dx-  ^ SecL,  x~iWix  <2  clX)  or 
• 1 

rr,=  .'P>j^.-x„  x=  ,'££>  - * stir  csd,"^  t Fir-k^z, 

Sedi^X  ^xcix  ~ ^ (setlix  'hh^K  ) (S&ck  X 

^ '~c*  ^ n*  ^ 

-kwh*'1  \)  e . ^ ^ + ~ 

+(£-nS*  % 

— ca  ~°Ci 

4 S' s^JiVk^*' x , i^.%aS^X\-te<ili  * e 

■s*  . , . O 

. sa-l^/  -h?nl')c  xe  11  dx^($-i)^^JiS<  "Z^1  ^ 

^ .j.  i l^L  . (\  'sf^oec^x  -kuk^  * 

e ^ Vl'  V ' 

(Dr,  Xjfe^XiM.  +‘J^± 

In  contrasting  the  weak-field  integration  required  for  hyperbolic 
function  coupling  and  level  shifting  potentials  with  the  R-n  potentials 
occurring  in  the  actual  atom-atom  system,  we  note  that  in  the  former  case. 
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the  singularity  in  the  complex  plane  is  a pole  for  y = 0,  and  an 

_3 

essential  singularity  otherwise.  In  the  R case,  there  is  a branch 
point  for  zero  level  shift  while  when  level  shifts  are  present,  the 
integrand  decomposes  into  one  term  which  possesses  an  essential 
singularity  and  a second  which  consists  of  a sum  of  functions,  each  of 
which  possesses  a branch  point.  The  relatively  milder  pathology  for 
the  hyperbolic  secant  case  leads  one  to  suspect  that  it  is  possible  to 
successfully  formulate  it  in  the  strong  field  case  as  a generalization 
of  the  Rosen-Zener  problem.  This  proves  to  be  the  case. 

We  now  consider  the  strong  field  generalization  of  this  problem. 

In  the  notation  of  Eqs.  (3),  V(t)  = {Ksech  nt/T}exp(-iytanhTrt/T) , so 
that  Eqs.  (U)  become 


Using  the  Rosen-Zener  problem  as  a model,  we  make  the  change  of  variable 
z = (tanh  irt/T  + l)/?,  so  that 


20-kW-  L^-V)-2i/u  f (12a) 

^ Z~  lL|  o ifO.  (12b) 

If  we  cast  this  into  a form  analogous  to  the  hypergeometric  equation 

*L  I =0>7 
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\-~TzkiL 

p x^- 


1 — T 1 fc  z 

/OiTTt 


bz  ~ ~ ) 


di  = iyu,  dl  = -cl, , e^e^o. 


The  hypergeometric  equation  (HGE),  is,  by  comparison 

20-23  f"  +Dg+t+i'ia-cj  f'  4QbF=o, 

The  HGE  has  regular  points  at  z = 0,  1,  and  m.  Our  equation  retains 
the  regular  singular  points  at  1 and  0,  but  has  an  irregular  singular 
point  at  That  is,  the  general  solution  of  Eos.  (13)  will  have  an 
essential  singularity  at  z = ~.  Comparing  the  z and  t planes,  it  is 
obvious  that  the  essential  singularity  that  appeared  in  the  weak  field 
integrand  has  been  mapped  to  infinity  in  the  new  variable.  The  general 
solution  at  the  finite  singular  points  will  have  either  a pole  or  a 
branch  point,  however,  a situation  that  persists  even  for  p ^ 0. 

We  attempt  to  solve  Eqs.  (13)  by  the  power  series  method,  looking 
particularly  for  one  solution  that  goes  to  unity  at  z = 0 (t  = - «). 
Writing  (the  subscripts  1,  2 are  suppressed) 
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3=0 
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we  obtain  the  three  term  recursion  relation 


oh 


C.  } 


As+ 1 - As 


(ou-sXi'Vs) 


As -i 


eAid 

<$+c\lHO 


2 z.. 


This  gives  a solution  of  the  relevant  differential  equation,  analytic  at 

z = 0,  and,  since  the  nearest  singular  point  of  the  differential  equation 

is  at  z = 1,  must  also  be  convergent  throughout  the  circle  Jz|  < 1.  We 

require  the  second  solution,  looking  again  at  the  HGE  for  guidance.  We 

c 

make  the  transformation  a = z P.  The  equation  for  P is  again  of  the 
same  form,  to  wit, 

(Y'-f  [Ja+I  -hY*  ~c  + JV3  3 -r  (al+-e 

(Hi) 

where  — O.  2.-  2_c_  ) Q b n ob  4-  I )0  'O  *t*^(C~i)/ 

c = j-c)  e-e  + 3-(  \-od. 

Clearly,  an  analogous  series  for  P may  be  written,  with  the  same  re- 
cursion formulae.  Thus  a complete  solution  has  been  found  in  a convenient 
form.  We  seek  a second  generalization  tractable  at  z = 1. 


To  accomplish  this,  we  make  use  of  the  symmetry  inherent  in  the 
equation  between  the  singular  points  at  z » 0 and  1.  Making  the  change 
of  variable  u = 1 - z,  we  obtain  a differential  equation  in  the  variable 
u in  which  the  roles  the  two  finite  singular  points  are  interchanged 

y (o-OoC ", t < **£ * 0*. - c%  lO'w'Jo/' \ E af+  e Wj oi  - G, 

(15) 

where  primes  denote  differential  with  respect  to  u,  and 

UJj  -ab  -te.  Z ^ <x^lr4\ -c-Zd 

A->  a /vj 


The  two  solutions  in  the  new  representation  may  be  immediately  written 
down  from  our  original  pair.  These  are  convergent  (in  terms  of 

z)  in  the  range  0 < z < 1.  Since  one  of  the  new  solutions  vanishes 
at  z = 1,  and  the  other  goes  to  1,  all  we  require  for  the  transition 
amplitude  is  the  coefficient  of  the  latter.  We  thus  proceed  at  follows. 
Using  the  series  centered  around  z = Q,  and  the  initial  conditions,  we 
write  down  a solution.  For  the  point  z = 1/2,  we  write  also  the 
solution  in  terms  of  the  two  series  centered  around  z = 1.  The  two 
forms  of  the  solution  must  be  analytic  continuations  of  each  other,  and 
the  values  of  functions  and  first  derivatives  must  be  the  same  for  the 
two  forms  in  their  common  domain  of  analyticity.  We  use  this  information 
to  evaluate  the  coefficients  of  the  two  independent  solutions  evaluated 
in  series  form  relative  to  z = 1. 

Preparation  for  numerical  analysis  of  this  model  have  begun.  It 
also  seems  likely  that  an  analogous  differential  eouation  may  be  written 
for  the  optical  collision  case.  This  will  be  investigated. 
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4.  Laser  Assisted  Heating  or  Cooling  - (P.R.  Berman  with  S.  Stenholm). 

21* 

We  have  proposed  a plan  to  use  collisions  in  the  presence 

of  a radiation  field  to  either  take  up  or  provide  kinetic  energy  in  the 

absorptive  process  and,  consequently,  cool  or  heat  the  sample.  A recent 

22 

successful  experiment  using  ions  in  a trap  uses  a method  similar  in 
spirit  to  our  proposal,  although  it  is  the  trap  rather  than  the  collision 
partners  which  absorb  the  energy. 

5.  Optical  Ramsey  Fringes  and  Photon  Echoes  - (J.  Le  Gouet  and  P.R.  Berman). 

We  have  begun  a study  of  photon  echoes  in  a standing  wave 

laser  field.  It  turns  out  that,  using  two  strong  standing  wave  pulses 

separated  by  a time  t,  one  obtains  a true  photon  echo  at  t = 2t  plus 

additional  signals  at  t = nx  (n  = integer  > 2).  These  subsequent  signals 

are  equivalent  to  the  Optical  Ramsey  Fringes  recently  discussed  by 
23 

several  authors.  An  interesting  feature  of  the  strong  field  calculation 
is  that  population  pulsations  appear  at  exp(2inkvt)  for  atoms  moving 
with  velocity  v.  Thus,  small  changes  in  velocity  owing  to  velocity- 
changing collisions  can  be  magnified  by  a factor  n so  that  standing 
wave  echoes  may  prove  useful  for  studying  weak  velocity-changing 
collisions.  Our  theoretical  results  are  at  a preliminary  stage. 
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ABSTRACT 


An  analytic  expression  is  derived  for  the 
transition  probability  during  a radiative  collision 
for  a model  problem  in  which  induced  level  shifts 
are  treated  in  the  impulse  approximation.  The 
formulae  predict  a behavior  that  qualitatively 
resembles  that  obtained  for  realistic  potentials , 
and  is  in  fair  agreement  with  the  experimental 
lineshape  obtained  by  Harris  and  his  colleagues. 


1 


Figure  Captions 

Figure  1:  Energy  Level  diagram  of  composite  system. 

System  absorbs  a photon  and  undergoes  a 
transition  from  A2-B1  to  A1-B3. 


Figure  2:  Radiative  collision  spectrum  for  Ca-Sr. 

Cross  section  for  absorption  as  a function 
of  detuning. 

a;  Impact  approximation  calculation,  this 
work. 

3;  Smoothed  graph  based  on  experimental  curve 
of  Harris  and  White. 

X;  Numerical  calculation  of  Harris  and  White. 
(Note:  each  curve  is  normalized  to  the 
same  maximum.  The  ordinates,  obviously,  are 
in  arbitrary  units.) 
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Study  of  velocity-changing  collisions  in  excited  Kr  using  saturation  spectroscopy 
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A saturated-absorption  experiment  is  used  to  measure  the  effects  of  velocity-changing  collisions  between 
Kr*  metastable  atoms  (♦p*5j[3/2)j)  and  He  and  Ar  perturbers.  Comparison  between  experimental  profiles 
and  profiles  calculated  assuming  hard-sphere  collisions  between  Kr*  and  perturber  atoms  confirms  the 
calculations  of  Borenstein  and  Lamb  concerning  the  change  of  velocity  associated  with  a collision. 
Furthermore  it  is  found  that  the  rate  of  velocity-changing  collisions  is  consistent  with  the  predictions  of 
kinetic  theory. 


I.  INTRODUCTION 

We  report  the  first  systematic  saturation  spec- 
troscopy study  of  velocity-changing  collisions 
(VCC)  between  electronically  excited  active  atoms 
and  ground-state  perturbers.  By  choosing  a per- 
turber-active-atom  mass  ratio  less  than  unity, 
we  are  able  to  clearly  isolate  the  effects  of  close- 
encounter  collisions.  As  such,  our  analysis  re- 
presents an  initial  step  in  obtaining  specific  in- 
formation on  the  short-range  part  of  excited- 
state  interatomic  potentials  by  monitoring  VCC 
using  saturation  spectroscopy,  and  it  indicates 
the  potential  of  this  method  for  studying  scattering 
in  short-lived  excited  states. 

The  study  of  velocity-changing  collisions  (VCC) 
in  the  ground  and  excited  states  of  atoms  or  mol- 
ecules using  laser  spectroscopic  techni  ;ues  has 
received  increased  attention  over  the  p„st  few 
years.1  *•  In  a typical  experiment  one  uses  a 
narrow-band  laser  source  of  frequency  O tuned 
near  an  atomic-transition  frequency  at  to  selec- 
tively excite  atoms  with  an  axial  velocity  v, 
-{n-ai/b,  (*  = n/c);  only  atoms  with  this  axial 
velocity  will  see  a Doppler-shifted  laser  frequency 
in  resonance  with  their  own  transition  frequency. 

A probe  laser  beam  is  then  used  to  monitor  any 
velocity  change  undergone  by  these  excited  atoms. 
Although  only  one  velocity  component  is  studied 
by  this  technique,  information  on  the  nature  of  the 
collision  interaction  can  be  inferred  from  the  re- 
sults of  such  experiments.  Theoretical  aspects 
of  the  problem  are  examined  in  Refs.  10-  12. 

To  date,  experiments  have  been  analyzed  as- 
suming collisions  are  ol  a “weak”  or  "strong” 
nature.  Weak  collisions  involve  relatively  small 
velocity  changes  Am  (am  « u = width  of  the  equilib- 
rium velocity  distribution)  and  presumably 
characterize  collisions  with  large  impact  param- 


eters. A strong  collision  usually  refers  to  a 
collision  which,  on  average,  completely  therm- 
alizes  the  atom’s  velocity  distribution  (a m = m) 
as  might  occur  following  a close  collision  between 
a heavy  perturber  and  a light  active  atom.  In  a 
previous  publication,  we  reported  on  the  influence 
of  collisions  on  saturated  absorption  profiles  of 
the  557-nm  line  of  Kr  I Kr*-Kr  and  Kr*-Xe 
collisions  were  considered  and  interpreted  in  terms 
of  a strong-collision  model. 

In  this  work,  we  study  the  effect  of  close- 
encounter  collisions  between  Kr*  atoms  and  He 
or  Ar  perturbers,  for  which  both  the  strong-and 
weak -collision  models  are  inadequate.  Saturated- 


V.  100  MK, 

A*/ys 


FIG.  1.  Typical  recording  of  the  557-nm  line  of 
Kri  in  the  presence  of  He  and  Ar  perturbers.  Upper 
trace:  P*,,  8 mTorr;  Ph.,  260  mTorr.  Lower  trace: 
PKr,  8 mTorr;  P^,  220  mTorr.  The  horizontal  line 
represents  the  Doppler  width  /v"2.  In  this  experiment 
the  frequency  scale  is  determined  by  Fabry- Perot 
fringes  spaced  by  83  MHz. 
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absorption  line  shapes  for  the  557-nm  transition 
(4p55s[  |]2  - 4p55p'[  |],)  in  KrI  were  obtained  using 
an  experimental  arrangement  that  has  been  de- 
scribed elsewhere.®  The  pressure  of  krypton  was 
kept  at  a fixed  value  of  8 mTorr  and  the  pressures 
of  He  and  Ar  perturbers  were  varied  between  0 
and  700  and  0 and  250  mTorr,  respectively. 
Typical  recordings  of  saturated-absorption  pro- 
files are  shown  on  Fig.  1.  Saturated-absorption 
profiles  have  shapes  and  widths  varying  with  mass 
and  pressure  of  perturbers;  they  could  be  con- 
sistently fit  over  the  entire  pressure  ranges  as- 
suming they  are  the  sum  of  a narrow  resonance  of 
width  2y  [full  width  at  half  maximum  (FWHM)] 
and  a broad  background.  The  narrow  component 
corresponds  to  atoms  that  have  not  undergone 
VCC,  and  the  broad  one  corresponds  to  atoms 
that  have  experienced  VCC.  The  shape  of  the 
background  signals  provides  direct  information 
characterizing  the  strength  and  cross-section 
for  Kr*-Kr,  Kr*-He,  and  Kr*-Ar  collisions.  In 
this  paper  we  use  a hard-sphere  model  to  de- 
scribe the  collisions,  but  more  detailed  models 
could  be  used  in  the  future. 

II.  THEORETICAL  ASPECTS 

The  following  assumptions,  consistent  with  the 
experimental  results,  were  used  to  arrive  at  a 
theoretical  line  shape. 

(i)  The  narrow  resonance  in  the  line  shape,  as 
modified  by  phase-interrupting  collisions,  satura- 
tion broadening,  weak  VCC,“•7•9■1,  beam  misalign- 
ment, and  laser  jitter  can  still  be  well  approxi- 
mated by  a Lorentzian  with  some  appropriate  width 
2y. 

(ii)  The  VCC  redistribute  the  velocity  of  Kr*  in  the 
4/>,5sf|]2  metastable  state;  effects  of  VCC  on  the 
rapidly  decaying  4/>55 />'[}  }l  state  can  be  neglected. 


ku  r0(r0  + r,) 


(iii)  Close-encounter  Kr*-He  or  Kr*-Ar  collisions 
can  be  characterized  as  hard-sphere  collisions 
with  a collision  kernel  of  the  Keilson-Storer  type14 

!V(  v ' - v)  = r2  [ tr(  Am)*  ] 's 

xexp[-(v- av'f /(&!<?] , (1) 

where  W'(v'- v ) is  the  probability  density  per  unit 
time  to  have  a velocity  change  from  v'  to  v, 

A«  = (l-a  2)l/2u  (2) 

is  VT times  the  rms  velocity  change  per  collision, 
u is  the  most  probable  speed  of  the  equilibrium 
distribution,  or  is  a parameter  related  to  the 
strength  of  a collision,  and  ra  is  the  rate  of  VCC. 
For  hard-sphere  collisions,  Borenstein  and  Lamb15 
have  calculated  a as  a function  of  perturber-active- 
atom  mass  ratio  using  a computer  simulation  of 
the  collisions. 

(iv)  A single  VCC  is  sufficient  to  remove  atoms 
from  the  velocity  holes  created  by  the  field,  i.e., 

Am  > 2y/fe.  This  condition  is  satisfied  in  our  ex- 
periment. 

(v)  The  Kr*-Kr  collisions  can  be  approximated 
as  “strong”  collisions  leading  to  a thermal 
distribution  after  a single  average  collision. 
Hard-sphere  collisions  between  equal-mass 
atoms  or  excitation-exchange  collisions  are 
approximately  described  by  such  a collision 
model.15  The  data  of  Fig.  1 contain  both  the 
Kr*-Kr  and  Kr*-He  or  Ar  contributions  to  the 
broad  background. 

Assuming  that  the  field  strengths  are  weak 
enough  to  justify  a rate-equation  approximation 
to  the  atom- field  interaction  and  using  the  above 
collision  model,  one  arrives  at  the  saturated- 
absorption  line  shape10*11  for  equal  pump  and  probe 
frequencies, 


(3) 


where  A = 0 - u>  is  the  detuning,  r,  is  the  Kr*- 
Kr  VCC  rate,  r2  is  the  Kr*-He  or  Kr*-Ar  VCC 
rate, 

r =t 

1 0 1 f 

where  t is  the  lifetime  of  Kr*  atoms  in  the  laser 
beams  (experimentally,  since  the  natural  life- 
time of  the  metastable  is  very  long,  r is  de- 
termined by  the  transit  time  in  the  beams,  which 
will  be  pressure  dependent)  and 


r'  =r0+r,+r1. 

The  line  shape  (3)  has  a simple  interpretation. 
The  first  term  arises  from  atoms  undergoing  no 
VCC  during  the  time  r,  the  second  term  arises 
from  atoms  undergoing  at  least  one  Kr*-Kr 
collision  in  the  time  r leading  to  a complete 
thermalization  and  the  final  term  corresponds  to 
atoms  undergoing  no  Kr*-Kr  collisions  and  some 
Kr*-He  or  Kr*-Ar  collisions  in  the  time  t lead- 
ing to  a partial  thermalization  of  the  sample. 
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III.  ANALYSIS  OF  THE  DATA 

C-.tl' 

Experimental  profiles  were  compared  to  the- 
oretical line  shapes  corresponding  to  expression 
(3)  by  making  the  following  assumptions:  For 
the  fixed  Kr  pressure  of  8 mTorr  chosen  in  this 
experiment  we  have  taken  the  rate  Fj  for  ther- 
malizing  Kr*-Kr  collisions  as  r,  = 1.4  x 105/sec 
determined  from  our  previous  experiment.8  We 
have  used  a hard- sphere  model  for  the  determination 
of  a and  ra;  a was  chosen  from  the  results  of 
Borenstein  and  Lambls:  a =0.936  for  Kr«-He  and 
<*  =0.6  for  Kr*-Ar.  {A  larger  value  of  a implies 
a smaller  as  fcf.  Eq.  (2));  for  perturber-active 
atom  mass  ratios  »1,  a =0,  and  for  perturber- 
active-atom  mass  ratios  «1,  a»l.}  The  rate 
I\  was  chosen  from  the  kinetic  theory  of  gases 
with  the  "radii”  of  atoms  as  indicated  by  Allen, 16 
ra  = 1.0x10*1)  sec'1  for  Kr*-He  and  r2  = 0.53xl0*p 
sec"1  for  Kr*-Ar,  where  p is  the  perturber  pres- 
sure in  mTorr.  At  each  perturber  pressure,  the 


FIG.  3.  Variations  of  the  transit  time  t vs  pertur- 
ber pressures.  (Dots  correspond  to  He-Kr*  and 
crosses  to  Ar- Kr*  data.)  The  broken  lines  correspond 
to  the  variations  given  by  the  classical  law  of  diffusion. 
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FIG.  2.  Typical  fits  for  He  and  Ar  perturbers.  Upper 
trace:  P^,  8 mTorr;  pHt,  110  mTorr.  Lower  trace: 
Pxr>  8 mTorr;  P u,  120  mTorr.  The  solid  line  corres- 
|)onds  to  the  recorded  profile  and  large  dots  represent 
the  calculated  profiles  corresponding  to  the  best  fit. 

The  three  dotted  lines  represent  the  contributions  from 
the  three  terms  of  Eq.  (3):  Curve  1 corresponds  to  the 
narrow  resonance  (first  term),  curve  2 corresponds  to 
the  Gaussian  background  (second  term)  arising  from 
Kr*-Kr  collisions,  and  curve  1 corresponds  to  the  con- 
tribution of  VCC  (third  term)  arising  from  Kr*-He 
and  Kr*-Ar  collisions. 


values  of  the  two  parameters  y and  r0  were  varied 
so  as  to  minimize  the  rms  error  of  the  fit  of 
Eq.  (3)  to  the  data.  Typical  examples  for  He  and 
Ar  are  shown  on  Fig.  2.  For  the  entire  pressure 
range  studied  the  data  could  be  fit  very  well  ( i.e . 
with  a normalized  rms  error  -0.05)  by  Eq.  (3). 

The  values  of  r = ro'  obtained  for  the  various 
perturber  pressures  are  plotted  on  Fig.  3.  There 
is  an  increase  of  r with  pressure  similar  to  the 
one  observed  for  Kr*-Kr  collisions8;  accordingly, 
we  conclude  that  r is  determined  by  the  transit 
times  of  atoms  in  the  beams.  The  findings  are 
consistent  with  the  premise  that  heavier  perturbers 
are  more  effective  in  increasing  the  transit  time 
than  lighter  ones.  The  two  curves  converge  to  the 
same  transit  time  at  low  pressure  (r0  = 4 ps)  in 
accordance  with  the  value  measured  in  pure  Kr 
at  8 mTorr;  for  the  highest  pressures  they  tend 
roughly  to  the  linear  variations  given  by  the 
classical  law  of  diffusion.8 

The  FWHM  of  the  VCC  contributions  (third  term 
in  Eq.  3)  as  a function  of  perturber  pressures  is 
shown  on  Fig.  4.  As  may  be  seen,  both  widths^ 
approach  the  thermal-equilibrium  value  af/vT 
= 510  MHz  as  the  perturber  pressure  is  increased. 
The  heavier  Ar  atoms  produce  thermalization  at 
a lower  pressure  than  the  lighter  He  ones  (see 
also  Fig.  1)  consistent  with  the  hard- sphere  col- 
lision model.  A value  for  these  widths  obtained 
at  very  low  pressure  can  be  deduced  from  Eq.  (3) 
by  taking  the  n = l term  in  the  sum.  These  values 
are  indicated  in  Fig.  4 by  short-dashed  lines. 

We  have  tried  various  values  of  a around  the 
values  mentioned  above,  o =0.93  for  Kr*-He  and 
a =0.6  for  Kr*-Ar.  If  a is  chosen  too  small  the 
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FIG.  4.  Variations  of  the  width  of  the  VCC  contri- 
butions [third  term  in  Eq.  (3)  of  text)  vs  perturber 
pressure  (dots  correspond  to  He  and  crosses  to  Ar 
perturbers).  The  horizontal  dashed  line  corresponds 
to  the  thermal-equilibrium  limit.  The  short  horizontal 
dashed  lines  correspond  to  the  width  of  the  n = 1 contri- 
bution in  Eq.  (3). 


resulting  am  is  too  large  to  fit  the  data  (for  example, 
a =0.85  for  He  is  a lower  limit);  on  the  other  hand 
if  a is  chosen  too  large,  the  data  could  be  fit  only 
with  unacceptably  small  values  for  2r  (for  He,  a 
= 0.98  is  an  upper  limit).  We  conclude  that 

a =0.93  ±0.03  for  Kr*-He, 

a =0.6  ±0.2  for  Kr*-Ar, 

where  the  error  limits  are  obtained  by  those  values 
of  a leading  to  a normalized  rms  error  of  0.2  in 
the  fitting  procedure. 

The  values  of  r2  that  we  have  chosen  have  been 
calculated  using  the  ground- state  kinetic- theory 
radii  tabulated  by  Allen.14  Values  of  r2  determined 
in  the  same  manner  but  by  choosing  for  Kr*  meta- 
stables a radius  corresponding  to  the  excited 
electron17  have  also  been  tried.  Results  are  less 
consistent;  first,  the  transit  time  takes  on  a 


nearly  constant  value  (=3  ps)  over  the  range  of 
pressure  and,  second,  the  width  of  the  narrow 
resonance  is  too  small.  Additional  experiments 
with  various  perturber-active-atom  mass  ratios 
may  help  to  clarify  the  manner  in  which  the  effective 
excited-state  radius  should  be  chosen. 

We  have  tried  also  to  fit  the  profiles  by  leaving 
free  the  three  parameters  y,  r0,  and  r3.  For  He 
the  variations  of  r2  with  pressure  showed,  as 
expected,  a linear  increase  at  low  pressure,  but  a 
decrease  at  high  pressure  (p>200  mTorr).  This 
deviation  from  the  linear  increase  was  attributed 
to  the  fitting  procedure  used  since  at  high  pressures 
VCC  contribute  well  into  the  wings  of  the  thermal 
distribution  and  it  is  difficult  to  distinguish  between 
Kr*-Kr  and  Kr*-He  contributions.  The  low- 
pressure  regime  gave  r2  =0.95x  10*p  sec"1  con- 
firming our  choice  for  r2  determined  from  kinetic 
theory  and  Allen’s  values. 

IV.  CONCLUSION 

We  have  arrived  at  a consistent  picture  of  VCC 
in  the  metastable  4p55s[  | )2  state  of  krypton  caused 
by  argon  and  helium  perturbers,  based  on  a hard- 
sphere  collision  model.  Future  experiments  of 
this  type  in  two-  or  three- level  systems  may  be 
expected  to  provide  a more  sensitive  probe  of 
the  short-range  collisional  interaction,  while 
providing  a new  method  for  studying  excited-state 
scattering. 
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Resume.  — Une  experience  d’absorption  saturee  a deux  lasers  conduite  sur  la  raiedeKrla 
A = 557  ran  a permis  I'etude  des  collisions  elastiques  Kr*-He.  Pour  une  classe  donnee  d’atomes 
pompes  hors  resonance  (n,  # 0)  on  peut  mettre  en  evidence  la  persistance  de  la  vitesse  initiate  apres 
n collisions  elastiques. 

Abstract.  — Elastic  Kr*-He  collisions  have  been  studied  using  a two  laser  saturated-absorption 
experiment  on  the  A = 557  nm  line  of  Krl.  With  the  pump  laser  detuned  from  the  atomic  resonance, 
one  can  follow  the  collisional  relaxation  of  Kr*  atoms  from  a non-thermal  distribution  centred  at 
vr  # 0 back  towards  a thermal  distribution. 


! 
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In  saturated-absorption  experiments,  a probe  beam 
monitors  the  longitudinal  velocity  distribution  of 
active  atoms  that  are  selectively  excited  by  a pump 
beam.  Velocity-changing  collisions  (V.C.C.)  alter  this 
longitudinal  velocity  distribution  and  produce  a 
corresponding  modification  of  saturated -absorption 
profiles.  Strong  or  thermal izing  collisions  give  rise  to  a 
broad  Gaussian  background  superimposed  on  the 
narrow  resonance,  i.e.,  following  a collision,  on  ave- 
rage. the  active  atom's  velocity  is  independent  of  its 
initial  velocity,  v2.  and  is  determined  by  the  thermal 
distribution  [1-4].  Weak  collisions  give  rise  to  a dis- 
tortion of  the  shape  of  the  narrow  resonance ; they  are 
collisions  with  large  impact  parameters  producing 
velocity  changes  smaller  than  the  width  of  the  hole 
burned  in  the  velocity  distribution  [5-8].  Elastic  hard 
sphere  collisions  between  light  perturbers  and  heavy 
active  atoms  lead  to  a non-Gaussian  background  signal 
owing  to  the  velocity  redistribution  that  they  induce 
[9,  10].  The  degree  of  thermalization  following  an 
average  collision  is  sometimes  referred  to  as  the 
persistence  of  velocity. 

In  a previous  publication  [9],  we  studied  the 
effect  of  V.C.C.  on  saturated-absorption  profiles 
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(A  = 557  nm  of  K rl)  as  a function  of  the  active  atom 
to  perturber  mass  ratio.  In  that  work  the  same  laser 
was  used  to  saturate  and  to  probe  the  transition. 
During  the  scanning  of  the  laser  frequency,  the  satura- 
tor selects  successive  classes  of  velocity  vz ; V.C.C. 
drift  and  broaden  the  corresponding  distributions  of 
active  atoms  but  their  effect  is  necessarily  detected  by 
the  probe  at  — v2  only,  hence  symmetrical  profiles 
centred  around  v2  = 0.  To  directly  view  the  drift  of 
active  atoms,  it  is  necessary  to  select  than  at  a given 
v.  # 0 and  to  tune  the  probe  over  the  velocity  dis- 
tribution as  modified  by  V.C.C. 

To  this  end.  we  have  performed  a two  laser  satu- 
rated-absorption experiment  on  the  same  transition  of 
Krf4p5  5s[3/2]2-4p5  5p'[l/2],)  with  the  V.C.C.  pro- 
vided by  the  perturbers.  The  first  dye  laser,  at  fre- 
quency v, , is  used  to  select  by  saturation  a given  class  of 
longitudinal  velocities  far  from  zero;  to  accomplish 
this,  it  is  detuned  from  the  atomic  resonance  at  v0 
by  an  amount  comparable  to  the  Doppler  width 
(Avd  = 720  MHz).  The  second-dye  laser  is  used  to 
probe  both  saturated -absorption  (counter  propagat- 
ing waves)  and  linear  absorption.  Fabry-Perot  fringes 
are  used  to  calibrate  the  frequency  scale  of  the  probe. 
A typical  recording  is  shown  on  figure  1.  Trace  I is 
the  linear  absorption  profile,  centred  at  v0.  Trace  II 
is  a saturated-absorption  profile  obtained  in  pure 
krypton,  for  a detuning  A = v,  - v0  = 600  MHz 
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satisfies  the  equation  : 
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with 


G(p,  -*  0)  = 6(v,  - 


vely.  The  quantity  T is  the  total  rate  of  population 
density  loss,  T = T0  + T,  + r2  where  T0  is  the 
inverse  lifetime  of  the  state  population  under  consi- 
deration, which  is  here  the  time  during  which  the 
atoms  stay  in  the  fields  [9].  fVt(v2  -*  v')  is  the  pro- 
bability density  for  a thermalizing  event  to  occur  per 
unit  time; 
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W'lU’z  - V')  = tv,(v;)  = — - =r-  exp  - (t>;V) 
sJGu 


Fig.  1.  — Trace  I : Linear-absorption  profile,  centred  at  v0. 
Trace  II  : Saturated-absorption  profile,  obtained  in  pure  krypton 
for  a detuning  A = 600  MHz  of  the  saturator.  Trace  III  : Saturated- 
absorption  profile,  obtained  in  a mixture  Kr-He  (P(Kr)  ^ 10  tntorr. 
/’(He)  - 450  mtorr)  for  a defining  A = 600  MHz.  The  squares 
represent  a profile  calculated  from  expression  (2)  with  the  values  of 
the  parameters  indicated  in  the  text.  Trace  (a)  is  the  Lorentzian 
centred  at  - A : Trace  (b)  is  the  contribution  due  to  thermalizing 
Kr’-Kr  collisions,  centred  at  A ' = 0;  Trace  (c)  is  the  contribution 
due  to  Kr*-He  V.C.C 


where  u is  the  most  probable  speed  of  the  thermal 
distribution.  fV2(v'i  -*  K)  is  the  probability  density 
for  non-thermalizing  events ; we  choose  the  model  of 
Keilson  and  Storer  [12]  to  express  : 


W2(v’;  - v’2)  = 


1 


Jn\u 


exp  - [(v'  - otv?)2/(Au)2] 


of  the  saturator ; it  shows  a narrow  resonance  centred 
at  A ' = — A and  a Gaussian  background  centred  at 
A ’ = 0.  and  attributed  to  Kr*-Kr  collisions.  Trace  III 
shows  a saturated-absorption  profile  obtained  in  a 
mixture  Kr-He,  for  the  same  detuning.  The  difference 
in  shape  between  trace  II  and  III  is  attributed  to 
Kr*-He  V.C.C.  whose  main  contribution  lies  between 
— A and  0. 


where  a is  the  mean  strength  of  V.C.C.  related  to  the 
r.m.s.  change  of  velocity  A u by  (Au)2  = (1  - a2)  u 2. 
In  this  model  r2  can  be  fairly-well  approximated  by 
the  rate  of  collisions  in  a hard  sphere  model  and  for  a 
small  ratio  of  perturber  to  active  atom  mass  m/M  « 1 . 
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Equation  (1)  can  be  solved  by  iteration  and  one 
obtains  [10]  : 


G(vz  -» t-;)  = 


S(vz  - >0  rt  wtuo 


r0(r0  + r.) 


Interpretation  of  profiles.  — To  calculate  the 
resulting  profile  in  the  presence  of  elastic  collisions, 
one  considers  their  effect  on  the  phase  of  active  dipoles 
(introduced  via  classical  broadening  constants)  and 
their  effect  on  the  velocity  associated  with  the  atomic 
state  populations  (introduced  via  a collision  propa- 
gator [11]).  The  propagator  G which  propagates  a 
population  density  n(vt,  0)  (i\  = longitudinal  velocity) 
to  the  value  n(i>',  t)  via 


~ x cxp  - \- — 

r .Ir/v/nAu.  L Au.  J 


where 


(Au.)2  = (I  — a2")  u 2 


and 
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Saturated-absorption  profiles  are  given,  in  a two 
laser  experiment  (counter  propagating  beams)  by  the 
following  expression  (at  very  low  saturation)  [10] 


» j | f 


+ (A  - kv,) 
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f + (A ' + kv't) 
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(I) 


Both  thermalizing  and  non  thermalizing  collisions 
are  assumed  to  occur  with  rates  r,  and  T 2 respecti- 


where  ? is  the  homogeneous  width  (H.W.H.M.),  A the 
detuning  of  the  saturator  and  A ' the  detuning  of  the 
probe. 

G(v,  -»  r')  is  the  sum  of  three  terms,  hence  the  total 
profile  is  also  a sum  of  three  terms ; it  can  be  expressed 
in  the  Doppler  limit  (f  < ku)  as  : 
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where  2 f.  = 7(1  + a"). 

The  first  term  in  (2)  is  a Loren tzian  corresponding  to 
atoms  which  have  not  undergone  velocity  changes ; it 
is  centred  at  - A with  total  width  (F.W.H.M.)  of  4 7. 
The  second  term  is  a Gaussian  corresponding  to 
atoms  which  have  undergone  at  least  one  thermalizing 
collision  ; it  is  centred  at  A ' = 0 and  has  the  Doppler 
width.  The  third  term  is  a sum  of  Voigt  integrals  and 
corresponds  to  atoms  which  have  undergone  velocity 
changes  but  no  thermalizing  event ; it  is  asymmetrical 
and  expresses  the  progressive  loss  of  memory  of  the 
initial  velocity  following  n collisions. 

The  profile  shown  on  figure  1 has  been  compared  to 
a theoretical  one  (expression  (2))  by  use  of  a least 
square  fit.  a has  been  chosen  equal  to  0.94  correspond- 
ing to  the  predicted  value  {1 1] ; 7 has  been  evaluated  at 
18  MHz  according  to  the  experimental  situation. 
Thus,  profiles  depend  on  two  parameters  only.  T 2/r 
and  T,  r/r0(r0  + r,)  which  determine  the  inten- 
sities of  the  two  components  due  to  V.C.C.  compared 
to  the  intensity  of  the  Lorentzian.  Here  we  have 
chosen  to  fix  f2/r  and  we  have  looked  for  a value  of 
T,  r/r0(r0  + r,)  which  gives  account  of  the  shape 
of  the  profile ; the  same  procedure  has  been  used  for 
several  different  values  of  T 2/r.  n being  limited  to 
sixty.  The  best  fit.  determined  by  the  minimum  value 
of  the  mean  square  deviation,  has  been  obtained  for  a 
well-defined  couple  of  parameters : in  particular,  the 
shape  of  the  profile  is  very  sensitive  to  the  value  of 

r2/r. 

An  example  of  the  fit  is  shown  on  figure  I ; on 
trace  HI.  dots  represent  the  calculated  profile,  sum  of 
the  three  terms  of  equation  (2),  each  of  them  being 
represented  by  a thin  curve  (trace  a,  b,  c).  The  asym- 
metry of  the  component  due  to  Kr*-He  V.C.C.  shows 
the  non-thermalizing  nature  of  these  collisions;  this 
is  better  shown  on  figure  2,  where  different  Voigt 
profiles  related  to  different  values  otn  are  represented. 


A profile  associated  to  n represents  the  contribution  of 
atoms  which  have  undergone  a mean  velocity  change 
A u„  after  n collisions. 


Fig.  2.  — Contribution  of  different  Voigt  profiles  with  different  n 
values  to  the  component  associated  to  Kr*-He  V.C.C.  (Trace  (c)  of 
figure  I).  Positions  of  these  different  components  are  marked  by 
dashes.  Note  that  the  profile  for  n = I is  already  shifted  of  an 
amount  k \u 


The  fit  yields  the  following  values  for  the  para- 
meters : r i/r2  = 0.04  and  r jr2  =1.6  in  close 
agreement  with  those  obtained  in  [9]  on  Kr  saturated- 
absorption  with  one  laser.  However  we  can  notice 
that  the  fit  of  the  asymmetrical  profile  is  much  more 
sensitive  to  certain  parameters;  in  particular,  the 
spread  of  a values  around  0.94  indicates  that,  actually. 
a is  defined  in  a narrow  range  : a = 0.94  ± 0.01 . 
Thus,  our  results  indicate  that  a simple  model  of  hard 
sphere  collisions  can  consistently  explain  the  effects  of 
V.C.C.  on  saturated-absorption  profiles. 

D.  Merle  is  warmly  acknowledged  for  her  contribu- 
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Abstract : 

A fundamental  physical  problem  is  the  determination  of  atom  atom,  atom  molecule  and  molecule  molecule  differential  and 
total  scattering  cross  sections.  In  this  work,  a technique  for  studying  atomic  and  molecular  collisions  using  spectroscopic  line  shape 
analysis  is  discussed.  Collisions  occurring  within  an  atomic  or  molecular  sample  influence  the  sample's  absorptive  or  emissive  properties. 
Consequently  the  line  shapes  associated  with  the  linear  or  non-linear  absorption  of  external  fields  by  an  atomic  system  reflect  the 
collisional  processes  occurring  in  the  gas.  Explicit  line  shape  expressions  are  derived  characterizing  linear  or  saturated  absorption 
by  two-  or  three-level  "active"  atoms  which  are  undergoing  collisions  with  perturber  atoms.  The  line  shapes  may  be  broadened,  shifted, 
narrowed,  or  distorted  as  a result  of  collisions  which  may  be  "phase-interrupting"  or  “velocity-changing"  in  nature.  Systematic  line 
shape  studies  can  he  used  to  obtain  information  on  both  the  differential  and  total  active  atom-perturher  scattering  cross  sections. 
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1.  Introduction 


I had  the  opportunity  of  working  with  Professor  Willis  Lamb  for  a period  of  five  years,  1966- 
1971.  For  three  of  those  years,  I was  a graduate  student  and  for  the  remaining  two  years  an  Instruc- 
tor at  Yale  University.  Professor  Lamb  not  only  guided  me  in  research,  but  also  left  me  with  an 
imprint  of  his  philosophy  of  physics.  It  is  with  sincere  gratitude  that  1 dedicate  this  article  to  him. 

There  is  a certain  mythology  that  has  arisen  concerning  Professor  Lamb.  It  has  been  rumored 
that  one  should  never  utter  the  word  "photon"  or  words  "Green’s  function"  in  his  presence.  In  my 
relationship  with  Professor  Lamb,  I never  found  this  to  be  the  case.  As  long  as  the  terminology 
clearly  described  the  physical  processes  under  discussion,  there  was  never  an  objection  to  one 
word  or  another.  However,  1 think  it  is  fair  to  say  that  Professor  Lamb  is  somewhat  impatient 
with  those  who  use  terminology  or  formalism  at  the  expense  of  physics. 

Stressing  the  physical  aspects  of  a given  problem  was  always  foremost  in  Professor  Lamb’s 
mind.  He  would  often  use  classical  or  well  known  quantum-mechanical  analogs  to  try  to  explain 
a new  result,  rather  than  simply  accept  the  fact  that  "it  follows  from  the  mathematics".  I should 
say  that  Professor  Lamb’s  greatest  influence  on  me  was  to  convey  this  need  to  get  “inside"  a problem 
to  understand  the  underlying  physical  principles. 

If  there  is  one  interaction  I shall  never  forget,  it  is  the  eight  or  nine  month  period  that  we 
spent  going  over  my  thesis.  Each  sentence  was  subjected  to  careful  scrutiny  and  it  was  not  uncom- 
mon to  change  a phrase  two  or  three  times  before  Professor  Lamb  found  it  acceptable.  The  general 
theme  of  his  critical  reading  of  the  manuscript  was  “If  it’s  worth  putting  in,  it’s  worth  explaining 
it  clearly”.  While  I admit  that  I have  been  guilty  of  not  following  this  maxim  from  time  to  time, 
it  is  the  type  of  philosophy  I try  to  pass  along  to  my  students. 

I have  tried  to  incorporate  Professor  Lamb's  philosophy  (or,  more  precisely,  my  perception  of 
Professor  Lamb’s  philosophy)  in  writing  this  article.  The  article  contains  a summary  of  the  research 
I have  carried  out  over  the  past  several  years  related  to  spectroscopic  studies  of  collision  effects 
in  atomic  and  molecular  systems.  It  is  not  a review  article,  but  rather  a personal  view  of  the  subject 
matter  in  which  the  physical  aspects  of  the  various  problems  are  stressed.  Extensive  bibliographies 
may  be  found  in  two  recent  reviews  of  this  subject  area  [1,  2]. 

A fundamental  problem  in  atomic  and  molecular  physics  is  the  determination  of  differential  and 
total  cross  sections  for  atom-atom,  atom-molecule  or  molecule-molecule  scattering.  The  tradi- 
tional method  for  studying  scattering  phenomena  involves  the  use  of  crossed  atomic  or  molecular 
beams.  Crossed  beam  experiments  are  usually  restricted  to  ground  or  metastable  levels;  excited 
state  densities  are  too  small  to  give  rise  to  detectable  scattering  signals.  Recently,  however,  signifi- 
cant excited  state  population  has  been  achieved  by  using  a laser  to  excite  atoms  in  the  scattering 
region. 

In  this  paper,  I discuss  an  alternative  method  for  studying  atomic  and  molecular  collisions. 
Consider  a bulb  containing  a gas  of  classical  oscillators  (“active”  atoms)  and  perturbers.  The 
oscillators  absorb  or  emit  radiation.  There  is  a characteristic  profile  associated  with  these  processes 
having  a width  that  reflects  both  the  Doppler  broadening  and  natural  damping  of  the  oscillators. 
When  an  oscillator  collides  with  a perturber,  its  phase  and  velocity  are  altered.  These  effects  give 
rise  to  a modification  of  the  absorption  or  emission  profile.  Consequently,  line  shape  studies  can 
be  used  to  implicitly  monitor  collision  processes  in  the  gas.  One  can  determine  the  active  atom- 
perturber  collision  rate  and  total  cross  section  by  this  method. 

Actually  the  situation  is  a bit  different  if  one  uses  two-level  atoms  rather  than  classical  atoms 
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to  represent  the  active  atoms.  The  linear  superposition  of  the  two  states  induced  by  an  external 
field  or  the  vacuum  field  produces  a dipole  moment  in  these  atoms,  in  analogy  with  the  dipole 
moment  of  the  classical  oscillators.  However,  the  active  atom-perturber  scattering  interaction 
may  be  state  dependent  for  two-level  atoms  (i.e.  dependent  on  which  state  the  two-level  atom  is  in) 
whereas  it  was  state-independent  for  classical  oscillators.  When  a two-level  atomic  dipole  enters 
a collision,  which  trajectory  should  it  follow?  There  is  no  classical  answer  to  this  question;  the 
quantum-mechanical  answer  is  discussed  in  section  2.  The  study  of  the  linear  absorption  or  emis- 
sion profiles  associated  with  two-level  atoms  can  provide  values  for  the  state-dependent  active 
atom-perturber  collision  rates  and  total  cross  sections. 

Linear  spectroscopy  can  provide  information  on  total  cross  sections  only.  To  learn  something 
about  differential  scattering  cross  sections,  one  may  use  saturation  spectroscopy.  A monochromatic 
pump  field  excites  only  those  atoms  with  a specified  velocity  component  in  the  direction  of  the 
field's  propagation  vector.  The  excited  atoms  undergo  collisions  in  which  their  velocities  are 
altered,  and  these  velocity  changes  can  be  observed  by  subjecting  the  atoms  to  a probe  field  acting 
on  the  same  or  a coupled  transition.  One  does  not  obtain  a true  differential  scattering  cross  section 
in  such  experiments,  but  rather  a differential  scattering  cross  section  averaged  over  perturber 
velocities  and  those  active  atom  velocities  not  selected  by  the  pump  field.  Still,  some  semi-quanti- 
tative data  on  the  excited  state  differential  scattering  cross  section  may  be  inferred  from  the  line 
shape. 

In  light  of  the  comments  given  above,  it  is  reasonable  to  believe  that  linear  and  saturation 
spectroscopy  should  be  effective  probes  of  atomic  and  molecular  collisions.  The  remainder  of  this 
paper  outlines  a systematic  method  that  permits  one  to  extract  appropriate  collision  rates  and 
cross  sections  from  the  line  shapes  associated  with  atomic  and  molecular  systems.  In  section  2, 
an  equation  for  the  collisional  time  rate  of  change  of  atomic  density  matrix  elements  is  derived. 
It  is  shown  that  for  state  dependent  collision  interactions,  it  is  not  possible  to  use  a classical  velocity 
to  describe  the  atomic  center-of-mass  motion.  In  section  3 the  linear  absorption  line  shape  for  a 
two-level  system  is  calculated  and  in  section  4,  the  saturation  spectroscopy  profile  for  a three- 
level  system  is  derived.  Phenomena  such  as  collision  broadening,  shifting,  and  narrowing  are 
discussed  in  these  sections.  Section  4 also  illustrates  the  way  in  which  saturation  spectroscopy  can 
be  used  to  provide  information  on  excited  state  differential  scattering  cross  sections.  In  section  5 
additional  problems  relating  to  inelastic  collisions,  level  degeneracy  and  resonant  exchange  colli- 
sions are  discussed  briefly. 


2.  Collisions 

To  study  collisions  between  “A”  and  “P”  atoms  using  spectroscopy  as  a collision  probe,  one 
starts  with  a bulb  containing  a mixture  of  A and  P atoms.  The  A (“active”)  atoms  are  allowed 
to  interact  with  some  external  fields  while  undergoing  collisions  with  the  P (perturber)  atoms. 
One  monitors  the  active  atom-field  interaction  (e.g.  absorption  of  an  applied  field)  as  a function 
of  perturber  pressure  to  obtain  information  on  A-P  collisions.  It  is  assumed  that  the  perturbers 
interact  negligibly  with  the  external  fields. 

The  line  shapes  associated  with  the  active  atom-field  interactions  reflect  all  collisional  processes 
occurring  in  the  gas.  It  is  possible,  however,  to  choose  the  experimental  parameters  such  that 
binary  A-P  collisions  are  the  only  type  of  collisions  that  contribute  significantly  in  line  shape 
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formation.  First,  a dilute  (<  500  Torr)  gas  must  be  used  to  insure  that  three-body  collisions  are 
not  important.  In  this  binary  collision  pressure  region,  the  time  between  collisions  is  much  greater 
than  the  duration  of  a collision.  Second,  a large  ratio  of  perturbers  to  active  atoms  must  be  main- 
tained to  guarantee  that  it  is  A-P  rather  than  A-A  collisions  which  are  dominant. 

When  the  ratio  of  perturber  to  active  atoms  is  large,  there  is  an  additional  simplification;  the 
perturbers'  velocity  distribution  is  essentially  unaffected  by  the  relatively  few  collisions  they  under- 
go with  active  atoms.  The  perturbers  may  be  thought  to  constitute  a thermal  bath  that  is  unchanged 
by  collisions.  In  the  presence  of  this  thermal  bath  the  evolution  of  the  active  atoms’  velocity 
distribution  may  be  described  by  a linear  transport  equation  rather  than  the  more  general  Boltz- 
mann equation.  The  appropriate  transport  equation  is  derived  below. 

The  situation  can  be  more  complicated  if  one  wishes  to  study  A-A  collisions.  In  that  case,  it  is 
possible  for  atoms  to  resonantly  exchange  excitation  during  a collision.  Resonant  collisions  are 
discussed  briefly  in  section  5.  In  sections  2-4,  the  discussion  is  limited  to  binary,  elastic,  foreign 
gas  (non-resonant)  collisions. 

2.1.  One-level  system 

The  simplest  system  one  can  envision  involves  one-level  active  atoms.  Physically,  these  one-level 
active  atoms  could  correspond  to  ground  state  atoms  which  are  injected  with  some  non-equilibrium 
velocity  distribution  into  the  perturber  bath.  For  these  one-level  active  atoms,  collisions  result 
solely  in  a change  of  velocity  (fig.  1).  In  fact,  the  system  may  be  though  of  classically  as  a thermal 
classical  gas  of  P atoms  into  which  few  non-thermal  A atoms  have  been  introduced.  The  velocity 
of  A atoms  is  then  governed  by  the  linear  transport  equation 

= - rA(r)/A(r,  t)  + J d v'W^v'  - v)/A(»>',  t),  (21) 

where /A(r,  f)  is  the  distribution  function  of  the  A atoms,  ITa(i>'  -♦  v)  is  the  collision  kernel  giving 
the  probability  density  per  unit  time  that  an  A atom  undergoes  a velocity  change  from  v'  to  v and 

rA(r)  = Jdr'W>-r)  (2.2) 

is  the  rate  of  collisions.  Eq.  (2.1)  multiplied  by  dv  describes  the  flow  of  A atoms  out  of  and  into  the 
velocity  range  between  v and  v + dv. 

It  is  not  too  difficult  to  arrive  at  a quantum-mechanical  transport  equation  analogous  to  eq.  (2.1). 
Provided  that  the  range  of  the  collision  interaction  is  much  less  than  the  collisional  mean  free 
path  (binary  collision  approximation),  collisions  result  solely  in  a change  in  the  average  velocity  v 
of  a wave  packet  associated  with  an  active  atom.  The  quantum  mechanical  collision  kernel  is 


Fig  I Classical  picture  of  a collision  between  a “one-level"  active  atom  A and  a perturber  P.  The  atom's  velocity  is  changed  from 
r to  r as  a result  of  the  collision. 
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O l’r|  J\(v,  Ol2)  where  is  the  perturber  density,  vr  the  active  atom-perturber  relative 
velocity,  /A( iv  -*  fr)  the  scattering  amplitude  for  active  atom-perturber  collisions,  and  the  average 
is  over  all  possible  perturber  velocities  consistent  with  conservation  of  momentum  and  energy. 
Explicitly 


MA<r'  -»*’)=•  Jdlr  d,’rM/p<*’P^r  ’<*  v<  + v' 


V.  = V 


V,  = V - P_, 


m 

— v + 


- v,)fK(v, 


vr)|2, 

(2.3) 


where  r (or  r')  is  the  active  atom  velocity,  >’p  (or  rp)  the  perturber  velocity,  H^Pp)  the  perturber 
velocity  distribution,  m the  active  atom  mass,  mp  the  perturber  mass  and  p the  reduced  mass. 

Similarly,  the  collision  rate  is  given  by  <./F i>r<xA(i»r)>  where  a is  the  total  elastic  cross  section 
for  active  atom-perturber  scattering.  As  in  eq.  (2.2)  one  has 


r» 


= Jdr'W^r  - r’)  = .4'  jdvpWp(vp)vr<r, 


K)- 


(2.4) 


The  quantum  mechanical  transport  equation  has  the  same  form  as  eq.  (2.1),  with  the  kernel  and 
rate  now  given  by  eqs.  (2.3)  and  (2.4),  respectively,  and  the  distribution  function  now  describing 
a wave  packet  with  average  velocity 

Since  it  will  aid  in  one’s  understanding  of  collisions  in  two-level  systems,  I should  like  to  outline 
a simple,  non-rigorous  quantum-mechanical  derivation  of  eq.  (2.1).  Imagine  an  active  atom 
wave  packet  of  cross-sectional  area  <j  and  average  momentum  p'  approaching  a stationary  pertur- 
ber (<r  square  of  the  collision  interaction  range).  This  packet  could  be  described  by  a wave 
function 


ip(R.t-)  = e<i/*,p  *.4(ff.r). 


(2.5) 


where  A(R.  t ) properly  localizes  the  packet.  Following  the  collision,  the  packet  consists  of  an 
unscattered  component  ,4(/?.  t* ) plus  a spherically  scattered  wave  S[R.  r + ) (fig.  2) 


>J/(R.  f + ) = e'1’"'’  *A(R,r+)  + ^ V)e(i/*>p  *S(/?,  r + ), 


(2.6) 


where 


MR.  t + ) = exp[-(i«ir'2/2)i)(t+  - r")]A[/?  - v'{t*  - t'\l  ]. 


(2.7a) 


2 Quantum -mechanical  picture  of  a collision  between  an  active  atom  A and  a perturber  P The  active  atom  is  represented  by 
a localized  wave  packet  Following  a collision,  the  active  atom  consists  of  an  unscattered  component  plus  a spherical  shell  representing 
the  scattered  wave. 
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S(R.t  + ) = exp  [— (imr'2/2ft)(f+  - t~)]A[Rv  - y'(f  + - r),f]. 


(2.7b) 


again  give  the  proper  localization  of  the  packet.  The  change  in  probability  density  resulting  from 
this  collision  is 

Sp(R,t\r)  = \MR,t+)\2  - \<HR,r)\2 

= [|/4(/?,t  + )|2  - | /!(/?,  r~)|2]  + 14-9.  *-p*'A(R,  t+)S{R,t  + )* 

R 


/A(y  -»  y)  iilt|lrll-|, 


p R)A(R.t  + )*SiR,t  + ) + 


+x  t \Ja(p'^v)\2 


|S(IU+)|2.  (2.8) 


In  the  limit  that 
St  — t+  — r 

goes  to  zero,  the  lead  term  in  (2.8)  will  vanish,  reflecting  the  fact  that  the  atomic  velocity  is  unchang- 
ed if  no  collision  occurs. 

The  remaining  terms  specify  the  change  in  velocity  from  an  initial  value  v'  = p'/m  to  a final 
value  y = v'R  For  a perturbcr  density  Jf.  the  density  of  particles  scattered  per  unit  time  into  a 
solid  angle  d£2,  is 

[<(p(y,  y',  f,  dfl/'dr]  dQ,  = (.  Vv'aSx){R2  dfi,./<r)<5p(/?,  f + , t~)/Sx 
= jre ^ e<i/»><y  * - P- *)A{ R , * )S(R  t * ^ 

+ — V*  V)  ei!h)^R-pmA(R.  t+)*S(R,t+)  + ~Vd7  |S(tf.f  + )|2l  R2dn...  (2.9) 
R R j 

It  is  necessary  to  integrate  eq.  (2.9)  over  all  possible  incident  velocity  directions  dQ,...  In  performing 
this  integration,  one  uses  the  fact  that  R > h/p'  (i.e.  in  the  binary  collision  pressure  region,  the 
(collisional  mean  free  path)  > = h/p'  - it  follows  that,  for  all  but  a negligible  range  of  R.  R > /B) 

to  set 

JdD,  Z(y'  - y) emur  *-pK)  * ~ /(y  - y).  (2.10) 

The  rapid  oscillation  of  the  phase  factor  leads  to  destructive  interference  in  all  but  the  forward 
direction.  Dropping  the  dCl,.  and  integrating  (2.9)  over  dQ,.  making  use  of  (2.10).  one  obtains 


n [,/A(y  - v)M(*,  f + )S(l?,  f + )*  — /A(y  - yM(tf.  r + )*S(tf,  f + )] 

ox  \ ip  ) 

+ JT\ r Jd!!,.  | /A(y'  - y)|2  |S(lf,  f + )|2.  (2.11) 

Using  (2.7)  and  taking  the  limit  St  -*  0,  Sx  -*  0,  one  arrives  at  the  quantum  mechanical  transport 


= .Ur  I- 
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equation 

dp(v,  t)/dt  | co„  = - TA(v)p(i\  t)  + ('dr' v'  -*•  v)p(v',  t), 

j 

where 

WA(v'  -+  v)  = Jfv  | /A(v'  -*  v)|2iT2<5(i’  - v'), 

rA(»-)  = jr [/A(r  -*■ v)  -fA(v  - *)*]• 

The  optical  theorem  can  be  used  to  transform  (2.14)  into 

rA(^)  = eV  t><xA(v). 


(2.12) 

(2.13) 

(2.14) 

(2.15) 


Eq.  (2.12)  with  eqs.  (2.13)  and  (2.15)  generalized  to  allow  for  moving  perturbers  is  seen  to  be 
equivalent  to  eq.  (2. 1 ) with  the  kernel  and  rate  given  by  eqs.  (2.3)  and  (2.4),  respectively.  In  a classical 
or  quantum-mechanical  picture,  collisions  for  one  level  systems  can  be  described  by  a transport 
equation.  Naturally,  some  measurement  technique  would  be  needed  to  monitor  the  active  atom 
velocity  distribution.  It  is  implicitly  assumed  above  that  the  time  scale  of  the  measurement  is 
much  greater  than  the  collision  time. 


2.2.  T wo-level  system 

In  contrast  to  beam  experiments  where  scattering  from  one-level  atoms  is  typical  (i.e..  ground 
state  ground  state  or  ground  state-metastable  state  collisions),  bulk  experiments  generally  involve 
scattering  of  active  atoms  which  are  in  a linear  superposition  of  internal  states.  In  the  beam 
experiments  one  directly  measures  the  scattered  particles;  however  in  laser  spectroscopic  col- 
lisions studies,  one  uses  absorption  or  emission  profiles  to  indirectly  monitor  collisions  in  the 
gas.  The  radiation  fields  that  are  either  absorbed  or  emitted  by  the  gas  produce  a linear  super- 
position state  in  the  active  atoms  (i.e.,  an  oscillating  atomic  dipole).  Consequently,  collisions  in- 
volving two-level  atoms  must  be  analyzed  to  properly  interpret  spectroscopic  line  shapes. 

Consider  the  two  level  atom  of  fig.  3 in  which  the  frequency  separation  of  levels  1 and  2 is  to. 
The  atom  is  assumed  to  be  in  a coherent  superposition  state  produced  by  a radiation  field  and 
undergoes  collisions  with  ground  state  perturber  atoms.  For  the  present,  I assume  that  collisions 

• re  ,l,isth  (inelastic  collisions  are  discussed  in  section  5),  not  possessing  sufficient  frequency 
•mponcnls  or  energy  to  induce  transitions  between  states  1 and  2.  Since  the  collision  duration  rc 

* p>.  ,IK  * 10  1 •’  m.*c.  collisions  contain  frequency  components  up  to  toQ  = (rc)~ 1 % 1012  sec- 

t N»  ’ no  ■ corresponding  to  thermal  energies  is  % 1013  sec" '.  Thus,  the  elastic  or  adiabatic 

W >ulll)WH»H  h 

Mt»i  ■ >,»  ii  * I0n  «c  '.  (2.16) 


l 
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Electronic,  vibrational  and  some  rotational  energy  level  spacings  are  sufficiently  large  to  ensure 
that  the  adiabatic  approximation  (2.16)  is  valid. 

The  atomic  system  is  described  by  density  matrix  elements  and  we  seek  a transport  equation 
describing  their  time  evolution.  Before  a collision,  the  wave  function  is  given  by 

MR,r)  = + A2(R,t-)<l>2(r)l  (2.17) 


where  R is  the  center-of-mass  coordinate,  r the  electronic  coordinate,  A,\R,  f)  the  state  i probability 
amplitude  and  i/>,(r)  th.  state  i eigenfunction.  Following  the  elastic  collision,  the  wave  function 
is  given 


[e<i /*)/>'« 

e<i RAl(R,t+)+—£—f1(v'  -+  v)S,(/?,f  + ) Mr,(r) 

[e0/*)P'R  "1 

e(>»>p'  *A2(R,t  + ) + — ^ — /2(v'  - r)S2(/?,  f + ) U2(r), 


(2.18) 


where  /!,(/?,  r+)  and  S,(R,t+)  are  obvious  generalizations  of  (2.7a)  and  (2.7b),  respectively.  The 
change  in  the  density  matrix  element  ptj  resulting  from  a collision  is 

[e<l/*)p'K  ~1 

cm'*  '*AtR,t+)  + -j-fA*'  - »0W  f+) 

e(i  lh)pR  "1* 

x e(*/«p  * + + ) - AlR,C)AfiR,C)*.  (2.19) 


Following  the  same  procedure  that  was  used  in  going  from  (2.8)  to  (2.12),  one  arrives  at  the  quan- 
tum-mechanical transport  equation 

dpiji f,  t)/dt lei,  = -Tj /Wpi/r,  t)  + |d v'W^v'  -*  v)piJ{v\  t),  (2.20) 

Wijiv'  -*  k)  = . Vim/p)*  j*d Vp  d vr >Vp( Vp)^- 1 <5 [ »*r  + - ( m/p)v  + r;] 

x 6(vt  - V')f(v't  -*  v,)ft( -*  vr)*,  (2.21) 

r(/r)  = \JT  Jdrpl*yrp)»f^^  [fUr  - rr)  - //*,  - r,)*].  (2.22) 

The  contributions  to  dp,/r,  t)/dt|con  are  easily  understood  in  terms  of  quantum-mechanical 
scattering  processes.  The  two  terms  appearing  in  (2.22)  correspond  to  the  interference  of  the 
unscattered  j or  i wave  with  the  spherically  scattered  i or  j wave,  respectively.  As  in  eq.  (2.10), 
these  terms  contribute  only  in  the  forward  direction.  The  contribution  (2.21)  represents  the  overlap 
of  the  i and  j spherically  scattered  waves.  The  collision  “kernel”  (2.21)  and  “rate”  (2.22)  are  now 
complex  quantities  so  that  a simple  classical  interpretation  of  the  transport  equation  is  no  longer 
possible.  The  transport  equation  is  best  discussed  by  considering  diagonal  and  ofT-diagonal  ele- 
ments separately. 


no 
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1.  Diagonal  elements  i = j.  For  i = j,  the  transport  equation  is 


<>i,(  »\  t)/ct  U = -rfi(r)pi{(i\r)  + j'dv'tt'jv'  - v)pn(v',t), 

(2.23) 

with  a real  collision  kernel 

Wjjfv'  ->  v)  = . v j*d»pd»>rWp(Vp)pr"  ‘^»-r  + ~ v'  ~ ~~  v + •’pj 

* d(rr  - v'r)\flv'r  -*  vr)|2 

(2.24) 

and  a real  rate 

r„(»’)  = JT  j'd»-pITp(»-p)tr(T1(»>r)  = Jdv'Hjjr  -*  i>'). 

(2.25) 

The  first  part  of  eq.  (2.25)  follows  from  eq.  (2.22)  and  the  optical  theorem,  while  the  second  part  of 
eq.  (2.25)  may  be  verified  by  direct  integration  of  eq.  (2.24).  Thus,  each  population  density  ph{ v,  t ) 
obeys  a transport  equation  equivalent  to  that  for  the  one-level  problem.  Since  the  kernel  and  rate 
are  state  dependent,  the  total  population  density  £iPji(*\  t)  does  not  obey  a transport  equation. 
In  this  sense  one  can  not  obtain  a classical  picture  for  scattering  of  the  atom  as  a whole.  Rather, 
it  is  possible  to  envision  a classical  limit  for  the  scattering  in  each  of  the  levels  i.  (This  result  is 
analogous  to  that  of  a Stern  Gerlach  experiment  where  one  may  assign  classical  trajectories  to 
each  of  the  spin  states.) 

Additional  simplifications  of  eqs.  (2.23  2.24)  are  sometimes  possible.  If  the  ratio  of  perturber 
to  active  atom  mass  is  small  {mjm  <$  I),  collisions  produce  a relatively  small  change  in  the  active 
atom's  velocity.  Mathematically,  this  situation  leads  to  a kernel  Hj((»-'  -►  »>)  that  is  sharply  peaked 
about  v'  = v and  enables  one  to  expand  pu(v'.  f)  appearing  in  eq.  (2.23)  about  r'  = v.  The  limiting 
form  of  eq.  (2.23)  for  this  weak  collision  case  is  a Fokker-Planck  equation 


dp,^v,t)/dt\eM  * [ — r„( v)  + Aoivflpaiv,  t)  + + A2{v)~  P^i' 


(2.26) 


with 


n! 


dr'(r'  — p)"Hj1(»’'  -»  r). 


(2.27) 


The  Fokker-Planck  limit  of  the  transport  equation  is  valid  only  if  pj r,  /)  has  a velocity  width  less 
than  the  kernel  width,  justifying  the  expansion  of  pj  v\  t)  about  r'  = »>.* 

The  opposite  limit,  mjm  P t.  corresponds  to  scattering  by  stationary  perturbers  with  the  kernel 
and  rale  given  by  eqs.  (2.13)  and  (2.14),  respectively. 

2.  Off-diagonal  elements  i * j.  A picture  of  the  scattering  process  is  shown  in  fig.  4.  As  in  a 
Stem-Gerlach  experiment,  one  can  associate  a trajectory  with  pu  and  p22  for  an  active  atom 
that  enters  the  collision  in  a linear  superposition  state.  However,  associating  a classical  trajectory 


• For  example,  one  may  have  a case  of  weak  collisions  Avjv  < I in  which  a radiation  field  excites  a very  narrow  velocity  distribu- 
tion Ac,.  If  Arf  « Apt,  the  Fokker  Planck  equation  is  not  valid,  even  though  the  collisions  arc  weak. 
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Fig.  4.  An  attempt  to  give  a classical  picture  of  a collision  between  a two-level  active  atom  A and  a perturber  P.  Before  the  collision, 
the  atom  is  in  a linear  superposition  of  states  I and  2.  Following  the  collision,  it  is  possible  to  assign  classical  trajectories  to  p, , and 
p22  but  not  to  pl2.  The  quantum  mechanical  picture  of  this  collision  is  again  given  by  fig.  2,  with  an  unscattered  and  a scattered  com- 
ponent for  each  state  amplitude. 


with  pi2  is  impossible  just  as  it  is  impossible  in  the  Stern-Gerlach  case.  A quantum-mechanical 
calculation  of  the  scattering  process  is  necessary  and  it  is  such  a calculation  that  led  to  the  quantum- 
mechanical  transport  equation  (2.20). 

The  need  for  quantization  of  the  atomic  center-of-mass  motion  escaped  notice  for  many  years. 
The  development  oflasers  and  laser  spectroscopy  served  as  an  impetus  fora  reevaluation  of  theories 
of  pressure  effects  in  atomic  systems  and  led  to  the  realization  that  a classical  description  of  the 
atomic  motion  was  inadequate.  1 remember  when  I came  upon  this  idea  in  1969  and  discussed 
it  with  Professor  Lamb.  1 was  having  some  trouble  formulating  the  problem  quantum-mechanically 
and  he  suggested  a reference  to  a 1930  article  by  Bethe  [3].  1 mention  this  incident  because  it  is 
typical  of  Professor  Lamb's  ability  to  pinpoint  appropriate  references  to  the  literature,  even  though 
the  article  may  be  40  years  old!  I was  able  to  generalize  Bethe’s  method  to  solve  the  dilemma  of 
what  to  do  with  the  off-diagonal  density  matrix  elements.  The  two  articles  emerging  from  that 
work  [4.  5],  while  not  presenting  the  results  in  the  form  of  a transport  equation,  were  straight- 
forward. quantum-mechanical  calculations  of  the  type  that  Professor  Lamb  favors.  I believe  that 
he  was  pleased  with  these  papers,  owing  to  the  basic  nature  oflhe  calculations  and  the  consequences 
of  the  results. 

Although  eqs.  (2.20-2.22)  are  difficult  to  interpret  owing  to  the  presence  of  a complex  kernel 
and  rate,  several  limiting  cases  may  be  discussed: 

Equal  collision  interactions.  In  the  limit  of  equal  scattering  amplitudes  for  levels  i and  j, 
f( -*  rr)  = //•’'  -»  vr),  the  state  independent  kernel  fT( v'  -»  »•)  and  rate  F(r)  are  real  and  given 
by  eqs.  (2.24)  and  (2.25),  respectively.  With  regards  to  collisions,  the  atom  has  no  internal  structure 
since  the  collision  interaction  is  state  independent.  In  this  case,  one  can  assign  a velocity  to  the 
atom  as  a whole  rather  than  to  individual  density  matrix  elements.  The  only  effect  of  collisions 
on  off-diagonal  density  matrix  elements  is  to  change  the  velocity  associated  with  those  elements. 
This  velocity  change  leads  to  a subsequent  change  in  the  atom's  Doppler  phase  factor  and  to 
narrowing  of  spectral  profiles  in  linear  spectroscopy  (clearly,  it  is  not  appropriate  to  speak  of  pres- 
sure broadening  in  this  limit)  and  to  either  broadening  or  narrowing  in  non-linear  spectroscopic 
profiles.  The  limit  of  state-independent  scattering  may  be  found  if  le  els  l and  2 belong  to  the  same 
hyperfine  multiplet.  fine  structure  multiplet,  vibrational  or  electronic  multiplet. 

Light  perturbers  - mjm  4 1.  In  this  limit,  the  active  atom's  velocity  is  unchanged  in  collisions, 
to  a first  approximation.  The  expression  for  the  kernel  (2.21)  reduces  to 


Wjji  *'  -*  r)  = A'6(v 


- o J. 


d*-pdfi,,i:rI Vp{vp)f(V'  - *•;)/,<»- r -*  <)* 


(2.28) 


leading  to  a corresponding  change  of  eqs.  42  20  2 22) 

Sp,^f,t)/St\coH  = -y 


I 


(2.29) 
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where 


?i/r)  = Jdt-pJVp(»>p)t)r  ^2nh/ipvt)  [f(vt  - vr)  - f/v,  ->  »-r)*] 

- Jd^/iK  -♦  v'r)f{vr  - r;)*J.  (2.30) 

The  quantity  yt/v)  is  a complex  decay  parameter  that  appears  in  the  quantum-mechanical  pressure 
broadening  theory  of  Baranger  [6],  in  which  all  active  atom  velocity  changes  are  neglected. 

To  interpret  eq.  (2.29)  one  need  only  consider  the  effect  of  a collision  on  p i;-  when  active  atom 
velocity  changes  are  neglected.  For  an  elastic  collision,  only  the  phase  of  the  atomic  state  probability 
amplitudes  can  change,  leading  to  a change  in  p(j  for  a single  collision 

Spij  = (eie  - \)PiJ,  (2.31) 

where  0 is  a relative  phase  change  (0,  - l)j)  of  the  probability  amplitudes.  If  this  equation  is  aver- 
aged over  collision  histories  one  arrives  at  (2.29)  with 


W 


ao 

= ^ j*dvp  Wp(vp)t>r  j* 


dvpWp(vp)i;r  j 2nbdb[]  - ew'"'], 


(2.32) 


where  b is  the  collision  impact  parameter.  The  correspondence  between  eqs.  (2.30)  and  (2.32)  is 
obtained  if  the  scattering  amplitudes  in  eq.  (2.30)  are  expanded  in  partial  waves  and  the  partial 
wave  phase  shifts  are  evaluated  in  the  WKB  or  eikonal  approximations. 

With  the  neglect  of  active-atom  velocity  changes,  collisions  are  solely  phase-interrupting  in 
their  effect  on  off-diagonal  density  matrix  elements.  Each  collision  produces  an  abrupt  phase 
change  for  theatomicoscillatorand  this  effect  is  translated  into  a shift  and  broadening  in  the  spectral 
profile  associated  with  the  transition.  Traditionally,  one  has  used  an  equation  of  the  form  (2.29) 
to  describe  collision  effects  on  spectral  profiles;  it  is  understandable  that  one  talks  of  pressure 
broadening  in  this  limit. 

Scattering  in  one  state  only.  If  the  collision  interaction  is  much  stronger  in  one  of  the  levels 
than  the  other  (|/,(r'  -►  vr)|  1/2(1-;  -►  rr)|  or  | f2(v'r  -*  vr)|  |/,(r;  -►  rr)|),  the  scattering  inter- 
action in  the  weakly  interacting  state  can  be  neglected  to  first  approximation.  In  that  limit. 


W^v  -*  r)  « 0 and  (2  20)  becomes 

Icon  = -f./Wp./M), 

(2.33) 

with 

W = i^|drpH;(rp)pr(g)  [/(„,  - yr)Sit  -f/vr  - ,r)*<*,J 

(2.34) 

and  x indicating  the  strongly  interacting  level. 

Eq.  (2.33)  has  the  same  form  as  (2.29),  implying  that  collisions  are  phase-interrupting  in  their 
effect  on  p,/r,  t ).  However  velocity-changing  effects  are  absent  in  eqs.  (2.29)  and  (2.33)  for  different 
reasons.  In  eq.  (2.29),  the  presence  of  light  perturbers  led  to  negligible  active  atom  velocity  changes 
in  collisions.  In  eq.  (2.33)  only  the  forward  scattering  direction  contributed  to  the  change  in  p./v,  t) 
owing  to  quantum-mechanical  interference  effects.  Whereas  active  atom  velocity  changes  will  be 
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small  for  the  atom  as  a whole  when  eq.  (2.29)  is  valid,  it  is  possible  to  have  significant  velocity 
changes  associated  with  the  population  p„(v,  f)  of  the  strongly  scattered  x state  when  eq.  (2.33) 
is  applicable. 

One  might  expect  a strong  state  dependence  of  the  scattering  interaction  if  levels  1 and  2 corres- 
pond to  different  electronic  or  possibly  to  different  vibrational  states. 

Although  the  origin  of  eqs.  (2.29)  and  (2.33)  are  different  as  stressed  above,  they  both  lead  to 
the  same  type  of  equation,  which  can  be  interpreted  in  terms  of  an  effective  phase  interruption 
following  each  collision.  In  this  spirit,  I shall  write 

dPi,iv'  O/r’t  Icon  = (2.36) 


= W 


mjm  <Z  1, 

single-state-scattering. 


to  represent  both  cases  where  velocity-changing  collision  effects  on  ptJ{  v,  r)  are  absent.  The  para- 
meter yf/jv)  is  the  phase-interrupting  collision  parameter. 

Further  approximations.  It  may  be  possible  to  carry  out  perturbation  theory  with  one  of  the 
above  limits  as  the  zeroth  order  approximation.  For  example,  if  |/Kv'  -*  vt)|  > |/j(r'  -»  v,)|, 
a solution  to  first  order  in  /,<?'  -»  vr)  might  be  attempted.  On  the  other  hand,  if f * fj,  collisions 
are  velocity-changing  in  first  approximation  and  pu  follows  a prescribed  trajectory;  the  slight 
difference  in  f and  f produces  a correlated  phase  interruption  calculated  along  this  trajectory 
[7],  In  still  another  scheme  [8],  larger  impact  parameter  collisions  are  assumed  to  produce  small 
velocity  changes  while  close-in  collisions  are  taken  to  be  phase-interrupting  in  nature. 

General  case.  If  none  of  the  above  limits  are  realized,  the  “kerner  and  “rate"  are  given  by  (2.21) 
and  (2.22).  respectively.  Little  progress  has  been  made  in  their  evaluation,  although  one  might 
imagine  that  some  semi-classical  or  eikonal  approximations  could  be  tried.  A simple  interpretation 
in  terms  of  a velocity  change  and  a phase-interruption  occurring  in  each  collision  does  not  seem 
appropriate  since  the  two  effects  are  completely  intermingled.  One  would  expect  asymmetric  and 
distorted  profiles  when  both  velocity-changing  and  phase-interrupting  effects  are  important. 

2.3.  Impact  approximation 

The  quantum  mechanical  transport  equation  (2.20)  provides  a description  of  elastic  scattering 
for  any  n level  atomic  or  molecular  system.  This  equation  was  derived  on  the  implicit  assumption 
that  any  external  fields  which  are  interacting  with  the  system  do  not  appreciably  alter  it  during 
the  collision  duration  rc  % 10" 12  sec.  In  other  words,  collisions  occur  “instantaneously”  (impact 
approximation)  with  regards  to  other  time  scales  in  the  problem. 

The  impact  approximation  will  be  valid  (1)  if  the  atom  does  not  decay  during  a collision  (yrc  < 1, 
where  y is  some  decay  rate),  (2)  if  two  collisions  do  not  occur  during  a collision  time  (binary  collision 
approximation),  (3)  if  the  field  is  not  strong  enough  to  induce  transitions  in  a time  xc  (valid  for  field 
strengths  $ 107  W/cm2),  and  (4)  if  the  atom-field  detuning  A satisfies  |A|rc  « 1 (the  atom-field 
interaction  effectively  takes  place  on  a time  scale  rd  $ | A|"‘,  as  required  by  the  uncertainty  prin- 
ciple - one  must  have  rd  > rc  for  the  impact  approximation  to  be  valid). 

The  calculations  are  restricted  to  situations  where  the  impact  approximation  is  valid.  For  the 
case  of  a dilute  gas  interacting  with  weak  external  fields  considered  below,  the  only  restriction 
from  the  impact  approximation  is  that  the  detunings  | A|  must  be  <:  1012  sec  ~ 
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3.  Linear  spectroscopy 

Both  linear  and  non-linear  spectroscopy  can  be  used  to  probe  collision  effects  in  atomic  and 
molecular  systems.  In  section  4.  the  manner  in  which  non-linear  spectroscopy  can  be  used  to 
provide  differential  as  well  as  total  cross-section  data  is  explored.  In  this  section,  it  is  shown  that 
values  for  total  collision  cross  sections  may  be  extracted  from  linear  spectroscopy  line  profiles. 

In  linear  spectroscopy,  one  measures  the  absorption  or  emission  profile  associated  with  a given 
transition.  The  two-level  system  shown  in  fig.  5 serves  to  illustrate  the  essential  features  of  linear 
spectroscopy.  Atoms  are  incoherently  pumped  into  level  1 with  a rate  density  k(v)  assumed  to  be 
independent  of  position.  The  population  in  levels  1 and  2 decay  with  rates  y,  and  y2 , respectively, 
owing  to  spontaneous  emission.  A monochromatic  field  (propagation  vector  k = kZ.  frequency 
fl  = ck) 

E(R.t)  - »£cos(/c2  - Qt)  (3.1) 

is  applied  to  the  atom  and  drives  the  1-2  transition  which  has  a natural  frequency  oj.  The  population 
of  level  2 is  a measure  of  the  absorption  of  the  external  field  by  the  atom.  For  weak  external  fields, 
the  population  />22  as  a function  of  ft  yields  the  linear  absorption  spectral  profile. 


E vw— * 

a 

i 


a 


r. 

Az  £l~u> 


Fig.  5.  The  two-level  system  considered  in  this  work. 


To  calculate  p12,  one  finds  its  equation  of  motion  in  the  absence  of  collisions  and  adds  to  that 
the  collisional  contribution.  The  results  depend  critically  on  the  velocity  distribution  in  level  1 
which,  in  turn,  is  a function  of  the  excitation  rate  density  /.(►’).  To  simplify  the  discussion,  I shall 
assume  that  level  1 is  a ground  state  with  a thermal  velocity  distribution.  This  assumption  corres- 
ponds to  taking 


/.(*)  ~ 0,  y,  ~ o,  /( v)/v , ~ N(v). 

(3.2) 

with 

N{r)  = 

(3.3) 

IV0(r)  = (tim2)'3/2  e-1'2'"2. 

(3.4) 

where  u is  the  most  probable  speed  of  the  thermal  distribution 
of  the  theory  that  arise  when  the  pumping  density  A(r)  can  not 
equilibrium  velocity  distribution  (as  might  occur  in  emission 
are  discussed  elsewhere  [1]. 

Wo( m).  The  necessary  modifications 
be  satisfactorily  represented  by  an 
or  absorption  from  excited  states) 

3.1.  Line  shape  - no  collisions 

In  the  absence  of  collisions,  the  Hamiltonian  for  the  nth  atom  (neglecting  pumping  and  decay 
terms  which  will  be  added  later  on)  is 
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H(r„,  Rn . 0 = H0(r„ ) - (h2/2m)V2n  + V(r„.  R„,  t),  (3.5) 

where  r„  represents  all  electronic  coordinates  of  atom  n . /?„  is  the  center-of-mass  coordinate  of 
atom  n,  H0(r„)  is  the  free  atom  electronic  Hamiltonian  and  K(rH.  R„ , t)  is  the  atom-field  interation 

V(rn,  Rn , r)  = -/in  • £(/?„,  f),  (3.6) 

where  /i„  is  the  atomic  dipole  moment  operator.  The  wave  function  of  the  nth  atom  can  be  expanded 
in  terms  of  eigenfunctions  i/>,(rn)  of  H0(O  as 

= (3.7) 

I 

Density  matrix  elements  (diagonal  in  coordinate  space)  are  defined  by 

pUR„,t)  = A^Rn.t)A-j(R„,t)\  (3.8) 

which  obey  an  equation  of  motion 

dp",jlR„  m = -i ^jrilR',  t ) - V • r^Rt,  t)  + (i hr'[V(Rn,  t ),  p'(Rn,  0]y,  (3.9) 

where 

toij  = (E,  - Ej)/h;  t o = (02t.  (3.10) 

£,  is  the  state  i eigenenergy, 

f!j(Rn,  t)  = [h/(2mi)]  [AjlR„.  t)*\AURm,  t ) - ylR*.,  t)V/(J(E„.  f)*]  (3.1 1 ) 

is  a quantum-mechanical  current,  and  matrix  elements  of  V(Rn.  t)  are  given  by 

0 = (*dflli/'l(r1I)*K(r11,  Rn , f)iA,(rn)  = (1  - <5„)//0£  cos  (kZ  - ftf),  (3.12) 


with  = pjj  a matrix  element  of  the  x component  of  the  dipole  moment  operator. 

Provided  that  the  external  field  does  not  alter  the  atomic  center-of-mass  motion,  it  is  possible 
to  obtain  a classical  limit  for  eq.  (3.9)  (classical  with  regards  to  the  center-of-mass  motion  - the 
internal  states  are  always  treated  quantum-mechanically).  An  external  field  does  change  the  atomic 
velocity  slightly  owing  to  recoil  effects  resulting  from  the  emission  or  absorption  of  radiation. 
If  these  recoil  effects  are  neglected,  however,  one  can  associate  a classical  position  R„  and  velocity 
vK  to  the  average  position  and  velocity,  respectively,  of  a wave  packet  representing  the  atom. 
In  this  limit,  the  current  becomes 

■/?(!?„,  r)  — > • \p1j{R„,  v„,  t)  (3.13) 

and  eq.  (3.9)  goes  over  into 

ty'JR*.  V 0/df  = -i lOupl^R',  v„,  f)  - • VplfiRv  vn,  f)  + (\h)~  ‘[T(E„,  f),  p”(R„  V f)](J, 

(3.14) 

which  now  possesses  a characteristic  convective  flow  term  - »>B  • \p1j. 

Macroscopic  density  matrix  elements 

ph{R.  r,  f)  = X I &vnp1frRn.  ME  - Rn)S{v  - r„)  = X PUR'  0 (315) 

n J n 


.... 
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are  defined  such  that 

N(*,M)  = f P«(*,M)  (3.16) 

is  the  density  of  active  atoms  in  the  sample.  Density  matrix  elements  ptj{R,  v,  r)  also  satisfy  (3.14) 
since  each  atom  interacts  independently  with  the  field.  The  pumping  of  level  1 can  be  introduced 
by  the  addition  of  a term  2(v)<5(1<5,,  and  the  decay  by  the  addition  of  a term  — y.yP,//?,,.  r), 

Vij  = (7i  + y;)/2,  (3.17) 

to  the  rhs  of  eq.  (3.14).  The  resultant  equation  for  the  macroscopic  density  matrix  elements  is 

dpit<R,  v,  t)/dt  = - v ■ \pit{R , v,  f)  - i OijpijiR,  v , r)  - y^/R,  v , f)  + (i/i)-  ‘[(^(rt,  t),p(R,  v,  f)],j 

+ HWntjn  (3.18) 

which  describes  the  time  evolution  of  the  system  in  the  absence  of  collisions.  A term  corresponding 
to  the  spontaneous  emission  from  level  2 serving  to  populate  level  1 could  be  added  to  (3.18); 
however,  this  process  does  not  lead  to  any  modification  of  the  line  shapes  in  linear  spectroscopy. 

Although  transient  solutions  ofeqs.  (3. 18)  are  sometimes  useful  for  collision  studies,  only  steady- 
state  solutions  are  considered  in  this  work.  Introducing  the  field  interaction  representation 


p12(/f,  v,  t)  = p12(v)exp  [ — i(kZ  - JJf)].  (3.19a) 

p^R,  r,  r)  = pti{v),  (3.19b) 

employing  the  rotating  wave  approximation  (neglect  of  rapidly  varying  anti-resonance  terms), 
and  using  eq.  (3.12),  one  finds  that  the  steady-state  density  matrix  elements  p(j{v)  satisfy  the  follow- 
ing equations  for  the  two-level  system  of  fig.  5: 

yiPul*’)  = i*[P2i('’)  ~ PizW]  + ^(»0,  (3.20a) 

?2P22('')  = iz[Pi2<»’)  - ^2i(»’)].  (3.20b) 

Oh 2 - ik^)pi2(»’)  = “ Pi !(*)].  (3.20c) 

P2i(r)  = Pi2(»’)*.  (3.20d) 

where 

X = PtiE/lh,  (3.21) 

flu  = 7 1 2 + 'A.  (3.22) 

and  the  detuning  A is  given  by 

A = £2  - to.  (3.23) 


Eq.  (3.20)  may  be  solved  exactly,  but  a perturbation  solution  of  the  equations  is  all  that  is  needed 
in  linear  spectroscopy.  The  two  perturbation  chains  leading  to  a value  of  p22(»>)  to  second  order  in 
X are 


Pi  l P 12  Pl2< 


I 


Pi  1 P 21  P 22- 


(3.24a) 

(3.24b) 


M 
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These  chains  are  conveniently  pictured  in  the  Feynman  type  diagrams  of  fig.  6.* 

Each  horizontal  line  represents  a history  for  a state  amplitude,  while  density  matrix  elements 
are  formed  from  a product  of  the  amplitudes.  The  wavy  lines  represent  atom-field  interactions 
and  the  density  matrix  elements  for  each  part  of  a history  are  indicated  explicitly  between  the 
two  lines.  Figures  6a  and  6b  correspond  to  chains  (3.24a)  and  (3.24b),  respectively,  but  it  is  easy 
to  see  that  they  are  really  the  same  diagram  with  a different  “time-ordering”  for  the  fields  (since 
the  problem  is  time  independent,  the  words  “time-ordering”  are  used  solely  in  analogy  with 
Feynman  diagrams).  Diagrams  of  this  type  are  especially  useful  when  collisions  are  incorporated 
into  the  problem. 


p ' p i p 

'll  • rl  t • rlJ 


Fig.  6.  Diagrams  illustrating  linear  absorption  of  a field  £ (*  is  the  field  strength  in  frequency  units)  talcing  an  atom  from  state  1 to 
state  2.  Each  solid  line  represents  a state  amplitude  and  each  wavy  line  an  atom-field  interaction.  The  numbers  above  or  below  the 
lines  indicate  the  state  the  atom  is  in  and  the  product  of  the  two  amplitudes  gives  rise  to  the  density  matrix  elements  indicated  in  the 
diagrams.  Figures  (a)  and  (b)  correspond  to  chains  (3.24a)  and  (3.24b).  respectively. 


The  perturbation  solution  of  (3.20)  according  to  the  chains  (3.24)  is  straightforward  and  one 
obtains 

Pii(v)  = 2x2NWo(r)y12{y2[(yi;!)2  + (A  - fctf.)2]}-1,  (3  25) 

with  W^r)  given  by  (3.4).  The  absorption  line  shape  is  defined  as 


/(A)  = I dr p22(v)  = ^ 


00 

--  f 

2)3/2  J 


)’i2  exp  [-(d/m)2] 
(y12)2  4-  (A  - kv.)2 


and  is  the  convolution  of  the  Doppler  shifted  Lorentzian  absorption  profile  associated  with  atoms 
having  velocity  r with  the  Gaussian  distribution  of  atomic  velocities  )V0(r).  The  line  shape  (3.26) 
is  the  so-called  Voigt  profile  and  may  be  expressed  in  terms  of  the  plasma  dispersion  function 
Z(/4)  as 

i/av  _ 2)f2W  .7  (&  + iyi2\  t\-n\ 


/(A  ) = -±—Zi 
y2ku 


M + »y  i 

V *«  / 


where 


Z(n)  = Z,(pi)  + iZj(/i) 


00 

■ai 


H + x 


• A similar  idea  was  presented  earlier  [4,  5]  with  another  pictorial  format.  Professor  S.  Stenholm  suggested  the  more  compact 
representation  of  fig.  6. 
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with  Im(p)  > 0.  Since  y12  is  generally  much  less  than  the  Doppler  width*  fcu[ln2]l/2,  /(A)  is 
essentially  a Gaussian  function,  representing  the  well  known  “Doppler  broadening”  of  spectral 
profiles. 


3.2.  Collisions  - general  considerations 


In  the  impact  approximation,  the  external  fields  and  collisions  contribute  independently  to 
dptJ(R,  v,  t)/dt  allowing  one  to  simply  add  eqs.  (3.18)  and  (2.20).  The  resulting  quantum  mechanical 
transport  equation 

dp^R,  v,  t)/8t  = -v  S/p^R,  v,  t)  - ioiijPtjlR,  v.  t)  - y^p^R.  v,  t)  + (ift)~  *[K(/?,  t),  p(R , v,  f)](J 


- lyvjp,/*,  v,  t)  + Jdv’ W^v'  -*  v)piJ{R , v\  t)  + }{v)SnSjt  (3.29) 


characterizes  atoms  interacting  with  external  fields  while  undergoing  collisions.  The  contribution 
(2.20)  gives  the  average  effect  of  collisions  that  “instantaneously”  change  the  velocity  and/or  phase 
associated  with  density  matrix  elements.  In  contrast,  eq.  (3.18)  gives  the  change  in  p{J{R , v,  t)  for 
the  time  between  collisions,  in  which  the  atomic  center-of-mass  motion  associated  with  each  density 
matrix  element  is  that  of  a free  particle. 

Steady-state  equations  for  the  p^v)  defined  by  eq.  (3.19)  may  be  obtained  from  (3.29)  in  the 
rotating  wave  approximation  and  are  given  by 


01  + rn(r)]pu(r)  - Jdv'W', ,(»'  -*  v)pn(i>')  = i*[p2,(r)  - p,2(W]  + Hv), 

02  + r22(r)]p22(r)  - Jdr'W22(v'  -*  v)p22(r')  = ix[p,2(v)  - p2, (»>)], 

[»fi 2 “ + r,2(i-)]p,2(»')  - j*dp'W/12(»’'  - v)p12( v')  = ix[p22(v)  ~ Pnf*)], 


(3.30a) 

(3.30b) 


(3.30c) 

(3.30d) 


P2 l(^)  = Pl2(*’)*- 

Solutions  of  these  linear  coupled  integral  equations  are  very  difficult  to  obtain  for  arbitrary 
collision  kernels  W^v'  -*  r). 

However,  a formal  solution  of  eqs.  (3.30)  can  always  be  found  that  turns  out  to  be  specially 
useful  for  perturbation  (in  the  external  field)  calculations.  This  solution  to  eqs.  (3.30)  is  given  by 


Pul')  = j*dr'C,1(r'  - v)[ix{p2,(r')  - p12(r')}  + -*(•’')], 
P22(*')  = |dr'G22(r'  -►  »)  [iz  {Pi2<0  “ P2i(r')}], 

Pi 2(^1  = fd»'G12(r'  -♦  •’)[ix{P22(»'')  ~ Pn(r')}], 

Pu(')  = Jdr'G2,(r'  - r)  [ix  {pn(r')  - p22(r')}]  = P,2(»’)*. 


(3.31) 


* By  "width".  I always  refer  to  the  half  width  at  half  maximum  (HWHM). 


•<* 
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where  the  propagators  G0<v'  -+  v)  satisfy 


!>,>*)  + r^jGijiv'  - v) 


- J* 


Wijiv"  - v)GiJW  -*  v")  = S(v  - v'). 


nM  = y »/i2(t'z)  = hi  - >kvz  = r/21(uz)*.  (3.33) 

All  the  effects  of  collisions  are  contained  in  the  propagator  equation  (3.32)  which,  in  general,  is 
very  difficult  to  solve. 

Assuming  one  has  a solution  of  eq.  (3.32)  for  G./v  -*  »>),  he  can  easily  obtain  an  iterative  solution 
of  eq.  (3.3 1 ).  The  perturbation  chains  leading  to  a value  for  p22(v) to  second  order  in  x are  completely 
analogous  to  those  of  (3.24).  One  obtains 


ft  2(*)  = X2  j 


dv2  dr,  dr0G22(r2  -►  r)G12(v,  -»  v2)G,  ,(r0  -►  v,)A(v0), 


corresponding  to  chain  (3.24a)  and 


dr2  dr,  dr0G22(r2 


(3.34a) 


(3.34b) 


corresponding  to  chain  (3.24b).  These  results  have  a simple  interpretation  in  terms  of  the  diagrams 
of  fig.  6.  For  example,  in  fig.  6a  and  in  (3.34a),  one  can  think  of  the  incoherent  pumping  as  providing 
a A(r0).  The  quantity  G,  ,(r0  -*  r,)  propagates  this  solution  up  to  the  first  field  interaction.  Between 
the  first  and  second  field  interaction  the  propagator  G,2(r,  -+  r2)  acts,  and  following  the  second 
field  interaction  the  propagator  G22(r2  -*  r)  carries  along  the  solution. 

Thus,  to  solve  (3.30)  to  any  order  in  perturbation  theory  in  the  external  field,  one  first  draws 
all  me  appropriate  diagrams  of  the  type  shown  in  fig.  6.  At  each  field  interaction,  there  is  a multi- 
plicative coupling  constant  ( ± i*)  that  enters  the  calculation.  Between  two  successive  field  inter- 
actions the  solution  is  propagated  by  the  G0{ v'  v)  that  corresponds  to  the  ptl{v)  which  appears 
in  the  diagram  between  those  field  interactions  (i.e.,  p,2(v)  appears  between  the  first  and  second 
field  interactions  in  fig.  6a).  Before  the  first  field  interaction,  G,  ,(r0  -»  v,)  propagates  the  incoherent 
pumping  A(r0)  and  following  the  last  field  interaction  G22(rn  -*  v)  propagates  the  solution  to  some 
final  velocity  v.  Integrals  over  all  intermediate  velocities  must  be  taken.  This  diagrammatic 
technique  is  used  in  section  4 to  obtain  results  for  a three-level  problem. 

By  adding  (3.34a)  to  (3.34b)  and  integrating  the  sum  over  all  v , one  determines  the  line  shape 
in  linear  spectroscopy 


= x2ytN  J*' 


dt»dr2  dv,  dr0G22(r2  -»  r)[G,2(r,  -»  v2)  + G,2(v,  -♦  v2)*]G,,(r0 


ViWVol 

(3.35) 


where  eqs.  (3.2  3.4)  have  been  used.  This  line  shape  may  be  simplified  in  two  ways.  First,  collisions 
do  not  change  IF0(r)  since  it  is  already  on  equilibrium  distribution  so  that 

jdr0G(((r0  ^ r,)WM»0)  = Wr0(«',)/?«-  (3  36) 

Second,  diagonal  propagators  have  the  property,  easily  derived  by  integrating  eq.  (3.32)  over 
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velocity,  that 


I 


dvG.ilv' 


»)  = l/y*. 


(3.37) 


Eq.  (3.37)  expresses  the  fact  that  the  total  population  of  any  level,  J drpil(»'),  is  unaffected  by  elastic 
velocity-changing  collisions  (collisions  redistribute  the  velocities  but  do  not  change  the  total 
population).  In  terms  of  the  diagrammatic  solutions  of  fig.  6,  eqs.  (3.36)  and  (3.37)  eliminate  the 
need  to  consider  collision  effects  before  the  first  field  interaction  (when  an  equilibrium  velocity 
distribution  existed)  and  following  the  last  field  interaction  (since  the  total  population  only  of  the 
final  state  is  measured). 

Using  eqs.  (3.36)  and  (3.37)  in  eq.  (3.35),  one  obtains 


/(A) 


Vi  J 


dvG^i-o 


•O^oK)  + conjugate, 


(3.38) 


which  is  the  final  form  for  the  line  shape  in  linear  spectroscopy.  Only  collision  effects  on  p12(v) 
enter  in  eq.  (3.38)  since  the  effects  of  collisions  on  populations  have  been  eliminated  by  eqs.  (3.36) 
and  (3.37).  Limiting  forms  for  eq.  (3.38)  are  discussed  below  corresponding  to  G12(v'  ->  v)  deter- 
mined by  collisions  that  are  either  phase-interrupting  or  velocity-changing  in  their  effect  on  p12(v). 


3.3.  Line  shape  - phase-interrupting  collisions 

When  collisions  are  phase-interrupting  in  nature,  the  collisional  time  rate  of  change  of  p12(i\  0 
is  given  by  eq.  (2.36).  Consequently,  the  collision  terms  in  eq.  (3.29)  for  p12(/f,  v,  t)  are  replaced 
by  — yf2(r)p12(^,  v,  r)  and  this  change  is  reflected  in  eq.  (3.32)  as 

0/i2(t>z)  + V?2(',)]Gf2(i''  -»  v)  = <5(v  - »>')•  (3.39) 


Recalling  the  definition  (3.33)  for  iu2(vz),  one  obtains  the  propagator 


G?2(v'  - v)  = 


fly  - Q 

[v 1 2 + r?h2(v)]  + i[A  - kvz  + S?>)]' 


where 

rT»  = Re  [y&V)] ; S?h(y)  = Im  [tf»] 


(3.40) 

(3.41) 


are  the  phase-interrupting  broadening  and  shift  parameters,  respectively.  The  propagator 
G,2(v'  -»  v)  is  proportional  to  flv  — v')  for  phase-interrupting  collisions  since  there  is  no  velocity 
change  associated  with  off-diagonal  density  matrix  elements  for  reasons  discussed  in  subsection 
2.2.2.  The  line  shape  obtained  by  combining  eqs.  (3.38)  and  (3.40)  is 


2x2N  f [Vi2  + rft(r)]e-^ 

y2(jru2)3/2  J [y12  + r?h2(r)]2  + [A  - kvz  + Sf2(v)]2' 


(3.42) 


The  absorption  profile  is  broadened  and  shifted  from  the  no-collision  line  shape  (3.26);  one  can 
go  from  eq.  (3.26)  to  (3.42)  by  the  simple  replacements 


-V,2  + r?h2(v);  A-A  + S?h2(v). 


V I 2 


(3.43) 
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The  broadening  and  shift  can  be  understood  in  terms  of  the  phase-interrupting  collision  mecha- 
nism. Collisions  interrupt  the  oscillator’s  coherence  train  and  result  in  an  effective  decrease  of  the 
oscillator  lifetime  with  a corresponding  broadening  of  the  absorption  profile.  The  persistence 
of  some  phase  memory  following  a collision  (for  collisions  with  large  impact  parameters)  results 
in  a shift  of  the  line  center.  That  distant  collisions  contribute  to  the  shift  is  easily  seen  from  eq. 
(2.31)  which,  when  expanded  for  small  phase  shifts  |0|  <?  1,  leads  to 

dpl2(R,  v,  t)/3t|coll  = i <d>  p12(Jf,  v,  f),  (3.44a) 

with 

00 

<0>  = ^JdrpW>p)„r  f 2itbdb9{b,vr).  (3.44b) 

^mtn 

Thus,  distant  collisions  lead  to  a shift,  (d-*m  + <$>,  of  the  oscillator’s  resonance  frequency. 
For  strong  collisions  (|0|  > 1),  the  average  value  of  p12(v,  t ) after  a collision  is  zero;  these  collisions 
give  rise  to  a decay  of  p,2(v,  t)  (and  to  a broadening  of  the  absorption  profile),  but  do  not  contribute 
to  the  shift. 

The  line  shape  (3.42)  has  some  interesting  features.  The  broadening  and  shift  parameters  are 
given  as  functions  of  velocity,  but,  typically,  they  will  be  functions  of  speed  only.  This  result  follows 
from  eqs.  (2.37,  2.30,  2.34)  provided  that  /]( rr  -*  vr)  depends  only  on  the  magnitude  of  vr  and  that 
IV  (vp)  is  an  even  function.  Since  these  limits  are  generally  realized,  it  is  assumed  that  Ti*h  and 
S?y  are  functions  of  v only.  1 shall  call  the  resulting  line  shape  (3  42)  a Speed  Dependent  Voigt 
Profile  (SDVP).  The  conditions  under  which  the  speed  dependence  is  significant  is  discussed  below. 

The  presence  of  a speed  dependent  shift  leads  to  an  asymmetry  of  the  profile.  To  see  how  this 
asymmetry  arises,  replace  (r2  4-  i>2)  by  some  average  u2  and  consider  the  contribution  to  the  line 
shape  from  atoms  moving  with  ±vZo.  This  contribution, 

yj^fTS  exP  { dvZ0 

[>’i2  + r?h2(u|  + t>,?0)] 

[Via  + nft(t4  + u20)]2  + [A  — kvZ0  + S?$(i4  + i’2„)]2 

[v  1 2 + rf2(«T +»?„)] 

[y  1 2 + r?h2(«|  + v2Zo)Y  + [A  + kvt0  + S?h2(u|  + r20)]2 

is  a symmetric  function  of  A about  -S?5(f20  + “t)-  If  one  next  considers  the  contribution  from 
atoms  moving  with  ±vZl , he  finds  it  to  be  symmetric  about  -S?2(t;2  + u£)  / -S?2(t>?0  + u$). 
Thus  the  sum  of  the  two  contributions  is  asymmetric.  Since  !V0(v)  is  composed  of  ± vz  pairs,  the 
line  shape  (3.42)  is  asymmetric.  Furthermore,  since  the  asymmetry  will  be  pressure  dependent, 
the  shift  of  the  line  center  will  be  a non-linear  function  of  perturber  pressure.  It  may  also  be  noted 
that  erroneous  values  for  the  broadening  parameter  will  be  obtained  if  one  attempts  to  fit  experi- 
mental data  to  a normal  Voigt  profile  rather  than  a SDVP. 

These  effects  are  important  only  if  the  speed  dependence  of  y?2  = T'Jy  + iSfy1  is  significant. 
If  mp  <|  m,  v,  * ip  and  it  follows  from  eqs.  (2.37,  2.30,  2.34)  that  y^h  depends  on  the  most  probable 
perturber  speed  but  not  on  the  active  atom  speed  v.  Therefore,  significant  deviations  from  a Voigt 
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profile  can  occur  only  for  mp  > m.  Since  the  condition  mp  > m is  not  consistent  with  the  phase- 
interrupting  collision  limit  (2.30)  requiring  mp  m,  collisions  involving  negligible  active  atom 
velocity  changes  can  lead  only  to  normal  Voigt  profiles.  On  the  other  hand,  deviations  from  a 
Voigt  profile  could  occur  in  the  quantum  interference  phase-interrupting  collision  limit  (2.34), 
provided  mp  > m.  Unless  mp  > m,  however,  it  is  generally  a good  approximation  to  neglect  the 
speed  dependence  of  the  shift  and  broadening  parameters  [9, 10]. 

A line  shape  analysis  using  (3.42)  enables  one  to  determine  r?2  and  S?2  which  are  related  to  the 
cross-section  and  real  part  of  the  forward  scattering  amplitudes  in  states  1 and  2.  If  one  state  is  more 
strongly  scattered  than  the  other  the  total  elastic  scattering  cross-section  for  that  state  may  be 
found.  In  this  manner,  linear  spectroscopy  profiles  can  yield  excited  state  total  elastic  cross  section 
data. 

The  line  width  (HWHM)  of  (3.42)  as  a function  of  perturber  pressure  is  shown  in  curve  a of 
fig.  7 for  y,  2 < ku , neglecting  the  speed  dependence  of  T?2  and  A?2.  At  low  pressures  the  width  is 
dominated  by  the  Doppler  contribution,  but  as  the  pressure  is  increased  to  the  point  where 
r?2  > ku,  the  line  profile  becomes  Loren tzian  with  width  r?2.  If  the  speed  dependence  of  the  line 
shape  parameters  is  ignored,  eq.  (3.42)  becomes 

»/*,  _ tfN  „/(A  + Sfh2)  + i(y,2  + ryh2)\ 

( ' V zku  \ ku  / 

3.4.  Line  shape  - velocity-changing  collisions 

In  contrast  to  phase-interrupting  collisions,  collisions  which  are  velocity-changing  in  their 
effect  on  p12(»0  give  rise  to  a narrowing  of  the  linear  absorption  or  emission  profiles  associated 


I 


PERTURBER  PRESSURE  (TORR) 


Fig.  7.  Typical  behavior  of  the  line  width  (HWHM)  as  a function  of  pressure  in  linear  absorption.  At  low  pressures  the  width  is  ap- 
proximately equal  to  the  Doppler  width  Jku[ln2],/2  (assuming  y,2  4 ku).  (a)  For  phase-interrupting  collisions  the  width  increases 
monotonically  with  pressure,  becoming  linear  with  pressure  for  r?j  > ku.  (b)  For  velocity-changing  collisions,  the  line  width  decreases 
monctonically  with  pressure,  ultimately  reachinga  value  of  yu.  (c)  When  levels  1 and  2 experience  nearly  equal  collision  interactions 
the  width  decreases  and  then  increases  with  increasing  pressure. 
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with  the  1-2  transition.  While  it  is  not  difficult  to  understand  the  narrowing  effect,  it  will  aid  the 
discussion  to  first  introduce  the  concept  of  an  effective  collision  rate. 

Imagine  that  a velocity  distribution  p„(v.  t)  is  produced  by  some  means.  The  velocity  distribution 
is  centered  about  v0  x u with  a spread  (Au0).  Collisions  will  cause  the  velocity  distribution  to 
relax  back  towards  a Maxwellian  distribution  centered  at  r = 0 with  a width  (HWHM)  of  u(ln  2),/2. 
The  rate  at  which  this  relaxation  occurs  is  called  the  effective  collision  rate  and  is  denoted  by  rj,. 
The  effective  collision  rate  can  be  much  smaller  than  the  velocity-changing  rate  rH,  especially 
for  weak  collisions.  If  the  rms  velocity  change  per  collision  is  (Au),  there  must  be  at  least  u/Au 
collisions  to  restore  the  Maxwellian  distribution.  A more  realistic  estimate  for  the  required  number 
of  collisions  turns  but  to  be  [u/(Au)]2,  taking  into  account  that  random  nature  of  the  collision 
process  (that  is,  <v>  undergoes  something  similar  to  a random  walk  from  its  initial  value  r0  to  its 
final  value  of  zero).* 

Thus,  the  relaxation  time  is  [u/(Au)]2/rl(  giving  an  effective  collision  rate 

r„  = ri([(Au)/u]2.  (3.46) 

For  reasons  to  be  discussed  below,  linear  absorption  experiments  are  sensitive  mainly  to  TJ, 
whereas  non-linear  spectroscopy  may  be  used  to  determine  both  T u and  (Au). 

The  spectral  profile  is  given  by  eq.  (3.38)  with  the  propagator  G”2(r'  -»  v)  satisfying  eq.  (3.32) 


(n  12  - kvJGW  ->  r)  + r\%(v)G\W 


- Jdv"W7|(v" 


r)Gr2(r'  - r")  = S(v  - r'). 


The  collision  kernel  W72(v'  -»  r)  and  rate  r^2(v)  are  real  and  state  independent  since  the  collision 
interaction  is  state  independent  when  collisions  are  velocity-changing  in  their  effect  on  pl2  (see 
discussion  in  subsection  2.2.2).  Eq.  (3.47)  is  difficult  to  solve  for  all  but  very  simple  collision  kernels. 
However,  one  can  explain  most  line  shape  features  without  explicit  reference  to  a particular  kernel. 

To  write  the  spectral  profile  (3.38)  in  a more  suggestive  fashion,  I first  rewrite  the  no-collision 
line  shape  (3.26)  as 


/( A)  = 


s-H*--" 


exp  { -y12t  - i(A  - kvz)t}  + conjugate. 


When  collisions  are  present,  one  might  expect  that  exp  (ik  ■ vt)  is  replaced  by  < exp  [it-  J{,  r(f')  dr']  ) 
where  v(0)  = r0  and  the  average  is  over  collision  histories.  With  this  clue,  one  writes  (3.38)  in  the 
form 


H 


dtC(t,  v0)e  <i,,j  + ia>*w0(vo)  + conjugate, 


where  C(t,  r0),  as  obtained  by  iteration  from  eq.  (3.47),  is  given  by 


• The  factor  (u/Au)]  is  easily  derived  for  the  case  of  weak  collisions  by  calculating  the  decay  rate  of  /»)  [11], 
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C(t,  v0)  = exp  { - r , 2(  v0)t } exp  (i* 


i Ac  • r0t)  + J* dr  j*df 


1 exp  { -nV^i } exp(i*  • VohWK'o  v) 


x exp{-nc2(i>)(f  - f,)}exp{i*r(f  - f,)}  + 


|d,dv,  T 


X exp(-nc2(v0)f,)exp(i*  - v0f,)H7|(v0 -►  v,)exp  {- nc2(v,)(t2  - rt)} 
x exp{i*  - v,(f2  - r,)}^?,  -»  r)exp{— n*2(r)(t  - f2)}  exp  {i*  • v(t  - r2)}  + ... 


= < exp 


ji*jV)dtj 


>. 


(3.50) 


Eqs.  (3.49,  3.50)  provide  the  expected  generalization  of  eq.  (3.48)  to  include  velocity-changing 
collision  effects. 

The  collisional  narrowing  effect  can  now  be  understood  as  a randomization  of  the  Doppler 
phase.  In  the  absence  of  collisions,  C(t , v0)  = exp  (iA  • v0t)  and  the  line  width  is  on  the  order  of  ku 
(assuming  ku  P y,2).  To  estimate  C(r,  r0)  in  the  presence  of  collisions,  it  is  useful  to  express  it  as 


C(t,  p0)  = <exp  [iA  • *(f)]> 


(3.51) 


where  R(t ) is  the  displacement  undergone  by  an  active  atom  in  time  t having  a velocity  r = r0  at 
t = 0.  Collisions  result  in  active  atoms  altering  their  velocity  by  a significant  amount  in  a time 
on  the  order  of  [(H  2)']  “ 1 (P  is  the  effective  collision  rate);  consequently,  a typical  atom  undergoes 
a type  of  random  walk  with  a mean  free  path  » u/(r,2)'.  The  mean  square  displacement  of  such  a 
particle  in  time  t is 


<[w)]*>  = [«/(rr2y]2[(rr2)'0  = (V/ori)']* 


(3.52) 


valid  for  (n2)'r  > 1.  The  collisionally  induced  reduction  of  the  Doppler  phase  from  its  no-collision 
value  k • r0f  to  a value  * A<[/?(r)]2>,/2  provides  a narrowing  mechanism  for  the  line  profile. 

One  can  now  trace  the  development  of  the  line  width  as  a function  of  pressure.  For  (r,2)'  < ku. 


the  atom  does  not  undergo  an  effective  collision  in  its  “coherence”  time  fcoh  = (ku)'1  and  the  line 


width  remains  close  to  the  Doppler  width  [tcoh  is  the  largest  time  for  which  the  integrand  in  eq.  (3.49) 


contributes  effectively  to  the  integral].  For  (f"2)'  > ku  but  (Au)2/(r”2)'  > y12,  rcoh  is  determined 


by  the  <fc2[K(fCOi.)]2>  = 1 giving  a value  equal  to  rcoh  = (r"2)'/(fcu)2  and  a line  width  of  order 

W = (ku)[lcu/(  rr2)'].  (3.53) 


The  Doppler  width  is  reduced  by  a factor  [(fcu)/(r”2)']  at  pressures  where  the  effective  mean  free 
path  [u/(r”2)']  is  less  than  the  atomic  wavelength  (2n/k).  At  still  higher  pressures  where 
<fc2[K(f)]2>  yl2f,  the  Doppler  phase  may  be  totally  neglected,  and  the  line  profile  is  a Lorentzian 
of  width  yI2. 

The  line  width  as  a function  of  pressure  is  shown  schematically  in  curve  b of  fig.  7.  At  low 
pressures,  the  line  is  a Voigt  profile  with  a width  equal  to  slightly  more  than  the  Doppler  width 
(assuming  y,2  « ku).  The  width  decreases  monotonically  with  increasing  pressure,  ultimately 
reaching  the  value  y,2.  In  the  infra-red  or  optical  range,  collisional  narrowing  becomes  important 
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at  pressures  on  the  order  of  several  hundred  Torr  and  has  been  experimentally  observed  on 
vibrational  and  rotational  transitions.  Inelastic  collisions  may  tend  to  obscure  the  narrowing 
effect  by  broadening  the  line  with  increasing  pressure. 

If  the  collision  interaction  for  levels  1 and  2 differed  slightly,  one  might  still  expect  collisional 
narrowing  over  a wide  range  of  pressures.  For  high  enough  pressures,  however,  (lcn)2/(r”2)'  < r?2, 
and  collisional  broadening  of  the  line  profile  may  be  seen.  This  case  is  shown  in  curve  c of  fig.  7. 

To  this  point,  I have  not  discussed  the  detailed  nature  of  the  collision  kernel.  While  the  line 
profile  does  depend  to  some  extent  on  the  specific  form  of  the  kernel  [12],  it  is  not  overly  sensitive 
to  either  rvl2  or  Am  separately,  but  rather  to  the  effective  collision  rate  (T"2)'  = rj2[(Au)/u]2. 
The  semiquantitative  discussion  of  the  narrowing  effect  given  made  no  reference  to  T\2  or  Am; 
only  (n2)'  was  mentioned.  It  should  be  noted,  however,  that  a physically  acceptable  kernel  (i.e., 
one  that  obeys  detailed  balancing)  should  be  used  if  one  wants  to  be  assured  of  getting  the  narrowing 
effect.  For  example,  a difference  kernel  W(v'  ->»>)=  W(v  — »>'),  commonly  used  to  describe  weak 
collisions,  does  not  obey  detailed  balancing  and  does  not  lead  to  collisional  narrowing.  From  the 
symmetry  of  such  a kernel,  it  is  easy  to  conclude  that  <»>>  = v0,  implying  that  equilibrium  is  never 
reached.  Consequently,  (rj2)'  = 0 and  the  effective  mean  free  path  is  infinite. 


4.  Non-linear  spectroscopy 

At  low  pressures  ($  100  Torr),  the  Doppler  effect  provides  the  major  contribution  to  the  line 
width  of  absorption  or  emission  profiles.  The  relatively  large  Doppler  width  tends  to  mask  any 
collision  induced  modifications  of  the  line  profile,  and  it  is  therefore  desirable  to  devise  schemes 
in  which  the  Doppler  effect  plays  a minimal  role  in  line  shape  formation.  Any  such  schemes  are 
also  well  suited  for  high  resolution  spectroscopic  studies.  In  the  past,  the  Doppler  effect  was 
reduced  by  using  either  cooled  bulb  samples  or  well-collimated  atomic  beams.  The  last  few  years, 
however,  have  witnessed  the  emergence  of  a large  number  of  "Doppler-free”  techniques  involving 
the  use  of  non-linear  spectroscopy. 

The  manner  in  which  non-linear  spectroscopy  can  be  used  to  eliminate  the  Doppler  width  is 
best  illustrated  by  two  specific  examples.  First,  consider  two-photon  transitions  between  atomic 
levels  separated  by  a frequency  tu31.  If  an  atom  is  moving  with  velocity  v and  absorbs  one  phi  on 
each  from  fields  £,  E having  propagation  vectors  k,  k\  respectively,  the  resonance  condition  for 
two-photon  absorption  is 

Cl  + O'  + (A  + A')r  = co3t.  (4.1) 

If  A « eq.  (4.1)  is  independent  of  the  atomic  velocity  and  each  atom  contributes  equally  in 
the  absorption  process.  In  effect,  the  Doppler  phase  created  by  the  first  field  is  exactly  compensated 
by  the  second.  In  this  Doppler-free  technique,  all  active  atoms,  regardless  of  velocity,  contribute 
in  line  shape  formation. 

As  a second  example,  consider  any  of  the  three-level  configurations  in  fig.  8.  Two  photon  transi- 
tions of  the  type  discussed  above  are  still  possible  for  these  level  schemes.  However,  an  additional 
resonant  interaction  is  possible.  If  the  field  E with  propagation  vector  * = kst  is  tuned  to  be 
nearly  resonant  with  the  1-2  transition  ir.  fig.  8,  only  those  atoms  for  which 

kvt  = Cl  — to  (4.2) 
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Fig.  8.  The  three-level  systems  considered  in  this  work:  (a)  inverse  V,  (b)  V,  (c)  upward  cascade. 


will  be  significantly  excited  by  the  field.  A probe  field  E with  propagation  vector  k = ±kt 
and  frequency  ft'  = k'c  will  be  resonant  with  the  2-3  transition  provided 

ft'  = w'  ± k'v„  (4.3) 

where  of  is  the  2-3  transition  frequency.  Since  v.  has  been  selected  by  the  field  E according  to  (4.2), 
the  resonance  condition  (4.3)  becomes 

ft'  = co'  ± (*'/*)  (ft  - a>).  (4.4) 


The  atomic  velocity  does  not  appear  in  eq.  (4.4),  indicating  that  a measurement  of  the  absorption  of 
E as  a function  of  ft  is  a Doppler-free  technique.  In  contrast  to  the  two-photon  Doppler  free  reso- 
nance condition  (4.1),  the  Doppler  free  nature  of  the  stepwise  absorption  of  fields  £ and  E is  due 
to  the  selective  absorption  of  the  field  £ by  a narrow  velocity  subset  of  atoms.  Only  a fraction  of 
the  atoms  are  used  in  this  Doppler-free  technique  whereas  all  the  atoms  are  used  in  the  two-photon 
Doppler  free  case  in  which  a cancellation  of  the  Doppler  phase  leads  to  the  Doppler-free  resonance. 

Thus,  non-linear  spectroscopy  offers  the  potential  for  achieving  Doppler-free  line  shapes. 
Moreover,  the  modifications  of  these  Doppler-free  line  shapes  resulting  from  collision  effects  can 
be  dramatic.  As  an  example  of  the  role  collisions  can  play,  reconsider  the  stepwise  absorption 
chain  (4.2H4.4).  Atoms  with  a selected  value  of  vz  given  by  eq.  (4.2)  are  excited  into  state  2 by  the 
field  £.  Collisions  result  in  a relaxation  of  this  narrow  velocity  distribution  towards  a thermal 
equilibrium  one.  The  absorption  profile  of  E measures  the  amount  of  velocity  relaxation  that  has 
occurred.  Since  it  is  collision  induced  velocity  changes  that  are  directly  monitored,  one  obtains 
information  on  the  differential  scattering  cross  section  in  level  2.  (More  precisely,  one  determines 
the  differential  scattering  cross  section  averaged  over  perturber  velocities  and  transverse  active 
atom  velocities.)  A detailed  discussion  of  this  effect,  as  well  as  the  collision  induced  line  broadening, 
narrowing,  shifting  and  line  shape  distortion  that  may  occur  in  non-linear  spectroscopy,  is  given 
below. 

To  illustrate  the  various  aspects  of  collision  effects  in  non-linear  spectroscopy,  I consider  the 
three-level  systems  of  fig.  8.  The  calculation  is  carried  out  for  the  inverted  V of  fig.  8a,  with  corres- 
ponding results  for  the  V (fig.  8b)  and  upward  cascade  (fig.  8c)  given  at  the  end  of  this  section. 
The  field  £,  pumping  A(r)  and  excitation  density  N(v)  are  the  same  as  for  the  two-level  system 
[eqs.  (3.1  X (3  2)  and  (3.3),  respectively].  The  field  £ drives  only  the  1-2  transition  and  a field 

E = /£'  cos  [ek'Z  — ftr)  (4.5) 

either  copropagating  (e  = + 1)  or  counterpropagating  (e  = - 1)  with  the  field  £ drives  only  the 
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2-3  transition.  The  population  of  level  3 is  monitored  as  a function  of  probe  frequency  ft'  for  fixed 
ft  and  serves  to  measure  the  probe  field  absorption  in  the  presence  of  the  “pump”  field  E.  The 
calculations  will  be  carried  out  to  lowest  order  in  the  fields  £ and  E since  it  is  in  the  weak  field 
limit  that  collision  effects  are  most  transparent. 

4.1.  Line  shape  - no  collisions 

In  the  absence  of  collisions,  eq.  (3.18)  gives  the  equation  of  motion  for  density  matrix  elements 


with  an  atom-field  interaction 

V(r,R,t)=  -p(E  + E).  (4.6) 

Steady  state  solutions  of  eq.  (3.18)  are  obtained  by  introducing  the  field  interaction  representation, 

p12(£,  v,  f)  = p12(r)  exp  [ — i(kZ  - Or)],  (4.7a) 

p23(R,  v,  t)  = p23(r)  exp  [i (ek'Z  - ft'f)],  (4.7b) 

p , 3(R,  v,  r)  = p , 3(r)  exp  { - i[(k  - ek')Z  - (ft  - ft')f ] },  (4.7c) 

PtAR.'rt)  = p,A  r),  (4.7d) 

making  the  rotating  wave  approximation,  and  recalling  that  E drives  only  the  1-2  transition  and  E 
only  the  2-3  transition.  When  eqs.  (4.7)  are  substituted  into  eq.  (3.18),  one  finds  that  the  steady 
state  density  matrix  elements  plV(r)  satisfy 

yiPnW  = 'X[P2l(i’)  - Pi2(*)]  + (4.8a) 

y2p22(*)  = •Z[pI2(*r)  - P2l(^)]  + «Z'[P32(»’)  ~ PliWl  (4.8b) 

y3 P33(»)  = - P32(')].  (4.8c) 

(>fi2  - ifci7*)Pi2<»r)  = ix[P22(»’)  - PnW]  - ix'PisM  (4.8d) 

(•hs  + i£k'vz)p23(v)  = i*'[p33(r)  - p22(r)]  + i*p,3(r),  (4.8e) 

[>13  “ Kk  - ek>,]Pi3(r)  = i*p23(r)  - ijt'p^M  (480 

PjM  = [Pi/*')]*.  (4.8g) 

where 

X = Pi2E/2h;  x'  = Pis^'/lh,  (4.9) 

flu  = Vi2  + iA,  (4.10a) 

123  = 723  ~ iA,  (4.10b) 

»lt3  = yi3  + i(A  - A'),  (4.10c) 

and  the  detunings  A and  A'  are  given  by 

A = ft  - co;  A ' = ft'  - to'.  (4.1 1) 


Perturbations  solutions  of  eqs.  (4.8)  are  easily  obtained.  The  six  chains  leading  to  a population 
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p33(v)  of  order  (xx'f  are 


Pi  1 * P 12  Pll  **  PlS  *+  P33> 

(4.12a) 

Pll  ^Pll  + Pll  **  Pl3  *+  P33> 

(4.12b) 

Pll  * Pll  * Pll  **031*+  P33’ 

(4.12c) 

Pll  ^Pll  * Pll  *+  Pll  4 P33> 

(4.1 2d) 

Pi  1 ^ Pll  **  Pl3  ^ Pl3  ^ P33* 

(4.12e) 

, , 

Pll  * Pll  ■*»  P31  * P31  **  P33. 

(4.120 

and  the  corresponding  Feynman  type  diagrams  are  shown 
intermediate  state  population  p22(v)  and  may  be  grouped 

in  fig.  9.  Chains  (4.1 2a— d)  involve  the 
as 

Pu  ^ ^‘2  "*  P11  ' * \ 23  "P33  sw. 

''Pii'  ''P31'  

(4.13) 

I refer  to  the  sum  of  these  four  terms  as  a stepwise  (SW)  contribution  to  the  line  shape.  The  term 
“stepwise”  is  used  since  (4.13)  involves  the  sequential  absorption  of  fields  E and  E.  The  remaining 
two  terms  (4. 1 2e,  0 

- y Pn  Pis  Pis  ^ ~ -rr\ 

PuC  - - ~ ^ P33  T Q 

Pll  P 31  P31  S 

(4.14) 

contain  ^}(vj  [or  p3,(»’)]  rather  than  p22(v).  I refer  to  the  sum  of  these  two  terms  as  a two-quanta 
(TQ)  contribution  to  the  line  shape.  While  this  label  is  somewhat  arbitrary,  it  is  meant  to  imply 
that  TQ  contributions  are  identically  zero  unless  two  quanta  (one  from  each  field)  act  on  the  system 
at  the  same  time. 

The  break-up  into  TQ  and  SW  chains  will  be  particularly  useful  when  collisions  are  included. 
For  the  no-collision  steady-state  calculation,  it  is  artificial  to  distinguish  the  TQ  and  SW  chains. 
The  no-collision  calculation  can  be  carried  out  using  probability  amplitudes  rather  than  density 
matrix  elements.  In  the  amplitude  calculation,  there  is  nothing  to  distinguish  the  TQ  and  SW 
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Fig.  9.  Diagrams  illustrating  the  absorption  of  fields  E and  E taking  an  atom  from  level  I to  level  3.  Each  diagram  is  simply  a different 
’time-ordering"  of  the  same  diagram. 
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chains.  I shall  return  to  this  point  in  discussing  the  resonance  conditions  for  the  atom-field  inter- 
actions. Note  that  each  diagram  of  fig.  9 is  simply  a different  “time-ordering”  of  the  same  diagram. 

A straightforward  solution  of  eqs.  (4.8)  according  to  the  perturbation  chains  (4.13)  and  (4.14) 
yields  a population  density  p33(v) 

P33(r)  = P?»  + P3?W,  (4.15a) 

where 

_ 4 (XX’)2NW0(v)  yl2 723 , , . 


x r . • n a 77 — 7, 77 — n + conJu8ate  (4.15c) 

[yis  + i {(A  - A)  - (k  - ek)vx} ] 

are  the  contributions  from  eqs.  (4.13)  and  (4.14),  respectively.  The  SW  term  (4.15b)  is  simply  the 
product  of  absorption  cross  sections  for  the  fields  E and  E.  The  resonance  conditions  are  those 
for  the  individual  absorptions,* 

A = kvz  ± y2/2,  (4.16) 

A'  = ek'vt  ± y23,  (4.17) 

written  for  the  case  y,  = 0.  For  the  TQ  term  (4.15c)  there  is  a characteristic  two  photon  resonance 
condition 


(ft  - kvx)  - (ft'  - ek'v,)  = (o  - co'  ± y3/2.  (4.18b) 

There  are,  in  addition,  resonances  given  by  eqs.  (4.16)  and  (4.17).  At  first  glance,  it  seems  that, 
for  fixed  A and  v2  there  should  be  two  A'  resonances,  given  by  eqs.  (4.17)  and  (4.18),  respectively. 

However,  when  state  1 has  zero  width  (y,  =0),  eqs.  (4.15b)  and  (4.15c)  may  be  combined  to 
give 

- (u’)2NW0(v) (4.19) 

Pi3'  [(72/2)J  + (A  - kvt )2]  {(y3/2)2  + [(A  - A')  - (k  - efc> J2}  ’ 

showing  that  there  is  only  one  A'  resonance  (4.18a).  This  result  may  be  understood  by  considering 
the  excitation  as  a single  process  rather  than  breaking  it  up  into  TQ  and  SW  components.  Starting 
from  the  ground  state  in  fig.  10,  the  monochromatic  field  E adds  an  energy  Wl[l  - (uz/c)]  in  the 
atom's  rest  frame.  The  field  E must  take  the  atom  from  this  energy  to  the  final  state  of  the  system 
(see  fig.  10).  For  fixed  ft  and  v,  the  process  is  resonant  when  (ft  - kvt)  - (ft'  - ek'v.)  - w - co', 
which  is  just  condition  (4.18).  The  width  of  level  2 is  unimportant  for  the  resonance  condition  on 
ft  owing  to  the  fact  that  the  fixed  frequency  ft  brings  the  atom  to  a specific  energy  within  level 


A resonance  condition  written  in  the  form  A =»  A ± B is  a shorthand  notation  for  |A  - ,4|  5 B. 
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Fig.  10.  Schematic  representation  of  absorption  of  photons  from  fields  £ and  E'  taking  an  atom  from  a zero  width  ground  state  to 
a finite  width  final  state.  The  curves  alongside  levels  2 and  3 symbolize  the  natural  widths  of  those  levels.  For  a fixed  detuning  (O  - to) 
and  longitudinal  atomic  velocity  absorption  is  optimal  when  (Q  - kvj  - (O'  - ck'v,)  - oi  - o>  ± y,/2. 


2’s  width.  There  is  no  mechanism  (i.e.,  collisions)  to  allow  for  a redistribution  of  the  energy  in  level  2. 
The  existence  of  a single  resonance  condition  on  A'  is  evidence  that  the  break-up  into  TQ  and  SW 
chains  is  artificial  in  the  no-collision  case. 

The  line  shape  is  obtained  from  (4.19)  as 


v.i)-} 


( XX')2N  f J 

~ nt,2u  J dV: 


— exP  [—(»>>!] l (4  20) 

[(?2/2)2  + (A  — kvz)2]  {(y3/2)2  + [(A  - A')  - (k  - ek')^]2}  ‘ - 1 


While  it  is  possible  to  express  eq.  (4.20)  in  terms  of  plasma  dispersion  functions,  the  physics  con- 
tained in  eq.  (4.20)  becomes  clearer  if  various  limiting  cases  are  considered.  In  all  cases,  the  Doppler 
limit 

y tj  < ku  (4.21) 

is  assumed  and  the  results  are  interpreted  in  terms  of  the  resonance  conditions  (4.16)  and  (4.18). 

1.  | A | > ku.  For  pump  detunings  larger  than  the  Doppler  width,  eq.  (4.16)  is  equally  well  satis- 
fied for  all  v2  since  the  detuning  A is  large.  Thus,  the  allowed  range  of  contributing  v2  is  dictated  by 
a distribution  centered  at  v2  = 0 with  width  u(ln  2)1/2.  If  this  distribution  of  vt  is  inserted 
into  eq.  (4.18a)  one  predicts  a line  shape  centered  at  A'  = A with  a width  equal  to  the  convolution 
of  a Gaussian  of  width  ( k - ek')u  (In  2),/2  with  a Loren  tzian  of  width  y,3  [the  sum  {(k  -Ek')vt  ±y13 } 
is  represented  in  the  line  shape  as  a convolution].  The  limiting  form  of  eq.  (4.20)  for  |A|  > ku  is 


/(A',  A)  ~ 


N_  f d exp[-(e,/u)2] 

*2  J ‘ (Vj/2)2  + [(A  - A')  - (k  - ek')v ,]2 


n1  2uA2  J 

— CD 

2 (XX')2N 
y3  |k  - ek'|  uA 


^ /A'  - A + \iy3\ 
lZ,\\k-sk'\u  } 


nich  verifies  the  prediction. 

2.  | A|  < ku.  For  near  resonant  tuning  in  the  Doppler  limit  (4.2 IX  there  is  a velocity  subset  of 
atoms  centered  at  v,  = A/k  with  width  yl2/k  that  is  resonant  with  the  field  E and  fulfills  condition 
(4. 16).  Atoms  with  these  velocities  provide  the  major  contribution  to  the  line  shape.  Inserting 
this  distribution  of  v , in  the  resonance  condition  (4.18aX  one  predicts  a line  shape  centered  at 
A'  = £(k'/k)A  with  a width  { [|k  — ek'|/k]y2/2  + yj/2}.  The  limiting  form  of  (4.20)  for  |A|  < ku. 
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ft,  -4  ku  yields  a Lorentzian  line  shape  with  precisely  these  features 


r,A.  AV  47c,/2to')2^exp[-(A//cu)2]  r2 

y2y3*«r  r2  + a2’ 

(4.23a) 

= | k-  ek'\y2  y3 

k 2 2* 

(4.23b) 

5 = A'  - e(k'/k)A. 

(4.23c) 

3.  |A|  « ku.  For  | A|  % 3 ku.  the  resonance  condition  (4.16)  is  satisfied  in  two  ways.  First,  atoms 
having  v2  within  the  central  region  of  W0(v ,)  contribute.  These  atoms  are  detuned  by  |A|  > ku 
and  give  rise  to  a resonance  as  in  the  |A|  > ku  case.  Second,  there  is  a small  number  of  atoms 
(weighting  factor  of  exp  [— (A/Icu)2])  having  v,  = A /k  that  satisfy  (4.16)  and  lead  to  a resonance 
as  in  the  | A | < ku  case.  For  | A|  % (2-A)ku  the  two  contributions  are  about  equal  and  one  sees  two 
resonances  [13],  with  the  line  shape  given  approximately  by  the  sum  of  (4.22)  and  (4.23). 

In  fig.  1 1,  the  line  shape  for  the  three  cases  is  illustrated  for  copropagating  waves.  The  line  widths 
for  copropagating  waves  are  smaller  than  the  corresponding  widths  for  counterpropagating  waves, 
owing  to  the  partial  cancellation  of  Doppler  phase  shifts  in  the  copropagating  case. 

4.2.  Collisions  - general  considerations 


A formal  expression  for  the  line  shape  incorporating  collision  effects  can  be  written  by  inspection 
if  one  uses  the  techniques  developed  in  section  3,  along  with  the  diagrams  of  fig.  9.  In  each  diagram, 
one  starts  with  p,  |(r0)  = yiNW0(r0).  The  collision  propagator  G1,(v0-»p1)  defined  through 
eq.  (3.32),  propagates  y , Af M^>( r0)  to  the  first  field  interaction.  The  appropriate  -►  r2)  [G,2 
for  diagrams  9a,  c,  e;  G21  for  diagrams  9b,  d,  f]  propagates  the  solution  between  the  first  and  second 
field  interactions.  The  process  continues  through  all  the  field  interactions  with  G33(p'  -»  »>)  propa- 
gating the  solution  following  the  last  field  interaction.  The  G^v'  -*  p)  are  obtained  as  solutions 
of  eq.  (3.32)  with  the  r\tJ{vz)  given  by 

Hu  = ft;  tiulvj  = ti2i(v,)*  = nn-  i 

nislv,)  = nu(v,)*  = m3  + nu(v,)  = 'faiM*  = nn  - '(*  - (4.24) 

and  the  ^ defined  by  eq.  (4.10). 

<•1  k*2k'  y,  * y , * 0. 1 ku  y,  * 0 


Fig.  II.  Ntvcollision  line  shapes  for  copropagating  waves  in  a three-level  system.  For  large  detunings  A J 4iu.  the  line  shape  is  a Voigt 
profile  centered  near  A'  = A.  For  a detuning  A % (2  3Mcu,  the  line  shape  consists  of  the  sum  of  a Voigt  profile  centered  near  A'  - A 
plus  a Lorentzian  centered  at  A’  = *(*'/*) A.  For  A < ku.  the  line  shape  is  a Lorentzian  centered  at  A'  = e(/t7(c)A.  The  intensity  is  in 
arbitrary  units,  but  the  same  scale  is  used  for  the  three  line  shapes. 
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The  line  shape  is  obtained  by  associating  a coupling  constant  ( ± i*,  ± i*')  with  each  field  inter- 
action and  integrating  the  product  of  Gj/s  over  initial,  intermediate,  and  final  velocities.  Grouping 
the  results  according  to  eqs.  (4.13)  and  (4.14),  one  finds 

/(A',  A)  = pf£  + pH  (4.25a) 

with 


„sw 
P 33 


= (xx')2J 


dv  dr4  dr3  dv2  dr,  di’0G33(»’4  -»  f) 

x [GzaK  - *’4)  + G^s  -»  •,4)]G22(»'’  -»  *3) 
x [Gi 2(^1  - >-2)  + G 2,(1-,  -►  •'2)]G,,(i'o  - »’i)[>’i^W/0(p0)], 


(4.25b) 


nTQ 

P33 


= (zx')2  J' 


dv  dv4  d*>3  dv2  d»,  dv0  G33(p4  - p) 

X [G23(^3  -»  •’4)Gl3(»’2  -»  »’3)G,2(»’l  -»  *2) 

+ G32(p3  - r4)G3,(p2  -»  *3)G21(r,  -»  ^2)]G , ,(^0  -*•  •’t)  [y , IV0(v0)]- 


(4.25c) 


There  are  six  separate  terms,  one  for  each  diagram  of  fig.  9. 

By  using  eqs.  (3.36)  and  (3.37)  to  perform  the  integrals  over  Gu  and  G33,  respectively,  along 
with  the  fact  that 


Gj/p'  -»  p)  = Gjj(p'  -»  p)*, 
one  reduces  eq.  (4.25)  to  the  simplified  form 

/(A',A)  = pH  + pI? 


(4.26) 

(4.27a) 


„SW 

P 33 


Pi? 


= — [dvdvjdv,  dp0[G23(p2  -♦*’)  + G23(p2  “♦  »’)*] 

Vs  J 

x G22(»'i  -»  ^2)  [G 1 2(^0  -*  rj)  + GjjK  -♦  *’i)*]M/o(»'o)  (4.27b) 

{Yy')*N  r 

= dvdp2  dp,  dp0G23(P2  -4  *')G13(p,  -►  P2)G,2(p0  -*•  p,)M^0(p0)  + conjugate. 

73  J 

(4.27c) 

Equations  (4.27)  provide  a formal  expression  for  the  line  shape.  Recall  that  /(A',  A)  is  a measure 
of  the  probe  field  absorption  in  the  presence  of  a pump  field  when  collisions  are  occurring  in  the 
system.  Eq.  (4.27)  was  derived  assuming  that  atoms  start  in  a thermal  velocity  distribution  in  level  1 
and  that  the  pump  and  probe  fields  are  weak. 

The  line  shape  (4.27)  depends  on  the  propagators  G 1 2,  G22<  G23  and  Gi3.  The  diagonal  propa- 
gator G22  enters  the  SW  chain  (4.27b)  and  is  a measure  of  velocity-changing  collisions  that  active 
atoms  undergo  while  in  level  2.  The  stepwise  absorption  process  involving  G22  is  easy  to  visualize 
absorption  of  field  E taking  the  atom  from  state  1 to  2,  velocity-changing  collisions  in  level  2, 
absorption  of  field  E'  taking  the  atom  from  level  2 to  3. 

The  ofT-diagonal  propagators  G12  and  G23  appear  in  both  the  SW  and  TQ  chains,  while  G13 
appears  only  in  the  TQ  chain.  These  propagators  contain  the  effects  of  collisions  on  off-diagonal 
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density  matrix  elements.  The  line  shape  will  take  on  different  forms,  depending  upon  whether 
collisions  are  phase-interrupting,  velocity-changing,  or  both  phase-interrupting  and  velocity- 
changing (in  some  inseparable  combination)  in  their  effect  on  the  various  off-diagonal  density 
matrix  elements.  Rather  than  examine  all  the  different  combinations  of  either  phase-interrupting 
or  velocity-changing  effects  on  each  off-diagonal  propagator,  I consider  only  two  specific  collision 
models  to  illustrate  the  basic  physical  principles.  In  both  models,  collisions  are  assumed  to  be 
phase-interrupting  in  their  effect  on  pl2  and  p23,  as  would  be  the  case  if  state  2 scattering  is  sub- 
stantially different  from  state  1 and  state  3 scattering.  Consequently,  G12(v'  -*•  v)  and  C23(v'  ->•  v) 
are  given  by  equations  similar  to  eq.  (3.40).  In  collision  model  /,  Gl2(v'  -►  v)  is  also  given  by  an 
equation  similar  to  (3.40)  - collisions  are  phase-interrupting  in  their  effect  on  p13.  In  collision 
model  2,  collisions  are  velocity-changing  in  their  effect  on  p13  (i.e.,  equal  collision  interaction  for 
levels  1 and  3,  as  might  be  the  case  if  they  belonged  to  the  same  electronic  configuration).  In 
collision  model  2,  the  propagator  G13(v'  -*  v)  is  given  as  a solution  of  eq.  (3.32)  with  a real  kernel 
and  rate. 

One  additional  approximation  is  adopted  in  both  models.  As  a good  first  approximation, 
the  speed  dependence  of  the  T./s  is  neglected.  (Although  neglected  in  this  work,  the  speed  depen- 
dence of  the  T./s  can  play  an  interesting  role  in  saturation  spectroscopy.  For  near  resonant  pump 
fields  < ku,  atoms  with  a velocity  component  v.  = A /k  are  selected  by  the  field.  Since  |r.|  oc  A, 
the  energy  or  temperature  dependence  of  the  Tj/s  and  the  corresponding  cross  sections  may  be 
determined  by  varying  A.  Thus,  in  a gas  at  a given  temperature,  one  can  measure  cross  sections  as 
a function  of  temperature  by  choosing  hot  or  cold  atoms  using  a suitable  pump  detuning.  The 
analysis  is  complicated  by  the  fact  that  the  transverse  velocity  components  are  not  selected  by 
the  fields,  so  that  the  line  shapes  must  be  averaged  over  these  components  [14].) 


4.3.  Line  shape  - collision  model  1 

If  collisions  are  phase-interrupting  in  their  effect  on  off-diagonal  density  matrix  elements,  the 
off-diagonal  propagators  are  given  by 


%v  - v') 


G’V  9)  ( V.j  + Hf)  + I(Ay  + S$) 

A12  = A - kvx ; A23  = -A'  + ek'vt;  A13  = A12  + A23 


(4.28a) 

(4.28b) 


with  H*  and  Sff  defined  by  eq.  (3.41)  (the  speed  dependence  of  these  parameters  is  neglected). 
Inserting  eq.  (4.28)  into  (4.27),  one  obtains  the  line  shape 


/(A\A)  = pf?  + pJ? 


(4.29a) 


,s  - + <3 ">°l- 


(4.29b) 


TO  _ (XX  ) 

P33  - 


y3  J f23  - i(i'  - Bk'vt)f  13  + i[A  - S'  - (k  - £/c')dz]  f12  + i(^l  - kvt)W°^ 
+ conjugate,  (4.29c) 
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where 


fy  = Yu  + Hi*. 
5 = A + sf2, 
A'  = A'  - s?h3, 


A - A'  = A - A'  + S?h3. 


(4.30) 

(4.31) 

(4.32) 

(4.33) 


Rather  than  analyze  the  line  shape  for  the  most  general  case,  1 shall  consider  (4.29)  in  the  limits 
of  large  pumping  detuning  |A|  > Icu  and  near  resonant  tuning  | A|  < ku. 

In  discussing  the  no-collision  line  shape  for  y,  = 0, 1 stressed  that  there  was  only  one  resonance 
condition  for  A'  [(4.18)].  The  conservation  of  energy  argument  regarding  internal  states  used  to 
arrive  at  that  conclusion  can  no  longer  be  maintained  since  internal  energy  can  be  transformed 
into  translational  energy  during  a collision.  Thus,  there  are  two  resonance  conditions  for  A' 


(4  34) 

(4.35) 

(4.36) 


A'  = A — {k  — ek')vz  + S?3  ± F13, 

A'  = ek’v z + Sf3  ± r23, 
along  with  the  resonance  condition  on  A 
A = kvz  - S?2  ± f 12. 

Eqs.  (4.34)  and  (4.35)  may  be  thought  of  as  two-photon  and  single  photon  resonance  conditions, 
respectively,  although  eq.  (4.35)  also  enters  the  TQ  chain. 

1.  | A j P ku. 

SIT  chain.  The  SW  chain  consists  of  the  absorption  of  field  £,  velocity-changing  collisions 
in  level  2,  and  the  absorption  of  field  E.  For  large  detunings,  the  detunings  A — kvz  & A of  different 
velocity  subsets  are  approximately  equal  so  that  the  entire  Maxwellian  velocity  distribution  of 
active  atoms  is  equally  excited  into  level  2 (albeit  with  small  probability)  by  field  E.  As  previously 
noted,  collisions  do  not  alter  an  equilibrium  velocity  distribution;  consequently,  velocity-changing 
collisions  in  level  2 have  no  effect  on  the  line  shape.  The  absorption  of  £'  takes  place  from  a Max- 
wellian atomic  velocity  distribution.  All  these  features  are  contained  in  the  SW  contribution 
(4.29b)  which,  for  | A|  §>  ku  and  f 12  < |A|  becomes 

_sw  _ 4fxx')2;vr12  r r23 

J d (T^FTfS^eife7^ 


Pi  3 = 


YiYi& 


W0(v), 


(4.37) 


where  eq.  (3.36)  has  been  used. 

TQ  chain.  I can  not  offer  a simple  physical  interpretation  of  the  TQ  chain.  Since  it  is  not  positive 
definite,  the  TQ  chain  does  not  correspond  to  an  excitation  probability.  Instead,  it  represents  an 
interference  term  that  contributes  to  the  total  level  3 population.  The  TQ  contribution  (4.29c) 
is  a product  containing  all  three  resonance  conditions  (4.34H4.36).  For  large  detunings,  condition 
(4.36)  can  not  be  satisfied,  condition  (4.34)  can  be  satisfied  near  A'  = A and  condition  (4.35)  can 
be  satisfied  near  A'  = 0.  Thus  (4.29c)  exhibits  two  resonances  in  A',  the  one  near  A'  = 0 combining 
with  the  SW  contribution.  For  large  detunings,  eq.  (4.29c)  reduces  to 


pIq3 


= ( w)2n  r 

y3(iA)  J 


1 


1 


f23  - i(A'  - ek'vjfn  + i[A  - A'  - (k  - £*>,] 


W0(v) 


+ conjugate. 

(4.38) 
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Line  shape.  Combining  eqs.  (4.37)  and  (4.38),  one  obtains  the  line  shape 

„a,  „ (XX’)2N  LJ  4f12f2,(y2)-1 

. 4) — j d.  + (j.  _ ^ 

+ k - W-'I  ^f„  + .[A  - i - »-  + «*-»}]»«* 

(4.39) 

which  exhibits  two  resonances.  At  A'  a A,  only  the  TQ  chain  is  important  and  (4.39)  reduces  to 


/(A',  A)  ~ 


(XX')2N  „ f A' 

-£fc>A2  ‘V 


A'  - A - S?h3  + if  13> 


k — ek'\u 


for  A'  % A, 


while  in  the  vicinity  of  A'  = 0 both  the  TQ  and  SW  chains  combine  to  give 


/(A',A)~ 


2(xx')2Jvr  2(bz  + rfi) 


— - 1 


A'  - Sft  + if23> 
k'u  j 


A'  as  0. 


y,MlL  y,  J \ k'u  J 

Near  the  two-photon  resonance,  A'  a A,  the  line  profile  (4.40)  is  broadened  and  shifted  from 
the  corresponding  resonance  (4.22)  of  the  no-collision  case.  Line  shape  studies  in  the  region  of 
A'  « A allow  one  to  determine  r?3  and  S?3. 

As  predicted,  there  is  a new  resonance  at  A'  a 0.  In  the  no-collision  limit,  PJ2  = 0,  and  the  TQ 
and  SW  contributions  cancel  one  another.  With  increasing  perturber  pressure  (or,  equivalently, 
increasing  rft),  the  SW  component  grows  relative  to  the  TQ  one  leading  to  absorption  near 
A'  = 0.  For  a fixed  detuning  A,  the  increase  in  the  SW  contribution  reflects  an  additional  absorption 
of  the  field  £,  owing  to  the  collisional  increase  in  the  damping  associated  with  that  absorption. 
The  amplitude  of  the  resonance  provides  a measure  of  r?2  while  the  width  and  shift  can  be  used 
to  obtain  values  for  rft  and  respectively.  A typical  piofile  is  shown  in  fig.  12. 

It  is  interesting  to  note  that  the  integrated  profile  J /(A',  A)  d A'  increases  with  pressure,  owing 
to  the  presence  of  the  new  resonance  at  A'  « 0.  Thus,  if  the  probe  frequency  is  swept  through  both 
resonance  regions,  the  total  number  of  atoms  excited  to  level  3 will  increase  with  pressure.  More- 


Fig.  '2.  Line  shape  (4.40,  4.41)  for  A •>  4ku  and  copropagating  waves  when  collisions  are  phase-interrupting  in  their  effect  on  pts. 
The  profiles  are  drawn  assuming  y,  - 0,  |yf$|  « (k  - k')u  and  |yf5|  < k’u.  The  line  shape  consists  of  Voigt  profiles  oentered  near 
A’  - A and  A’  - 0.  Compare  this  profile  with  the  A - 4ku  no-collision  profile  in  fig.  1 1.  Collisions  give  rise  to  a new  resonance  centered 
near  A’  — 0. 
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over,  excitation  occurs  in  a new  region  (A'  & 0)  of  the  probe  frequency  spectrum.  Collisions  give 
rise  to  an  increased  range  of  probe  frequencies  for  which  level  3 excitation  is  possible. 

2.  | A | < ku. 

In  all  cases,  it  is  assumed  that 

Sj’h,  f < ku.  (4.42) 

SW  chain.  For  |A|  < ku,  there  is  always  a velocity  subset  of  atoms  resonant  with  field  E and 
this  field  excites  atoms  having 

vt  = (A  + Sf2)/Jc  ± (f  I2/fc).  (4.43) 

The  population  p22<uz)  excited  by  E is  shown  schematically  in  fig.  13a. 

Velocity-changing  collisions  result  in  a relaxation  of  this  distribution  back  towards  an  equili- 
brium one.  Regardless  of  the  detailed  nature  of  the  collision  kernel  W22(r'  -*■  v)  giving  rise  to  this 
relaxation,  the  degree  of  relaxation  is  determined  by  the  number  of  collisions 

<n>  = r22/y2,  (4.44) 

occurring  within  the  lifetime  of  level  2 and  the  root  mean  square  change  in  velocity  per  collision 
(Au).  The  distribution  p22(r.)  following  <m>  collisions  is  shown  schematically  in  fig.  13b  as  a 


Fig.  13.  Stepwise  contribution  to  the  line  shape  for  |A|  < ku  (or  total  line  shape  Tor  countcrpropagating  waves)  when  collisions  are 
phase-interrupting  in  their  effect  on  p,s.  (a)  The  Lorentzian  velocity  distribution  of  atoms  excited  by  field  E is  centered  at  e, « A/k 
with  width  T,,/*.  (b)  Redistribution  of  velocities  in  level  2 resulting  from  collisions.  These  profiles  were  calculated  using  a collision 
kernel  W(r'  -»  »)  » (OAiuJn)'*  exp  { -[(r  - 0.8»)/0.6u]J}  having  a (Au)  = 0.6u.  With  increasing  <n>  - r22/y2  (or.  equivalently, 
with  increasing  pressure),  the  velocity  distribution  tends  towards  thermal  equilibrium,  (c)  Probe  absorption  line  shapes  as  a function 
of  <a>.  The  line  shapes  represent  the  absorption  of  field  E'  from  the  level  2 population  distributions  shown  in  (b).  The  profile  is  a 
Lorentzian  centered  at  A'  = tik'lk)t A U is  taken  equal  to  - 1 in  the  figure)  at  low  pressures  and  tends  towards  a Voigt  profile  centered 
at  A'  * 0 at  high  pressures.  Vertical  scales  are  in  arbitray  units,  but  the  same  scale  is  used  for  each  diagram  in  (b)  and  in  (c). 
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m 


function  of  <n>  assuming  k{Au)  > f ,2.  (As  long  as  k(Au)  > f ,2>  each  collision  is  strong  enough 
to  he  distinguished  since  it  can  remove  an  atom  from  the  velocity  distribution  excited  by  field  E. 
In  this  limit,  it  is  the  true  rather  than  the  effective  collision  rate  that  is  pertinent  to  the  discussion.) 
For  low  pressures  (<n>  :$  IX  the  distribution  is  a narrow  Lorentzian,  but  with  increasing  pressure 
(<n>  P 1)  the  distribution  tends  towards  a Maxwellian.  (In  many  cases,  at  pressures  required  to 
thermalize  level  2,  the  Doppler  limit  condition  (4.42)  no  longer  holds  and  the  distributions  will  be 
somewhat  more  complicated  than  those  shown  in  fig.  13b.) 

The  probe  field  absorption  resonance  condition  is  given  by  (4.35)  with  v.  obtained  from  the 
velocity  distributions  of  fig.  13b.  The  corresponding  line  shapes  are  depicted  in  fig.  1 3c.  For  low 
pressures,  e.  is  determined  by  eq.  (4.43X  which,  when  substituted  in  eq.  (4.35),  shows  that  A is 
resonant  for 


A'  = **'/*)( A + Sft)  + S£h3  ± (f23  + (k'/k)tx2). 


(4.45) 


For  higher  pressures,  the  probe  absorption  profile  is  the  convolution  of  p22(v)  with  a Doppler 
shifted  Lorentzian  of  width  f23;  it  becomes  broader  and  its  center  shifts  back  towards  S'  = 0, 
ultimately  reaching  the  limit  of  a Voigt  profile.  There  is  no  difficulty  in  understanding  the  SW 
contribution  on  the  basis  of  simple  physical  arguments. 

The  form  of  the  SW  contribution  depends  on  the  specific  nature  of  the  kernel  W22(v'  ->  v). 
For  an  arbitrary  kernel,  however,  eq.  (3.32)  may  be  formally  solved  to  give  the  propagator 


G22(v'  - p)  = <5(v  - v')/(y2  + T22)  + <722(p'  - v ), 
where 

<722(„'  - „)  = -j—  fdr" W22(v"  - p)G22(p'  -»  v"). 
?2  + 1 22  J 


(4.46) 


(4.47) 


The  first  term  in  eq.  (4.46)  gives  the  contribution  from  atoms  which  have  not  undergone  velocity- 
changing collisions  in  level  2 before  being  excited  to  level  3.  The  second  term  contains  all  the  effects 
of  velocity-changing  collisions.  By  substituting  eq.  (4.46)  into  eq.  (4.29b)  and  using  eq.  (4.42)  one 
obtains  the  SW  contribution 


sw  = 4nvl(XX')zN  exp  [ -(Mm)2] 
P3i  (72  + r22)y3kura 


..  f <ra)2  «ya  + r22)r.f  f23 

|(ra)2  + (Aa)J  + nV  rdV°  JTT^W-ek'vy 


where 


r.  = r23  + (fc7fc)f,2;  A,  = S’  - e(k’/k)&. 


(4.48) 

(4.49) 


which  contains  a Lorentzian  contribution  from  atoms  that  have  not  undergone  velocity-changing 
collisions  in  level  2 plus  a term  reflecting  the  velocity  relaxation  in  level  2.  One  may  calculate 
(j22{*'  -»  v)  for  specific  collision  kernels,  but  such  calculations  are  not  presented  in  this  work. 
The  general  nature  of  the  solution  is  shown  in  fig.  13c. 


> 1^- — - — :■ - - 
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TQ  chain.  In  the  Doppler  limit  (4.42),  H^v)  appearing  in  eq.  (4.29c)  is  evaluated  at  v.  = Kjk 
and  the  resulting  TQ  contribution  to  the  line  shape  is 


tq  _ (XX')2N  cxp[-(K/ku)2~\ 


P 33  = 


7I1/2UV3 


00 

I 


dr. 


f23  ~ - ek'vz) 


f,3  + i [S^TA7  - (k  - cfc>r]  f12  + i(5  - kvz)  + conju8ate- 


(4.50) 


For  counterpropagating  waves  (e  = - 1)  all  the  poles  are  in  the  same  half  plane  and  the  integral 
vanishes.  For  copropagating  waves  (e  = +1)  and  k'  > k the  pole  at 

= (5  - if12)/*  (4.51) 

contributes  in  the  contour  integration  and  corresponds  to  satisfying  condition  (4.36).  When  this 
value  of  vz  is  substituted  in  the  remaining  conditions  (4.34),  (4.35),  one  finds  A'  resonances  at 


k'  k'  — k h . k'-k~  1 

A - k A + * s?h2  + sft  ± |^13  + nr~ri2J 

(4.52) 

+ + S5‘,±[r„  + ^r11]. 

(4.53) 

Similarly,  for  k > k\  the  pole  at 

= (5'  + ir23)/*'. 

(4.54) 

corresponding  to  condition  (4.35),  contributes  in  the  integration  and  leads  to  resonance  conditions 

A' " + nr-sfi  + ts*  ± [f  r.> + nr-r»} 

(4.55) 

A'_74+fS!5  + ss5±(r“  + fr'i) 

(4.56) 

following  from  (4.34),  (4.36).  For  each  case  (k  > k'  or  k < k')  there  are  two  resonances  with  some- 
what different  central  positions  and  widths. 

On  integrating  (4.50)  over  v2,  one  obtains  the  TQ  contribution 


tq  = 47tl/2(xx')2Nexp[— (A/M2]  (r,rb  - A,Ab)  . 

P33  y3ku  or.)2  + (A.)2]  orb)2  + (Ab)2]  U[e)' 

where 

r.  = r23  + j r12;  a.  = a-  - £(a  + s?$)  - s§h3, 

Fb  = £ [t  r>  j + k~ir  r“]0(fc " k']  + [fij  + nr fi2] ~ ki 


(4.57) 


(4.58a) 


(4.58b) 
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[k'  k'  - k 1 

A'  - j A - Sft  - ^y-^S?h2  0(/C>  " k)' 

•Hi  r;:- 


(4.58c) 


This  line  shape  possesses  the  resonance  positions  and  widths  predicted  in  eqs.  (4.52-4.56). 

To  get  some  insight  into  this  equation,  one  can  multiply  by  (y2  + T22)/(y2  + T22)  and  perform 
algebraic  manipulations  to  obtain 


Pl?  = 


4k1  2(XX')2N  exp  [-(5/fcu)2] 


(y2  + r22v/3fc«r. 


f__  <!Z_ 

1 (rj2  + ( Aa)2 


4-  * » + 

(Tb)2  + (Ab)2 


ra[(^r)(rarb  - AaAb)  + qssxta  - 
[(Ta)2  + (Aa)2]  [(rb)2  + (Ab) 


+ rbAJ)]{ 

>J]  f 


where 


(<*n  = -(T?h2  + rshj  - rft  - r22)^0(fc  - *')  + ©(*'  - k)J, 

(<«?)  = -(S?h2  + Sft  - S?h3)|j©(*  - k ’)  + ©<*'  - *)J. 


(4.61a) 


(4.61b) 


In  this  form,  there  are  three  contributions  to  the  TQ  contribution.  The  first  term  just  cancels 
that  part  of  the  SW  arising  from  atoms  which  have  not  undergone  velocity-changing  collisions  in 
level  2 [first  term  of  eq.  (4.48)].  The  second  term  is  a Lorentzian  with  a narrower  width  than  the 
first  term,  owing  to  a partial  cancellation  of  the  Doppler  phases  in  the  resonance  condition  (4.34) 
leading  to  this  term.  The  third  term  results  solely  from  collisions  since  (<5D  and  (<5S)  vanish  in  the 
no-collision  limit.  By  using  eqs.  (4.61,  3.41,  2.34)  and  the  optical  theorem,  one  can  estimate  that 
the  magnitude  of  this  term  is  relatively  small  and  that  it  vanishes  exactly  when  level  2 experiences 
the  strongest  collision  interaction  of  the  three  levels.  If  the  term  is  nonvanishingly  small,  it  adds 
an  asymmetry  to  the  TQ  contribution. 

Line  shape.  The  line  shape  obtained  by  combining  eqs.  (4.48)  and  (4.60)  is 


rrv  m 47t,/2u/)2Nexp[-(S/M2]  [ (ra)2  r„rb 

,(4'4’- to+f^-Srr--- ^ e,  + <r»r+WOW 

kifi  + r22)Ta  j*  . A (r  . 

+ Jdfdro(r2,)*+(l-«lk'r7  12(  0 


(Tl2)2 -ma  - kt>j2  L V «2  7J 

r,[(^r)(r,r„  - A,Ab)  + (<5S)(r,Ab  + Aarb|ai<.t| 
[(r.)2  + (A.)2]  [(Tb)2  + (Ab)J]  ^ 'j- 
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For  counterpropagating  waves  (e  = - 1),  the  TQ  chain  does  not  contribute  and  the  line  shape  is 
represented  schematically  in  fig.  1 3c  for  various  pressures.  There  is  a Lorentzian  contribution  from 
atoms  which  have  not  undergone  velocity-changing  collisions  in  level  2 plus  a “shoulder”  arising 
from  atoms  which  have  been  partially  thermalized  by  velocity-changing  collisions.  Studies  of  the  line 
shape  as  a function  of  pressure  can  provide  values  for  S?2  and  the  characteristic  kernel  parameters 
(Am)  and  <w>,  assuming  that  r53,  r^,  Sft  and  have  been  measured  from  large  detuning 

experiments.  There  has  been  recent  experimental  verification  of  these  line  profiles  [15]. 

Assuming  that  the  third  term  in  (4.62)  is  small,  the  line  shape  for  copropagating  waves  (e  = 1) 
has  the  same  general  features  as  that  for  counterpropagating  waves  except  that  the  copropagating 
Lorentzian  is  narrower  than  the  counterpropagating  one,  owing  to  a partial  cancellation  of 
Doppler  phases.  The  interpretation  of  the  line  shape  is  the  same  as  that  for  the  counterpropagating 
case  and  fig.  13c,  with  the  Lorentzian  components  having  somewhat  narrower  widths,  again  serves 
to  illustrate  the  line  shape  as  a function  of  pressure.  If  the  third  term  in  eq.  (4.62)  is  non-negligible, 
there  will  be  additional  asymmetry  added  to  the  line  shape. 

(As  an  aside,  I can  mention  that  the  line  shape  takes  on  a qualitatively  different  form  than  that 
shown  in  fig.  13c  if  velocity-changing  collisions  are  “weak”.  Velocity-changing  collisions  are 
“weak”  if  the  rms  velocity  change  per  collision  (Am)  and  the  collision  rate  T22  are  not  sufficient 
to  remove  atoms  from  the  velocity-profile  (of  width  f J2/k)  excited  by  field  £.  In  this  limit,  which 
can  occur  over  a certain  pressure  range,  the  SW  contribution  remains  Lorentzian  with  a width 
equal  to  r„  (4.49)  (as  if  velocity-changing  collisions  were  absent)  plus  a small  pressure  independent 
component  owing  to  velocity-changing  collisions.  For  counterpropagating  waves,  the  line  shape, 
as  given  by  the  SW  chain,  is  Lorentzian:  the  line  width  extrapolated  to  zero  pressure  will  yield 
an  anomalously  large  value  for  the  atomic  state  lifetimes,  owing  to  the  presence  of  the  pressure 
independent  velocity-changing  collision  contribution  to  the  width.  For  copropagating  waves, 
the  line  shape  consists  of  a narrow  Lorentzian  given  by  the  second  term  in  the  TQ  contribution 
(4.60)  plus  line  wings  representing  the  effects  of  velocity-changing  collisions  in  the  SW  chain. 

For  the  counterpropagating  case  (or,  equivalently,  the  SW  chain),  it  is  also  possible  to  measure 
a line  width  that  is  a non-linear  function  of  perturber  density.  To  see  how  this  effect  arises,  imagine 
that  at  low  pressures,  (Am)  > f 12/k  « y,2/k.  Velocity-changing  collisions  remove  atoms  from  the 
homogeneous  width,  but  the  total  line  shape  remains  nearly  Lorentzian  if  (A u)/u  <?  1.  The  width 
of  this  low  pressure  pseudolorentzian  is  the  sum  of  two  pressure  dependent  components  - one 
from  phase-interrupting  collisions  and  one  from  velocity-changing  collisions.  However,  as  the 
pressure  is  increased.  T , 2 grows  and  collisions  may  become  “weak”  [k( Am)  < f , 2],  In  that  pressure 
region,  the  only  pressure  dependent  contribution  to  the  line  width  comes  from  phase-inter- 
rupting collisions.  Consequently  the  slope  of  a line  width  versus  perturber  density  graph  will 
be  greater  in  the  low  pressure  region  than  at  higher  pressures;  i.e.,  the  line  width  is  a nonlinear 
function  of  perturber  density  [16,  2].  Such  nonlinearities  have  been  observed  experimentally 

[17]) 


4.4.  Line  shape  - collision  model  2 

In  collision  model  2,  the  propagators  G12(i>'  -*  p)  and  G23(r'  -►  r)  are  still  given  by  eq.  (4.28) 
and  the  SW  contribution  is  unchanged  from  (4.29b).  For  equal  collision  interactions  in  states  1 and 
3,  collisions  are  velocity-changing  in  their  effect  on  p|3(i>)  and  the  propagator  G\%  is  obtained  as 
a solution  of  [see  eq.  (3.32)] 
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{7,3  + nc3  + i[A  - A .'-(*-  ek')V'-]}G\%(v'  ■ 


with  a real  collision  kernel  W,7|(k'  -*  v)  and  rate  r^3.  The  line  shape  follows  from  (4.27)  as 

/(A,,A)-pS?  + p]?,  (4.64a) 


JW  MxxrN  f r23 


4 (XX')2N  j*d 

V3  J 

_ (XX')2  N f 

73  J' 


2 + (£'-  ek'v2)2  G22(V°  V*\ri2)2+(£-kvl0)2  W°(Vo)' 


drdr0 


‘23  - i(5'  - ek'vz) 


* - r)  + _ kv  ) Ho(»’o)  + conjugate. 


(4.64b) 


(4.64c) 


Again  the  discussion  is  limited  to  the  large  detuning  and  near  resonant  tuning  limits.  Only  the  TQ 
contribution  and  total  line  shape  need  be  discussed  since  the  SW  contribution  is  the  same  as  in 
collision  model  1. 

1.  J A|  P ku.  As  for  collision  model  I,  there  are  resonances  near  A'  = 0 and  A'  = A.  For  | A|  P ku 
and  r,j  4 |A|  the  propagator  G13(*>'  -*  •>)  in  the  vicinity  of  the  A'  = 0 resonance  is  independent 
of  collision  model,  owing  to  the  large  value  of  | A|;  for  either  phase-interrupting  (4.28)  or  velocity- 
changing (4.63)  collisions,  one  obtains  G13(r'  -►  v)  a £(i>  — »>')/iA.  Thus,  the  line  shape  near 
A'  = 0 is  the  same  as  in  collision  model  1 and  is  given  again  by  eq.  (4.41). 

Near  A'  = A,  the  TQ  contribution  dominates  and,  for  |A|  ku  and  fw  4 |A|,  (4.64)  reduces  to 


i.  ( XX)2N  f 

‘•4I~ 


drd»>0G"3(r0  -*•  r)lT0(»>0)  + conjugate. 


A'  % A, 


which  is  similar  in  form  to  (3.38)  with  G"3  replacing  Gv,2.  A discussion  of  the  line  shape  parallels 
that  of  section  3.4.  Collisions  which  are  velocity-changing  in  their  effect  on  p,3  give  rise  to  a narrow- 
ing of  the  two-photon  resonance.  The  narrowing  becomes  significant  when  (r,3)'  > |fc  — efc'|u, 
where  ( r\c3 )'  is  the  effective  collision  rate  discussed  in  section  3.4.  For  copropagating  waves  (e  = 1), 
the  onset  of  narrowing  occurs  at  lower  pressures  than  in  linear  spectroscopy  while  for  counter- 
propagating  waves  (e  = —1),  higher  pressures  are  required  to  see  the  narrowing. 

At  low  pressures,  the  line  shape  (4.65)  is  the  convolution  of  a Gaussian  of  width  |fc  -ek'|u  (In  2)1/2 
with  a Lorentzian  of  width  y13.  The  line  width  decreases  monotonically  with  increasing  pressure. 
For  pressures  where  (r"3)'  > \k  — efc'|2u2/y13,  the  line  becomes  a Lorentzian  with  width  y13. 
(Note  that  the  assumption  f(>  4 | A|  probably  fails  at  pressures  where  significant  narrowing  occurs; 
the  line  shape  is  modified  slightly  by  this  effect.)  At  all  pressures,  the  line  is  centered  at  A'  = A. 

The  large  detuning  limit  in  the  region  of  A'  % A serves  to  test  whether  collisions  are  phase- 
interrupting or  velocity-changing  in  their  effect  on  pl3.  If  the  two  photon  resonance  narrows  with 
increasing  pressure,  collisions  are  velocity-changing  in  nature  and  one  may  determine  the  effective 
collision  rate  (r”3)'.  If  the  line  broadens  and  shifts  with  increasing  pressure  (after  any  inelastic 
collision  effects  have  been  properly  subtracted  out),  collisions  are  phase-interrupting  in  nature  and 
one  may  determine  Tft  and  S?3. 
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2.  | A | < ku.  For  near  resonant  tuning,  the  line  shape  does  depend  on  the  specific  form  of  the 
collision  kernel  H7$(i»'  -►  v).  However,  it  is  still  possible  to  obtain  a semiquantitative  picture  of 
the  line  profile  without  going  into  specific  details  of  the  kernel.  The  discussion  is  limited  to  “low” 
and  “high”  pressure  regions. 

a.  Low  pressure  ((T^)'  <?  \k  — ek'|u;  f,,  ku).  As  in  eqs.  (4.46,  4.47)  it  is  useful  to  write  the 
propagator  G\c3  as 


G,a(r  V)  (7,3  + rr3)  + i[A  - A'  — (Ic  - ek')v J + ^ V)  (466) 


where 


(V”  - v)Cr3(v'  - r").  (4.67) 

The  lead  term  in  eq.  (4.66)  has  the  same  form  as  the  phase-interrupting  collision  propagator  (4.28) 
with  the  substitutions  T?3  Tv13,  S?3  <-+  0.  Thus,  taking  into  account  only  the  contribution  from 

the  first  term  of  eq.  (4.66),  one  obtains  the  same  line  shape  as  in  collision  model  1 [eq.  (4.62)] 
with  r?3  replaced  by  r^3  and  S?3  replaced  by  zero.  There  is  however,  an  additional  contribution 
to  the  line  profile  arising  from  the  second  term  of  (4.66).  The  total  line  shape  is  given  by 

//a'  A>  = 47t1/2to,)2iVexp[-(3/M2]  [ (ra)2  rar; 

(72  + r22)y3fcMra  |(ra)2  + (Aa)2  °(  £)  + (r;)2  + (a;)2  0(£) 

, f *(y2  + r22)ra  j*  1 

L 4rt2“2  J r°  f23  - i(A'  - ek'/v,) 

.exp  [-(i>2  + u2  )/u2]  . 1 

X Sl3iV°  V)~T72  + i(A  ~kvj  + C°njUga,eJ 

+ i72  +zj2)  a {drdro(rJi)2'+Ii'3-ek'pl)2 

. r.[(dr  )(rarb  - A.Ay  + (<5S)(r.Ab  + ) 


where 


[(T.)2  + (Aa)2]  [<r;)2  + (a;)2] 


r,>  ~ (7' 3 + nc3>  H £ — ^23^j ®(k  - k ) + £y13  + T*  + * k - k), 

(4.69a) 

= £ [A'  ~ I A ~ ’Lir  S”]0**  -^')  + [a'-|a-  SftJ©(fc'  - k ),  (4.69b) 
(<5D  = -(rft  + rsh,  - nc3  - ©(*  - k')  + ©(*'  - *)J,  (4.70a) 


(4.70a) 
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(<$S)  = -(Sft  + Sft)[y0(*  - *')  + &ik'  - ft)J 


(4.70b) 


and  ra  and  4,  are  defined  by  eq.  (4.58a).  The  third  term  in  eq.  (4.68)  represents  the  contribution 
from  collisions  which  are  velocity-changing  in  their  effect  on  p,3  while  the  fourth  term  is  the 
SW  velocity  redistribution  resulting  from  collisions  in  level  2.  The  fifth  term  is  assumed  to  be 
small. 

The  third  term  is  the  critical  one  for  the  present  discussion.  This  term  will  take  on  a qualitatively 
different  form  depending  upon  whether  collisions  are  “weak”  or  “strong”  (in  the  sense  described 
in  section  4.3). 

Copropagating  waves  (e  = 1).  Any  strong  collision  upsets  the  Doppler  phase  cancellation  that 
characterizes  the  £ = 1 TQ  chain.  Thus,  the  third  term  in  (4.68)  is  negligible  for  strong  velocity- 
changing collisions  in  the  low  pressure  region.  Strong  collisions  broaden  the  narrow  resonance  of 
the  TQ  contribution  [second  term  of  eq.  (4.68)];  these  collisions  have,  in  effect,  decreased  the 
coherence  time  associated  with  pti.  The  narrow  resonance  is  a Lorentzian  centered  near  A' =(k'/k)A 
with  a width 


_ fyj 
b [2 


(i r?i  + nr  nh3)0(k  “ k)  + (rr3  + nr^r?h2)0(*'-k) 


(4.71) 


For  weak  collisions,  the  Doppler  phase  cancellation  is  not  significantly  disturbed  to  render  the 
third  term  in  eq.  (4.68)  negligible.  This  term  combines  with  the  second  to  produce  a non-Lorentzian 
line  shape  that  broadens  with  increasing  pressure.  The  specific  form  of  the  line  shape  will  depend 
on  the  kernel  W7|(i»'  -»  »»). 

In  addition  to  the  narrow  resonance  of  the  TQ  chain  for  either  weak  or  strong  collisions,  there 
is  a background  [fourth  term  in  eq.  (4.68)]  resulting  from  velocity-changing  collisions  in  level  2. 

Counter  propagating  waves  (e  = —1).  In  the  absence  of  collisions,  the  TQ  chain  contribution 
vanishes.  Mathematically,  this  result  was  obtained  owing  to  the  fact  that  all  the  poles  in  the 
velocity  integration  were  in  the  same  half-plane  (physically,  one  might  say  that  a cancellation  of 
Doppler  phase  factors  does  not  occur  for  counterpropagating  waves).  If  weak  velocity-changing 
collisions  are  present,  the  TQ  chain  still  vanishes  to  first  approximation  since  the  velocity  change 
is  not  sufficient  to  change  the  general  nature  of  the  solution.  Even  for  strong  collisions,  the  third 
term  in  eq.  (4.68)  does  not  contribute  appreciably,  as  may  be  verified  by  using  an  iterative  solution 
for  (J”3  to  estimate  the  importance  of  the  term.  At  low  pressures,  the  TQ  contribution  for  counter- 
propagating  waves  remains  small. 

b.  High  pressure  (r?2  > ku,  rj*  > k'u , (I”,  3)'  > |fc  - £k'|u,  r'22  p y2).  In  this  limit  the  sample 
is  completely  thermalized  by  velocity-changing  collisions.  The  probability  that  an  atom  has  not 
undergone  a velocity-changing  collision  in  its  lifetime  is  negligible  and  one  can  easily  deduce  from 
eqs.  (3.32,  4.46,  4.47)  and  (4.63,  4.66,  4.67)  that 


G22(r'  - r)  ~ W0(v)/y2  * <?22(r'  -*  *>), 

(4.72) 

(4.73) 

fM: 


d 
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Moreover,  at  these  high  pressures,  it  follows  from  (4.28a)  that 

G,2(r'  ->  r)  % <5(v  - v')/f12;  G23(v'  - v)  * <5(v  - v')/t23. 


By  substituting  eqs.  (4.72-4.74)  in  (4.27),  one  obtains  the  high  pressure  line  shape 

KA-.  A)  * ,2fyN  \ — + , A?]  (4.75) 

*23r  i2y3y,3 1_  y2  (y13)2  + (A  - A)2J 

There  is  a nearly  constant  background  from  the  SW  chain  plus  a Lorentzian  of  width  y13 
centered  at  A'  = A from  the  TQ  chain.  Velocity-changing  collisions  have  narrowed  the  line  [it  is 
even  narrower  than  the  no-collision  line  width  F,  given  by  eq.  (4.23b)]  and  moved  its  center  from 
A'  * e(k'/k)A  to  A'  = A.  The  line  profile  is  the  same  for  copropagating  or  counterpropagating 
waves  (for  F,y  £ ku,  the  counterpropagating  wave  TQ  chain  begins  to  contribute)  and  is  shown 
schematically  in  fig.  14.  Such  a line  profile  has  yet  to  be  found  experimentally. 

Therefore,  the  line  width  of  the  narrow  resonance  for  copropagating  waves  first  broadens  and 


then  narrows  with  increasing  pressure  while  its  center  shifts  from  A'  = (k'/k) A to  A'  = A.  For 
counterpropagating  waves  the  Lorentzian  contribution  of  the  SW  chain  [first  term  of  eq.  (4.68)] 
centered  at  A'  = —{k'/k) A diminishes  with  increasing  pressure  while  a new  narrow  resonance 
arises  from  TQ  chain  [third  term  in  eq.  (4.68)],  with  the  line  shape  ultimately  reaching  the  form 
(4.75)  at  high  pressures. 


Fig.  14.  High  pressure  limit  of  the  collision  line  shape  when  collisions  are  velocity-changing  in  their  effect  on  p,  The  collision  narrowed 
profile  is  a Lorentzian  of  width  y,  3 centered  at  A'  = A superimposed  on  a broad  background  from  the  SW  chain. 


4.5.  Alternative  level  schemes 

The  results  of  this  section  may  be  taken  over  directly  to  describe  the  alternative  level  schemes  of 
fig.  8. 

For  the  V,  fig  ?b,  A -*  — A,  (A  — kv,)  -*  —(A  — kvt). 

For  the  upward  cascade,  fig.  8c,  A'  -*  - A’,  and  (A'  - ek'v,)  -*  -(A'  - Ek'vx). 
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S.  Inelastic  collisions,  degeneracy,  resonant  collisions 

The  modifications  of  the  theory  needed  to  account  for  inelastic  collisions,  level  degeneracy, 
and  resonant  excitation  exchange  collisions  are  discussed  below.  In  some  cases,  there  is  no  difficulty 
in  extending  the  theory,  in  others,  the  generalization  of  the  theory  is  by  no  means  trivial. 


5.1.  Inelastic  collisions 

Inelastic  collisions  can  be  separated  into  two  categories,  external  and  internal. 

External  inelastic  collisions.  Quite  often,  levels  1 or  2 in  fig.  5 or  levels  1,  2 or  3 in  fig.  8 are 
coliisionally  coupled  to  states  not  shown  in  these  diagrams.  An  external  inelastic  collision  is  one 
that  transfers  population  out  of  the  two  or  three-level  system  under  consideration.  (Inelastic 
collisions  that  transfer  population  into  the  system  can  be  incorporated  into  the  problem  by  an 
appropriate  generalization  of  the  incoherent  pumping  terms  /,(»>),  i = 1,  2,  3.) 

The  effects  of  external  inelastic  collisions  are  easily  calculated.  In  fact,  the  derivation  of 
cpijiv.  f)/f/|con  presented  in  section  2.2  is  unchanged  when  external  inelastic  collisions  are  present. 
Eqs.  (2.20—2.22),  with  elastic  scattering  amplitudes,  are  still  valid.  The  forward  elastic  scattering 
amplitudes  appearing  in  expression  (2.22)  for  ro(»')  implicitly  contain  the  effects  of  inelastic 
collisions,  since  these  elastic  amplitudes  are  related  to  the  total  (elastic  -I-  inelastic)  scattering 
cross  section  via  the  optical  theorem. 


-TOTAL, 


4nh . r , 

(yr)  = — Im[/i(i>r  - vr)]. 


(5.1) 


Since  the  form  of  Tiy(r)  is  unchanged  from  the  purely  elastic  collision  case,  it  is  clear  that  inelastic 
and  elastic  collisions  affect  the  parameter  T0fv)  in  the  same  way.  Thus,  inelastic  collisions  provide 
additional  channels  for  the  decay  of  level  population  densities  />„(»’),  with  a decay  rate  given  by 


f'11(*’)  = 


v)  + Jf  dv_)Vp(v  pK<T?(rr) 


(5.2) 


where  Wjv'  -»  v)  is  the  elastic  collision  kernel  and  <7?(i’)  is  the  inelastic  (or  quenching)  cross  section. 
For  off-diagonal  density  matrix  elements,  inelastic  collisions  add  a contribution 


Jd  rp  Wp( rp)  [<x?(  vr)  + <T?(vr)] 

to  the  real  part  of  T./v).  providing  a broadening  mechanism  in  line  shape  formation.  Since  the 
imaginary  part  of  rjy(r)  is  also  affected  by  inelastic  collisions,  they  give  rise  to  a shift  in  line  profiles. 
Inelastic  collision  broadening  can  often  obscure  the  collisional  narrowing  effect. 

Inelastic  collisions  can  be  studied  by  the  four-level  system  shown  in  fig.  15.  The  pump  field  E 
with  frequency  fi  velocity  selects  a population  p22(r.).  This  population  is  transferred  to  level  3 by 
inelastic  collisions  with  a corresponding  velocity  redistribution  and  the  velocity  profile  p3J(t;.) 
is  probed  by  E as  a function  of  Si'.  This  type  of  experiment  provides  information  on  the  differential 
and  total  inelastic  2-3  scattering  cross  section  and  can  be  used  to  study  the  collision  induced 
relaxation  of  fine  structure  or  rotational  levels.  Note  that  the  atom-probe  field  interaction  is 
negligible  in  the  absence  of  collisions. 


146 


HR  Bernuin , Effects  of  collisions  on  spectroscopic  line  shapes 
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Fig.  15.  A four-level  system  used  to  illustrate  inelastic  collision  effects.  Field  E pumps  the  12  transition  and  field  E'  probes  any  col- 
lisional  transfer  to  levels  3 and  4. 


Internal  inelastic  collisions.  It  may  happen  that  the  energy  separations  of  the  levels  shown  in 
fig.  5 or  8 are  small  enough  to  permit  inelastic  collisions  between  levels  within  the  system  under 
consideration.  For  the  dipole  allowed  1-2  and  2-3  transitions  of  figs.  5 and  8,  inelastic  collisions 
play  a role  only  for  transition  frequencies  < 1012  Hz.  Inelastic  collisions  between  levels  1 and  3 
of  figs.  8a  and  8b  are  likely  to  be  important  if  these  levels  lie  closely  spaced  in  the  same  electronic 
subspace. 

Incorporation  of  internal  inelastic  collision  effects  into  the  quantum  mechanical  transport 
equation  for  dptJ/dt\eM  is  difficult,  in  general,  owing  to  the  collisional  coupling  of  states  with 
different  total  (internal  plus  center-of-mass)  energy.  I do  not  know  of  any  treatment  of  this  problem 
for  arbitrary  level  separations.  However,  if  the  collisionally  coupled  states  have  frequency  separa- 
tions much  less  than  (tc)_i  (tc  = duration  time  of  a collision),  the  impact  approximation  is  still 
valid  and  the  problem  is  greatly  simplified.  In  the  impact  limit,  using  the  techniques  of  section  2, 
it  is  easy  to  derive  [18] 


dpij{v,t)/dt\coU  = - X rfjVlP/J*’-')  + X fdr'W^V  - *)PiJ*\t), 

kn  kn  J 


where 


= 2-^  fdVplFpIVp)  [fk(Vr  - ~ fjn(Vr  - vr)*^ik]', 

J \ * */ 


- -)  is  an  inelastic  “kernel"  given  by  (2.21)  with  the  inelastic  amplitudes  fik(v'  -*■  r)  and 
fjn(v'  -»  v)*  replacing  the  elastic  amplitudes  /|(v'  -»  v)  and  ffv'  -»  v)*,  respectively,  appearing  in  that 
equation. 

Inelastic  collisions  couple  different  density  matrix  elements,  in  contrast  to  elastic  collisions. 
I have  not  yet  analyzed  the  effects  of  internal  elastic  collisions  on  the  line  profiles  derived  in  sec- 
tions 3 and  4.  The  calculation  is  straightforward  to  any  order  of  perturbation  theory  in  the  external 
fields.  There  may  be  some  new  line  shape  features  introduced  by  the  collision  induced  coupling 
of  different  density  matrix  elements. 

The  four  level  system  of  fig.  15  may  also  be  used  to  probe  inelastic  collision  effects.  Imagine 
that  levels  1 and  2 are  collisionally  coupled  to  levels  4 and  3,  respectively.  The  field  E'  can  monitor 
population  transfer  from  levels  1 to  4 and  2 to  3 as  well  as  any  transfer  of  “coherence”  from  p12(r') 
to  Pj^v)  [19]. 
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5.2.  Degeneracy 

Typically,  each  of  the  levels  of  the  “two"  or  “three”-level  systems  contains  several  degenerate 
magnetic  substates,  leading  to  a corresponding  increase  in  the  number  of  states  in  the  system. 
Inelastic  collisions  between  magnetic  substates  of  a given  level  are  energetically  allowed  so  that 
internal  inelastic  collisions  can  occur.  For  all  but  very  strong  external  magnetic  fields,  the  magnetic 
sublevel  spacing  is  small  enough  for  the  impact  approximation  to  be  valid  and  the  collisional  time 
evoluation  of  density  matrix  elements  is  given  by  eq.  (5.3). 

Collisions  result  in  the  relaxation  of  both  the  populations  and  “coherences”  (off-diagonal  density 
matrix  elements)  of  the  magnetic  substates  within  a given  level.  Transfers  of  the  type 
Pim,2m ■ are  also  possible,  in  principle.  The  various  “rates”  and  “kernels” 

are  subject  to  certain  selection  rules  and  it  is  often  useful  to  expand  the  density  matrix 
in  terms  of  its  irreducible  tensor  components  in  order  to  exploit  any  symmetries  in  the  problem. 

Information  on  the  differential  as  well  as  total  cross  sections  associated  with  magnetic 
relaxation  can  be  obtained  by  the  saturation  spectroscopy  scheme  indicated  in  fig.  15,  in  which 
levels  2 and  3 represent  different  magnetic  substates  of  a given  level.  By  choosing  the  polarization 
of  fields  E and  E'  properly,  one  is  able  to  determine  the  rate  at  which  collisions  transfer  population 
from  level  2 to  3 and  the  associated  change  of  velocity  occurring  in  those  collisions.  It  would  be 
interesting  to  perform  an  analogous  experiment  to  study  the  velocity  changes  associated  with 
transfers  of  off-diagonal  density  matrix  elements  pim%im(v)  -»  pim (r)  (m  # m',  m"  / m'"),  since 
there  is  as  yet  no  simple  physical  picture  to  describe  this  process. 

5.3.  Resonant  collisions 

For  the  active  atom  and  perturber  level  schemes  of  fig.  16  resonant  exchange  collisions  may 
occur  in  which  the  perturber  and  active  atom  exchange  roles.  It  is  possible  to  transfer  level  popula- 
tion as  well  as  off-diagonal  density  matrix  elements  from  one  atom  to  another. 


Fig.  16.  Active  atom  and  perturber  energy  level  schemes  in  which  resonant  excitation  exchange  is  possible.  Collisions  can  transfer 
population  or  “coherence  " (off-diagonal  density  matrix  elements)  from  an  active  atom  to  a perturber. 


The  transfer  of  population  is  easy  to  understand.  A pump  field  produces  some  velocity  selected 
population  in  an  active  atom.  Resonant  exchange  collisions  can  transfer  this  population  to  another 
atom.  However,  since  the  atom  receiving  the  population  transfer  has  an  arbitrary  velocity,  the 
velocity  selectivity  created  by  the  pump  field  is  lost  in  the  transfer  process.  In  effect,  each  transfer 
collision  completely  thermalizes  the  active  atom  velocity  distribution.  Thus,  in  addition  to  the 
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probe  absorption  line  shape  features  shown  in  fig.  13c  arising  from  absorption  from  active  atoms 
that  have  not  transferred  population  to  other  atoms,  there  will  be  an  additional  broad  (Doppler 
width)  background  signal  owing  to  absorption  from  velocity  thermalized  atoms  that  have  received 
population  via  excitation  exchange  collisions. 

The  transfer  of  off-diagonal  density  matrix  elements  is  more  difficult  to  analyze,  since  one  is 
again  faced  with  the  problem  of  associating  a velocity  with  off-diagonal  density  matrix  elements. 
While  I have  not  analyzed  this  process  in  any  detail,  1 imagine  that  off-diagonal  density  matrix 
elements  can  not  be  transferred  in  strong  collisions,  owing  to  both  the  large  phase  shifts  and  quan- 
tum interference  effects  that  characterize  such  collisions.  There  may  be  some  possibility  to  observe 
this  type  of  exchange  in  weak  collisions.  It  is  an  interesting  effect  since  phase  information  is  trans- 
ferred in  the  collision. 

There  is  a complication  that  very  often  obscures  resonant  excitation  exchange.  For  levels 
radiatively  connected  to  the  ground  state,  the  exchange  of  real  photons  (radiation  trapping)  is 
possible.  For  example,  an  active  atom  may  spontaneously  emit  a photon  that  is  absorbed  by 
another  atom  in  the  sample.  This  process  gives  rise  to  the  same  type  of  line  shape  modifications 
as  does  resonant  excitation  exchange.  In  practice,  it  is  often  impossible  to  distinguish  the  two 
effects,  since  the  resonant  excitation  and  radiative  trapping  exchange  rates  are  both  proportional 
to  the  active  atom  density  over  a broad  pressure  range. 


6.  Conclusion 

I have  described  in  general  terms  the  way  in  which  spectroscopic  line  shapes  reflect  the  collisional 
processes  occurring  in  atomic  and  molecular  systems.  Additional  experiments  are  required  to  test 
the  potential  of  these  techniques  and  to  determine  whether  or  not  linear  and  non-linear  spectro- 
scopy can  provide  a practical  method  of  probing  atomic  and  molecular  collisions. 
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with  fundamental  problems  of  light-matter 
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applications  to  spectroscopy. 
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This  volume  contains  the  25  invited  lec- 
tures and  seminars  presented  at  the 
NATO  Advanced  Study  Course  on  Co- 
herent Optical  Engineering  (the  sixth 
course  of  the  International  School  of 
Quantum  Electronics,  affiliated  with  the 
Centre  for  Scientific  Culture  “E.  Majorana”, 
Erice,  Sicily).  The  aim  of  the  course  was 
to  provide  an  authoritative  expose  of 
several  optical  measurement  methods, 
including  holographic  and  speckle  inter- 
ferometry, laser  interferometry,  laser 
telemetry,  intensity  correlation  spectro- 
scopy and  Doppler  velocimetry,  white 
light  and  “moirG"  techniques,  and  to 
describe  the  state-of-the-art  in  the  thriving 
field  of  optical  information  processing. 
The  course  was  made  quite  self-contained 
by  this  series  of  lectures  which  eluci- 
dates the  basic  principles  underlying 
coherent  optical  engineering  and  re- 
views the  design  of  and  the  work  done 
using  the  most  important  optical  devices, 
such  as  laser  sources,  laser  deflectors 
and  modulators,  detectors  and  recording 
media.  Since  it  gives  the  most  up-to-date 
available  review  of  developments  in 
Coherent  Optical  Engineering,  this  book 
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Ellipsometry  is  a unique  optical  technique  of  great  sensitivity  for 
in  situ  non-destructive  characterization  of  surface  (inter-facial) 
phenomena  (reactions)  utilizing  the  change  in  the  state  of  polar- 
ization of  a light-wave  probe.  Although  known  for  almost  a century, 
the  use  of  ellipsometry  has  increased  rapidly  in  the  last  two  decades. 
Among  the  most  significant  recent  developments  are  new  applica- 
tions, novel  and  automated  instrumentation  and  techniques  for 
error-free  data  analysis. 

This  book  provides  the  necessary  analytical  and  experimental  tools 
needed  for  competent  understanding  and  use  of  these  develop- 
ments. It  is  directed  to  those  who  are  already  working  in  the  field 
and,  more  importantly,  to  the  newcomer  who  would  otherwise  have 
to  sift  through  several  hundred  published  papers.  The  authors  first 
present  a comprehensive  study  of  the  different  mathematical  repre- 
sentations of  polarized  light  and  how  such  light  is  processed  by 
optical  systems,  going  on  to  show  how  these  tools  are  applied  to  the 
analysis  of  ellipsometer  systems.  To  relate  ellipsometric  measure- 
ments to  surface  properties,  use  is  then  made  of  electromagnetic 
theory.  Experimental  techniques  and  apparatus  are  decribed  and 
the  many  interesting  applications  of  ellipsometry  to  surface  and 
thin-film  phenomena  are  reviewed. 

This  reference  work  is  addressed  to  researchers  and  students  with  a 
strong  interest  in  surface  and  thin-film  physics  and  optics  and  their 
applications.  The  book  should  be  on  the  shelf  of  any  library  in  the 
fields  of  solid  state  physics,  physical  chemistry,  electro-chemistry, 
metallurgy  and  optical  engineering. 
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la  this  work,  a theory  of  bier  nonlinear  spectroscopy  and  its  application  to  the  study  of  atomic  collisions  is 
developed.  A three-level  system  (TLS)  is  considered  in  low-pressure  gases.  The  effect  of  weak  velocity- 
chancing  collisions  (vcc)  upon  atomic  coherences  is  introduced  into  the  TLS  equations  established  by  Hansch 
and  Toschek.  It  appears  that  a TLS  is  suitable  for  the  observation  of  the  effect  of  vcc  upon  optical 
coherences. 


In  a three-level  system,  it  is  possible  to  have 
a saturated  absorption  signal  arising  solely  from 
terms  related  to  optical  coherences  and  not  de- 
pendent on  populations.  By  examining  the  effect 
of  collisions  on  these  profiles,  one  eliminates  any 
problems  arising  from  velocity- changing  collisions 
on  populations  which  would  tend  to  mask  their  ef- 
fects on  "coherences.” 

There  exists  a large  literature  of  calculations  of 
velocity-changing  collision  effects  in  both  linear 
and  nonlinear  spectroscopy.6'10  To  our  knowledge, 
no  calculation  exists  in  which  the  model  of  weak 
velocity-changing  collisions  affecting  internal  co- 
herences in  a three- level  system  (TLS)  has  been 
fully  explored . In  this  work,  such  a model  is 
adopted  and  may  prove  relevant  to  explain  recent 
experimental  results  involving  saturation  spectro- 
scopy in  a three-level  system.11 


The  fruitful  development  of  nonlinear  laser  spec- 
troscopic techniques  during  the  last  few  years  has 
provided  new  evidence  for  collision- induced  homo- 
geneous broadening,  shifts,  and  even  distortion 
of  spectral  lines  in  low-pressure  gases.  Saturated 
absorption1*4  and  two- photon6  spectroscopy  have 
been  used  for  such  studies . 

The  use  of  nonlinear  spectroscopy  to  study  col- 
lisional  effects  represents  a marked  improvement 
over  linear  spectroscopy  in  measuring  the  effects 
of  velocity- changing  collisions  (vcc).  In  linear 
spectroscopy,  one  starts  with  a thermalized  sam- 
ple. Consequently,  elastic  collisions  do  not  alter 
the  velocity  distribution  of  the  various  level  pop- 
ulations of  the  atom,  which  are  already  in  equili- 
brium. Collisions  will  affect  only  the  coherence 
or  off-diagonal  density-matrix  elements,  leading 
to  a disturbance  of  the  phase  of  the  atomic  oscil- 
lators. In  general  this  phase  disturbance  consists, 
of  an  inseparable  combination  of  velocity-changing 
and  phase- interrupting  effects,  although,  in  cer- 
tain limits  (see  below),  one  contribution  can  dom- 
inate. In  linear  spectroscopy,  effects  of  vcc  are 
easily  lost  in  the  large  widths  of  spectral  profiles 
arising  from  the  Doppler  effect,  in  some  cases, 
the  effect  of  vcc  on  atomic  coherences  can  be  de- 
tected by  a narrowing  of  the  Doppler  profile,  but, 
in  general,  vcc  are  difficult  to  detect  in  linear 
spectroscopy. 

Nonlinear  Doppler-free  spectroscopy  provides 
more  promise  for  studying  subtle  collisions)  pro- 
cesses. In  such  experiments  one  excites  a given 
velocity  group  of  atoms.  By  probing  the  system, 
one  can  determine  the  rethermalization  of  this 
population  velocity  group  resulting  from  collisions. 
Moreover,  there  are  systems  where  one  can  at- 
tempt to  study  the  effect  of  velocity- changing  col- 
lisions on  the  atomic  coherences. 


II.  TLS  UNIDIRECTIONAL  SPECTROSCOPY 

In  a typical  experiment,  two  cw  monochromatic 
collinear  laser  beams  interact  with  an  atomic  or 
molecular  gas  contained  in  a low-pressure  cell 
(Fig.  1).  The  field  E(tii,k)  is  resonant  with  and 
saturates  the  1-2  transition  alone,  selectively  ex- 
citing optical  dipoles  with  an  axial  velocity  close 
to  zero.  The  field  £'(<ii',Jfr')  with  a frequency  <*?' 
close  to  that  of  the  2-3  transition  monitors  the 
changes  in  the  velocity  distribution  and  phase  co- 
herences of  the  system,  induced  by  the  pump  field 
in  the  presence  or  absence  of  collisions. 

The  parameter  measured  in  typical  experiments11 
is  that  part  of  the  absorption  coefficient  of  the 
probe  beam  containing  saturation  effects;  the  most 
interesting  situation  occurs  with  unidirectional 
beams  and 

A theoretical  description  of  this  system  was 
given  by  H&nsch  and  Toschek.1*  Using  perturba- 
tion theory,  neglecting  collisions,  and  assuming 
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Colliskmally  aided  redistribution  of  scattered  laser  light  is  suggested  as  a method  to  cool  or  heat  gaseous  samples.  The 
efficiency  is  evaluated  and  restricting  conditions  are  considered.  Some  potential  applications  are  given. 
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Many  recent  treatments  oflaser  light  scattering 
(e.g.  refs.  [1  — S J ) imply  that  the  process  can  be  uti- 
lized to  heat  or  cool  an  atomic  sample.  To  our  knowl- 
edge, there  has  been  no  explicit  discussion  of  the 
phenomenon,  hence  we  shall  try  to  explain  the  mech- 
anism and  estimate  its  potentialities. 

The  effect  is  easily  explained  with  reference  to 
fig.l  . Weak  monochromatic  light  of  frequency  fi 
acts  on  a two-level  atomic  system  with  a stable  lower 
level  (i.e.  a zero  width  eigenstate).  For  the  prelimi- 
nary' discussion  the  atoms  are  assumed  stationary, 
and  no  Doppler  shifts  occur.  Without  collisions,  ener- 
gy conservation  requires,  that  any  fluorescence  (or 
Rayleigh  scattering  if  you  prefer)  has  to  occur  at  the 
frequency  fl  [6] . The  scattering  process  first  absorbs 
energy  to  a specific  position  inside  the  frequency 
band  of  the  natural  line  width,  and  rcemission  must 
start  from  this  point;  the  result  is  the  delta  function 
spectrum  of  fig.  lb. 

With  collision  events  present,  the  atoms  can  ex- 
change energy  during  excitation  and  fluorescence  ac- 
quires a new  frequency  component  centered  at  the 
atomic  level  spacing  (1-5)  u>  as  shown  in  fig.  Ic. 
Collisions  provide  a mechanism  to  redistribute  the 
excited  state  population  over  the  natural  width  of 

* Supported  by  the  U.S.  Office  of  Naval  Research  under  con- 
tract number  NOOOI4-77-C-0533. 


the  upper  level,  y say.  This  effect  we  shall  call  colli- 
sionally  aided  fluorescence  (CAF).  Since  each  scatter- 
ing event  transfers  a photon  from  the  frequency  £2 
to  another  frequency,  the  energy  difference  must  be 
compensated  by  the  translational  energy  of  the  col- 
liding atoms,  and  we  have  to  consider  moving  atoms, 
with  the  Doppler  width 


T\’i 


1<C-U 


Fig.  I . (a)  Two-level  system,  with  the  energy  difference  Tiw, 
coupled  by  the  field  at  frequency  n.  The  lower  level  is  sharp, 
the  upper  one  has  the  natural  width  y.  (b)  Without  collisions 
energy  conservation  requires  the  scattered  light  to  have  the  fre- 
quency, n'ow*  An,  centered  at  n.  (c) Collision  event  pro- 
vide Fourier  components  in  the  range  r "• , allowing  the  popu- 
lation to  settle  on  the  excited  level  before  decay.  Hence  there 
appears  a component  centered  at  O’  = w in  the  scattered  ligh  t. 
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«0~(n/c)(*Br/A/),/2. 

Tuning  with  fl  > to  produces  heating  and  that 
With  n<«  cooling  of  tire  gas.  The  magnitude  of  the 
effect  depends  strongly  on  the  ratio  of  the  detuning 
A * n - u)  to  the  Doppler  width  cuD,  the  total  ener- 
gy  per  particle  WT  « *„rand  the  duration  of  a colli- 
sipn  r( ; the  Rabi  frequency  x = v£/2h , the  upper 
state  lifetime  and  the  collision  rate  r.  Here  p is 
the  transition  dipole  and  E is  the  field  strength.  We 
Pnly  need  to  consider  phase-interrupting  collision 
effects.  We  want  to  choose  |A|  large  enough  to  have 
a considerable  energy  change  per  CAF  event,  but  too 
large  detunings  will  produce  no  upper  state  popula- 
tion to  act  as  the  source  for  the  fluorescence.  For 
common  values  of  the  parameters  of  a gas  sample, 
we  can  choose  the  tuning  such  that 

hy*>Tc,>1A|>wD>r*  {1) 

which  also  provides  the  advantage  that  collisions  can 
be  treated  by  the  impact  approximation  as  |A|  rc  < 1 . 

Each  collision  event  which  gives  rise  to  CAF 
changes  the  translation  energy  of  the  system  by  an 
amount  /i|A|.  If  we  want  spontaneous  emission  to 
dominate  we  must  assume  a small  population  prob- 
ability for  the  upper  state,  n2  < 1 . The  rate  of  CAF 
events  is  then  given  by  |4] 

"2  " “ 2X2r/A2.  (2) 

* +(r+7/2)2 

Hence  the  total  rate  is  the  product  of  a collisional 
rate  and  the  excitation  probability.  The  population 
of  the  upper  level  decays  at  the  rate  y and  in  steady- 
state  we  have 

nj  m 2x2r/yA2.  (3) 

Requiring  this  to  be  much  smaller  than  one  we  find 
the  condition 

2x2r/A2«v  (4) 

This  also  implies  that  each  excitation  event  will  im- 
mediately be  followed  by  spontaneous  emission, 
which  is  necessary  for  optimum  energy  exchange. 

The  fractional  rate  of  energy  exchange  per  active 
atom  is 

dHydf  //>a\  2v2r 

T W't"*  <5) 


Assume  we  have  a sample  of  jV  active  atoms  and  N' 
perturbers  per  unit  volume,  and  the  laser  illuminates 
a volume  l'T,  which  is  less  than  the  total  volume  V. 
Then  the  fractional  energy  (temperature)  change  per 
unit  time  is 

, ^ d log  r / \/d,*/T/d/\ 

;T"  dr  \(Ar  + iV’)K/\  WT  ) 

T.2X2r/tA 

a2  wtIn  + n')v 

Eq.  (6)  is  valid  only  when  cq.  (4)  is  satisfied.  With 
typical  values  Wy/h  = lO^s-1,?”1  = 10,2s_1  and 
«D  = 10>°  s-1 , we  can  choose  |A?  = 101 1 s-1  satis- 
fying eq.  (1).  With  moderate  power  requirements  (e.g. 
W/cm2)  one  can  achieve  values  of  x2/A2  — 10-5.  For 
simplicity  we  can  take  the  collision  rate  to  be 
(Af+jV  )uro  — {N+N')  10~ 10  cm2  s-1 , where  u,  is  the 
relative  perturber-active  atom  velocity  and  o is  the 
collision  cross  section.  Setting  FL  = V we  obtain 
from  eq.  (6)  an  attainable  value  of 

/T  = 10_,7yV  (cm3*-1)  (7) 

subject  to  the  restriction  (eq.  4) 

10-,s(/v+AO«7  (s  cm-3).  (8) 

An  increase  of  foreign  gas  perturbers  is  seen  to 
leave  the  effect  unchanged,  since  /T  depends  on  the 
active  atoms  only.  The  enhancement  of  the  rate  of 
CAF  events  is  counteracted  by  the  growth  of  the  to- 
tal energy  reservoir  to  be  cooled  or  heated.  Total 
compensation  occurs. 

Eq.  (7)  suggests,  that  potentially  large  heating  or 
cooling  effects  may  occur.  At  1 .0  torr  the  time  scale 
in  which  the  energy  of  the  gas  is  removed  or  doubled 
is  1.0  s(sic!).  The  condition  (8)  can,  however,  be  quite 
restrictive  when  radiation  trapping  occurs.  Then  y 
approaches  zero  for  densities  of  the  order  of  1013 
cm-3,  giving  an  upper  limit  of /T  = 10~4  s-1.  Even 
in  this  case  heating  or  cooling  is  bound  to  occur  for 
long  enough  interaction  times.  To  enhance  the  heat- 
ing or  cooling  one  should  try  to  avoid  radiation  trap- 
ping. If  the  upper  state  has  a cascade  decay  branch 
back  to  the  ground  state  no  problem  will  arise.  It  may 
be  possible  to  depopulate  the  upper  level  to  some 
high  lying  state  by  optical  pumping,  and  the  atom 
would  cascade  back  to  the  ground  state  by  different 
routes.  If,  however,  we  empty  the  lower  level,  each 
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